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Preface 



Although Charles Babbage laid down several ideas for computing "en- 
gines," the forerunners of today's computers were largely developed in the 
early 40's as part of the war effort. The Robinson series cryptoanalytic ma- 
chines developed in England in 1941 spawned many families of computers 
still in use today. The MIT differential analyzer and real-time aircraft simula- 
tion project led to the Whirlwind, and eventually to the immensely successful 
DEC family of PDP computers (Programmed Data Processors). And, of 
course, Eniac built at the University of Pennsylvania, was the guiding light 
behind Univac, IBM and many other successful manufacturers. 

Many people today poke fun at these early machines and regard them as 
dinosaur-like relics. However, it is interesting to consider that a large-scale 
computer of about 25 to 30 years ago had about the same amount of power as 
a typical personal computer of today. It was generally not as reliable or user- 
friendly as a personal computer, and, of course, cost tens of thousands times 
as much. Why bring this up? 

Because a computer of the 50's required that the programmer be very 
clever and resourceful to solve problems within the capabilities of the com- 
puter. He did not have vast gobs of memory available, blinding quick calcula- 
tion speed, or random disk access. In other words, he had about the same 
problem to face as you do with your personal computer. 

I do not mean to imply that your personal computer is not a full-fledged 
computer. It certainly is just as much a computer as a room-filling giant of to- 
day. However, because of the relatively small memory, it cannot store a large 
data base. Nor is it suitable for extensive word processing or massive calcula- 
tions. Highly detailed graphics are best left to other machines as well. 

What can we learn from the computing pioneers of the 50's that will help 
us today? Perhaps most important is the discipHne of thoroughly analyzing a 
problem, breaking it down into manageable steps, and solving it a step at a 
time. It is also important to determine what can be done "off line" and what 
must be done on the computer. 
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That is what this book is all about. While it has more than 50 ready-to-run 
programs, the main thing you should look for from the book is an approach 
to solving problems — ^big and small. Some of the problems demonstrate the 
capability of the computer; others identify its shortcomings. It is important to 
be famihar with both the strengths and weaknesses of your tools so you can 
recognize the types of jobs for which they are suitable (and not suitable). 

This book focuses primarily on mathematical and educational applications 
for the computer. There are many other excellent sources of information 
about other applications and for making best use of your personal computer. 
Using the approaches described in this book should enable you to easily con- 
vert programs and use applications from other books and magazines such as 
Creative Computing, 

This book is designed to be read with a working computer at hand. While 
there is textual material to be read, the most important things are the experi- 
ments and problems to be tried with your own computer. The book raises 
many questions for which you should try to find answers. There are no an- 
swers to these questions and problems in the back of the book; you should be 
able to discover the answers as you work the problems out on your computer. 

You will be able to incorporate many of the routines and approaches in the 
book into programs of your own as you use your computer to deal with "real 
world" problems. Other programs will simply point you in the right direc- 
tion. And some of the programs in the book are just plain fun. Learn. Experi- 
ment. Have fun! 

Morristown, New Jersey 

March 1983 David H. Ahl 
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Drill and Practice 



Throughout Hfe, there are certain things that simply must be memorized. 
Obvious things that fall into this category are the addition and multiplication 
tables, the spelling and meaning of words, how to tell time, and the monetary 
system. 

However, depending upon one's chosen profession, there are many other 
things to be memorized. A doctor must know what diseases match what 
symptoms. A chemist must know the gas laws, the properties of elements and 
so on. A pilot must instantaneously know the meaning of readings on scores 
of instruments. 

To memorize a set of facts, you must go over them again and again and 
keep trying different variations. Here is where the computer comes in. It is 
able to present randomly scores of different problems to you for as long as 
you wish. Some programs will automatically adjust to your level of com- 
petence and will grade you; other programs simply present the problems and 
leave the grading up to you. 

There are four programs in this chapter and two in the Science chapter 
which present material in a drill and practice format. Examine the methods 
used in these programs and then make up some drill programs of your own 
for subjects with which you are having trouble, or make up programs for 
other members of your family. 
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Addition Practice 



This program demonstrates a simplified addition drill and practice routine. 
This type of drill is sometimes called computer assisted instruction (CAI), al- 
though CAI can also apply to tutorial and other approaches as well. 

When the program is run, it will first ask "No. of digits?" You enter a 
number and each addend will contain that many or fewer digits. 

The program will present any number of practice problems that you spec- 
ify. The program presents each problem in turn. The program will not pro- 
ceed to the next problem until the current one has been answered correctly. 
After the last problem is answered correctly, the score is printed with an 
appropriate comment. 

There are many improvements and extensions possible in a program like 
this one. For example, you might want to modify the program so it tells the 
user the correct answer after a problem has been answered incorrectly two (or 
three) times. 

A more complicated modification would be to change the program to 
present different kinds of arithmetic problems such as subtraction, multiplica- 
tion, and division. 

Some of these modifications have been made in the next program. 
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190 
200 
210 
220 
230 
240 
250 
N E " 

260 GDTQ.SO 
270 PRINT : 
ORRECT THE" 
28 PRINT ^ „ 
290 IF U<.22*P 
300 PRINT "BUT 
3 10 STOP 
320 PRINT "GnOD 
33 STOP 



THEN 2 10 



TRY flGFlIN" 



IF Q>i 
W = W+1 

PRINT "UHflT? 
GOTO 100 
PRINT "RIGHT!"? 
IF f^>=rp THEN 270 
PRINT "HERE IS fiNDTHER 



"YOU GOT " J P-W 5 " C 
"FIRST TIME. " 



THEN 320 

YDU MISSED" ;W 

WORK ! " 



> R ij N 
fl II D IT I □ N 



PRACTICE 



NO. 
NO . 



OF DI6ITS^=3 „ 
□F PROBLEMS =5 



73 1 
761 



? 1492 
RIGHT ! HERE 

387 
+ 2S 



? 415 
R I 6HT ! HERE 



IS flNQTHER ONE. 



ANOTHER DNE- 



1 33 
4 08 



? 54 1 
RIGHT f HERE 

933 
+ 987 



? 1 S 1 
WHAT? TRY 



IS ANOTHER ONE. 



AGAIN 



93 3 
987 



? 1920 
R I GHT ! HERE 

229 
+ 103 



IS ANOTHER ONE. 



? 3 3 2 
R IGHT ! 

YDU GOT 4 CORRECT THE 
FIRST TIME. 

GOOD WORK! 
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Addition Practice, 
Adjusted by Grade Level 



One of the major disadvantages with many drill and practice exercises is 
that they tend to be either boring or frustrating, depending upon the ability of 
the user relative to the level of the material. To compensate for this, a method 
is needed which will adjust the difficulty of the problems to the ability of the 
user. 

Ideally, such a system would weigh the most recent performance most 
heavily but would not ignore previous performance. It should allow a user to 
advance to more difficult problems than his current mastery level. It should 
also continue to give some practice on problems already mastered. 

Some commercial software packages approach these goals along traditional 
lines, i.e., determine in which type of problems a student should receive prac- 
tice by using a complicated computer managed instruction score recording 
and adjustment system. 

The approach here is more innovative; it uses a single measure for each 
type of problem — call it "estimated grade level" — which meets all of the 
objectives stated above. 

How does it work? The most recent problem presented counts 10% of the 
overall score if it was answered correctly and was over the current user grade 
level, or if it was answered incorrectly and was under the current user grade 
level. Otherwise it is ignored. This may be easier to visualize in the form of a 
chart: 

Answer 
Right Wrong 



Higher than Raise student Ignore 
grade level grade level 



Problem 

grade level ^ grade level 



Lower than Ignore Lower student 



At first glance this might look complex and somewhat goofy, however, 
what it really means is that a student is rewarded for doing a problem beyond 
his grade level but he is not penalized if he cannot do it. On the other hand he 
is penalized if he cannot do a problem lower than his grade level, but is not 
rewarded for doing one lower. 

Each problem affects the estimated grade level a little bit, with the most 
recent problems being weighed the most heavily. If the current grade level of 
a student is L and the level of the most recent problem to be averaged in is P, 
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then the averaging formula is simply: 

L = .9L + .IP 

The remaining task before a program can be written is to assign a grade 
level to each problem presented. Unfortunately, this will vary depending 
upon the local school system, the textbook used, and the teaching method. 
Also a huge data base can not be stored in a small computer, so it is desirable 
to devise a simple method of determining grade level for different problems. 
One straightforward approach is to present problems up to one-half a grade 
level over and under where the student currently is. Thus the overall range of 
problems for a student at grade level 3.2 would be 2.7 to 3.7. 

How do we generate the right problems? Consider one type of skill, vertical 
addition. It is normally introduced in the first grade and continues through 
Grade 4 (actually 4.9). The simplest problem in this program is 1 + 1 and the 
most difficult is 999 + 999. Since learning is not a linear process (it is slow at 
first, and then progresses rapidly), an exponential formula can be used. For 
example: 

Addend = 1.73 x (Grade level)^ 



or 



Grade level = V Addend/1.73 

This means that students at various grade levels will be working with the 
following maximum addends: 

Grade level Addend 
1.0 1 
2.0 27 
3.0 140 
4.0 442 
4.9 997 
Now it is a relatively straightforward, although somewhat tedious, matter 
to tie all these elements together in a computer program. 

A few notes about the program. The variable G2 is the problem grade level 
that is always within one-half of a grade level of the current student level, Gl. 
The complicated mess in Statement 340 produces a moving grade level 
average. 

The recording of the grade level and carrying it over to the next lesson is a 
manual process. On a computer system with a permanent mass storage de- 
vice, this would be kept on the system. 

There are many possible changes and extensions to this program. For 
example it could present different types of problems such as horizontal addi- 
tion, vertical and horizontal subtraction, as well as multiplication, division 
and fraction problems. 
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20 5-0"°°"^^^ 

^iglg fls""^^' ^MPUT 

4 " P R I N T "YOUR Fl N S Ul E R . " 

^?,= f:§P,n "IJHRT IS YDUR GRRDE 

6 INPUT Gl 

70 62=G1-.5*RND 

ag R=INT<2 + 1 . 73*G2 *-4> 

yp Fl = I NT < 1 00*G2*RNIi ) 

!9P if'r'^'^'^ THEN 90 
110 B = R - fl 

115 PRINT 

l^^O PRINT TfiBC8-LENCSTRJ<B> > 

ll?fi5?J^S " "'TflE<8-LEN<ST 

140 PRINT " " 

150 R=R+B 

160 PRINT Tfli!<:7-LEN(STRS<R) > 

170 IMPUT G 

200 IF G=SSS9 THEN 500 

210 IF G=R THEN 310 

2 2 U = W + 1 

2 3 IF W > 1 THEN 260 

Ui %^vi\ir ''''''' ^"^^^ 

260 PRINT ''YDU MISSED THAT □ 
N E T W I 1.. E . 

270 PRINT "THE CORRECT flNSWE 

K IS ? R 
28 I.J = 

ii!S IStB^I^J 

3 10 W =5 

3 2 PRINT " C □ R R E C T I " 
'd30 IF G2<gl THEN 400 
3 4 g ^ G 1 = . 9 * G 1 + . 1 * S Q R S Q R C R . - 2 .'•• 
l « 7 ^1 ; 
400 PRINT 

410 PRINT "HERE IS RNDTHER.. 

420 GOTO 70 
5 PRINT 

"'tJjS..'^?^'^'^ ^'OKRY. SO LONG FOR 

gl = < INT< 100*G1 > >^ 100 
?|l5^J/'YaUR NEW GRADE LE 

V t. L 1 o .1 b 1 

999 END 
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HI. TO STDP^ INPUT 99 9 9 flS 
V □ U R R N S W E R . 

WHRT IS YDUR GRfllE LEVEL? 3 
5 

+ 87 



? 92 
CORRECT ! 



HERE IS RNOTHER. 
58 

+ 67 

? 1 1 4~ 
WRONG. TRV RGRIN 

58 

+ 67 



? 125 
CORRECT ! 



HERE IS RNOTHER 
4 

+ 104 



? 10 8 
CORRECT ! 



HERE IS ANOTHER 
28 

+ 50 



? 99 9 9 

□KRY. SO LONG FDR NOW. 
YOUR NEW GRRUE LEVEL IS 

* * HONE «• * 
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Time/Speed/Distance Problems 



As well as being able to present simple numerical drill and practice prob- 
lems, the computer can present word problems for solution as well. 

In this program the formula relating time, speed and distance was appUed 
to a problem involving both a car and train. 

The problem can be stated as follows. A car traveling C miles per hour 
(computer generates an integer 40 through 65) can make a certain trip in D 
hours (computer generates an integer 5 through 20) less than a train traveling 
at T mph (computer generates an integer 20 through 39). How long does the 
trip take by car? When the two simultaneous equations are solved they pro- 
duce the single equation for the answer shown in Line 110. 

Notice the calculation in Line 120. This calculates the percent difference 
between the actual answer and the one entered by the user. 

Notice also that the computer calculates the correct answer in Line 1 10 and 
prints it (on the screen) in Line 200. This answer may have many decimal 
places; as the program is written it is rounded off to two decimal places. If the 
0.5 was not added in Line 110, the number would be truncated and not 
rounded off. This is an important calculation and one which you will find in 
many other programs throughout the book. 

Consider other problems that can be used as the basis for this kind of drill 
and practice exercise. Teachers, for example, might wish to have students 
write drill and practice programs on their own. Different problems could be 
given to one or a small group of students to serve as the basis for a program. 

After the programs are written, students can try out the programs of other 
class members. This approach ensures that each student not only understands 
the type of problem assigned to him, but also gets practice in solving other 
problem types as well. This is an effective technique for stimulating interest as 
well as for learning how to solve word problems. 
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5 RRNDDMIZE 

10 C = I N T c: R t J D * 2 5 > + 4 G 

2 D = I N T -.: R D * 1 5 > + 5 

3 T = I N T R r^l D * 1 9 > + 2 

4 PRINT ^'R CAR TRAVELING" 5 C 

5 " MPH CRN " 

5 PRI^4T ''MRKE R CERTAIN TRI 
P I N " ? If 

6 PRINT "HOURS LESS THRN R 
TRAIN" 

7 PRINT "TRAVELING flT"liT5"l1 
P H " 

80 'print "HDW long DDES THE 

TRIP TAKE" 

90 PRINT "BY CAR"? 

100 INPUT R 

110 V = n-J T ( . 5 + 1 * < D * T / C C - T ::0 

> 1 

12 E = I r J T < R B S < V - R > * 1 R > + . 

5 > 

130 IF E>5 THEN 170 
140 PRINT "GOOD. YOUR ANSWER 
W AS " 

150 PRINT "WITHIN" 5 E; "PER C EN 
160 GOTO 190 

17 P R n4 T "SORRY. YOU WERE 
FF BY'' 

180 PRINT E5 ^'PERCENT. 
190.^PRINT "THE CORRECT RNSWE 

2 '"^ P R I N T V 5 " H □ U R S . " 
2 10 PRINT 

220 PRINT "HERE IS ANOTHER.. 
It 

230 GOTO 10 
9 9 9 END 



R U N 

R CRR TRRVELING 41 MPH CRN 
MAKE A CERTAIN TRIP IN 6 
HOURS LESS THRN R T R R I N 
TRRVELING RT 33 MPH. 
HDW LONG DOES THE TRIP TRKE 
BY CAR? 25 

GOOD. YOUR ANSWER URS 
W I T H I 1 PERCENT, 
THE CORRECT ANSWER IS 
.24. 75 H □ U S . 

HERE IS ANOTHER. . . 
A CAR TRAVELING 42 MPH CAN 
MAKE A CERTAIN TRIP IN 7 
HOURS LESS THAN A TRAIN 
TRAVELING AT 24 MPH. 
HDW LONG DDES THE TRIP TAKE 
B Y C R R ? '9 . 3 
GOOD. V D U R ANSWER WAS 
0./ 1 THIN PERCENT. 
THE CORRECT ANSWER IS 
9.33 HOURS. 
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Kinematics Problems 



Kinematics relates to the dynamics of motion of bodies apart from consid- 
erations of mass and force. Like many other types of problems, these can be 
generated and presented by the computer to provide practice in solving them. 

This program presents a simple kinematics problem for solution. The com- 
puter generates a new value for each problem. The problem is as follows. A 
ball (or any other object) is thrown up at velocity V meters per second (com- 
puter generates an integer between 5 and 40). The user must then calculate 
three factors about the resulting flight of the ball: maximum height, time until 
it returns, and velocity after T seconds (computer generates a time less than 
the total flight time). 

The key benefit in using a computer to present problems of this type is 
motivation. The calculations required of the user are no different than those 
in the back of a chapter in a book or those on homework assignments. How- 
ever, answering them when they are presented by the computer seems to 
make it more of a challenge and, frankly, more fun. 

The computer program checks each of your responses to see if it is within 
15% of the correct answer. If it is, your answer is considered correct. You 
may wish to change this percentage to require more or less accurate calcula- 
tions. You may also wish to change the computer calculations to round off to 
one or two decimal places. 




10 RFINBDMIZE 

2 Q = 

3 V = 6 + I N T <: 3 5 * R H D > 

40 PRINT "f\ BALL IS THROWN U 
P IJ fl R D S " 

50 PRINT "RT " J V; "METERS PER 
SECOND." , ^ 

70 PRiNf"''HDW HIGH WILL IT S 
□ < I N " 

80 PRINT M1ETERS.V5 
9 6 D S U B 5 

110 PRINT "HOW LONG UNTIL IT 

RETURNS" 
120 PRINT ^'<IN SECQNDS>"? 

13 G □ S U B 5 

14 T = 2+ I N T < 2 * V* R N D > 1 
145 IF T>R THEN 14 

15 A = V - 1 * T 

160 PRINT "UiHFiT WILL ITS VEL 
□CITY EE" 

170 PRINT " R FT ER" ? T SECONDS 

130 60 SUB 50 
190 PRINT 

20 PRINT Q; "RIGHT OUT OF 3. 

210 IF Q<:2 THEN 230 
220 PRINT " NOT BAD ! " 
230 STOP 
50 INPUT S 

5 10 IF Fl B S «C < G - R > V fl >< . 1 5 THEN 
sSo^PRINT "NOT EVEN CLOSE... 

IP 

5 3 GOTO 560 

540 PRINT "CLOSE ENOUGH-*' 

550 Q=Q+1 

560 PRINT "CORRECT RNSIJER IS 
" ? R 

57 PRINT 

580 RETURN 

9 9 9 END 



Ft BALL IS THROWN UPUIflRDS 

RT 12 METERS PER SECDNB. 

HDW HIGH WILL IT GiD CIN 

METERS)? 8 

CLOSE ENOUGH. 

CORRECT ANSWER IS 7.2 

HOW LONG UNTIL IT RETURNS 

c:iN SECONDS V? 2.5 

CLOSE ENOUGH. 

CORRECT RNSWER IS 2,4 

WHAT WILL ITS VELOCITY BE 
AFTER 2.4 SECONDS? -12 
CLOSE ENOUGH. 
CORRECT ANSWER IS-12 



RIGHT OUT DF 3. 
NOT BAD! 
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Problem Solving 



In many courses in school, the textbooks present a great variety of devices 
for solving the problems that have been neatly grouped together at the end of 
each chapter. Typically these devices consist of formulae, equations, rules and 
theorems. After a careful study of these devices, teachers give exams which 
test your ability to recall them. 

But what do you do if you are faced with the more realistic situation of not 
being told what device is likely to solve which problem or, worse yet, of hav- 
ing forgotten how to use a technique altogether? Is all hope lost? Of course 
not, although some people seem to believe that it is. 

In this chapter, several of these nasty devices mentioned above are pre- 
sented. For example, there are devices written into computer programs that 
will solve a quadratic or exponential equation and others that will calculate 
the roots and draw a plot of any function. 

However, you must remember that while these devices are useful in solving 
certain problems, the really important thing is to understand the underlying 
logic and approach. Then when it comes time to solve real world problems 
you will be better prepared to face them. Incidentally, many of the methods 
and devices presented in this chapter are used in later chapters to solve other 
kinds of problems. 
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How Many Tickets? 



Here is a problem. At a school raffle to raise money, the organizers have as 
a prize an electronic game for which they paid $18.00. To add interest to the 
raffle, the organizers have decided to sell tickets for an amount (in cents) 
equal to the number on the ticket for tickets numbered 1 to 50. For ticket 
numbers over 50, the price is 25 cents each. The organizers want to know 
how many tickets they must sell to exactly break even. 

It is probably easiest to visualize a problem of this sort with a flowchart. In 
the flowchart, T will equal the total money collected and will increase as 
more tickets are sold. The ticket number is N. When T equals or exceeds 
$18.00, N will be the answer. 

Note that the flowchart has two logical branching points (if statements in 
the program). The first compares the current ticket number to 50; if it is less, 
the ticket number is added to the total whereas if it is greater than 50, the 
total is increased by 25 cents. 




The second branch point compares the total amount collected, T, to $18.00 
(actually 1800 cents). If T is equal to or greater than 1800, the break even 
point has been reached and the values of N and T are printed (on the screen). 

A problem of this kind can be done by hand, however, because of the repet- 
itive additions it is quite tedious. Also, doing it by hand frequently leads to an 
answer of 72 rather than the correct answer of 71. Try it yourself and see 
what you get. 
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T = 

FOR N = i 

IF N<=50 

T=T+25 
c n T n "7 Q 



TO 2000 
THEN 60 




SOLD TO" 
? T / 1 



M 

1 TICKETS SOLD 
LLECT $ 18 



TD 



DONE 



Many problems can be solved quickly and correctly with a computer using 
logical analysis and a flowchart. More complex problems may have to be bro- 
ken down into additional steps and require a longer flowchart, but the ap- 
proach is fundamentally the same. 

Here are two problems for you to solve. 

The diameter of a long-playing record is 12 inches. The unused center has a 
diameter of 4 inches and there is a smooth outer edge 1/1 inch wide around 
the recording. If there are 91 grooves to the inch, how far does the needle 
move during the actual playing of the recording? 

A movie theater charges $2.50 for an adult admission and $1.00 for a child. 
At closing, the cashier counted 385 ticket stubs and had $626.50 in cash. 
How many children entered? 
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Drinking and High Blood Pressure 



This program illustrates how several simple equations can be put in a com- 
puter program to solve a more difficult overall problem. 

Here is the problem. In a survey of 1000 adults, it was found that 35 had 
high blood pressure. Of those with high blood pressure, 80% drink 15 oz. or 
more of alcohol per week. Of those without high blood pressure, 60% drink a 
similar amount. What percent of drinkers and non-drinkers have high blood 
pressure? 

This problem requires the solution of several simple equations. They could 
all be combined into one large equation, but it may be easier to understand 
the approach (and change variables later on) by writing a program with five 
separate equations. 

If H equals the number of people with high blood pressure, then H = 35 
(Line 10). Letting HI equal the number of people with high blood pressure 
who drink leads to HI = .8 x H (Line 20). 

Letting LI equal the number of people with low blood pressure who drink 
yields LI = .6 x (1000 - H). Then, the total number of drinkers, D = HI + 
LI. 

Finally, the percentage of drinkers with high blood pressure is X = HI x 
100 /D. 

The program solves the problem in a jiffy. The solution for this type of 
problem can be easily written directly in Basic without any need for a flow- 
chart or detailed analysis. Recognizing this type of problem readily will save 
a great deal of pencil pushing time. 

10 H = 3 5 

2 H 1 =: . S 3^ H 

:-iO Li=:. 6^c: 1000-H> 
40 D=H1+L1 

50 X=INT<; 10*<H*iOO.''D> >/10 
60 B=H*100/1 000 

7pF*RINT TfiBCe:- 5 ''PEOPLE WITH 
H I b H " 

3 ^ P R r N T T R B < 9 > ? " B L □ D P R E S 
90 PRINT ^'DRINKERS " ; X ? 

y. " 

100 PRINT ^'BLL PEOPLE "5E 
9 9 9 E hi D 

> 

> R LI M 



PEOPLE WITH HIGH 

BLDOD PRESSURE 
S - 

PLE 

DDhiE 



DRINKERS 5 - 7 

ALL PEOPLE 3.5 
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Here is a problem that doesn't require a single equation but makes use of 
the approaches discussed so far in this chapter. Can you solve it with three 
Basic statements? 

In early January, a shopkeeper marked down some calendars from $2.00 to 
a lower price. He sold his entire stock in one day for $603.77. How many did 
he have? 

Here's another easy one. A town in India has a population of 20,000 peo- 
ple. Five percent of them are one-legged and half of the others go barefoot. 
How many sandals are worn in the town? 




Two Simultaneous Equations 



So far in this chapter, only problems with linear equations have been 
considered. But the computer can be used to solve much more difficult equa- 
tions. In fact, it is in problems involving second and third degree equations, 
exponentials, and the like where the computer really starts to pay off. Con- 
sider the following two simultaneous equations: 

2X ^ 16y 3^ = 27y 

It is not at all easy to solve these two equations by hand. But a simple Basic 
program can be written to solve the equations using trial and error. This is 
sometimes referred to as a brute force approach because every possible 
combination of numbers between an upper and lower limit is tried until a 
solution is reached or until the program runs out of values. 



10 
20 
30 
40 
45 

» m 

50 
60 
70 
80 
90 
W " 
1 

IB 

99 



FDR X=l TO 100 

FOR = i T D 100 

IF 2->X<>16*V/'3 THEN 70 

IF 3-^K<>27*V THEN 70 

PRINT "FINRLLY SOLVED IT. 

II 

D D T KJ T = " X ? " V = " ; Y 



PRINT 
STOP 
NEXT Y 
N E X T X 

PRINT^^'NO INTEGER 
PRINT ^'BETWEEN 1 
9 END 



>RUN 

FINRLLY SOLVED IT. . . 
X= 4 Y= 3 



SOLUTIO 
AND 100 



In this particular case, a solution is reached rather quickly with x = 4 and 
y = 3. However, if in the second equation the y coefficient is changed shghtly 
from 27 to 28, the computer will try 10,000 possible solutions before finally 
concluding that no integer solution exists — at least within the range of to 
100. Warning: this will run for a very long time. 
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>40 IF 3-^X<>28^&Y THEN 70 
>RUN 

NO INTEGER SOLUTION 
BETWEEN 1 AND 100. 

** DONE ** 



Although the trial and error (brute force) approach is widely used, it is 
highly inefficient. In general, a systematic or guided trial and error approach 
is preferable to one that simply tries every possible solution. However, for 
some problems the simple "try every value" approach may be appropriate. (A 
comprehensive discussion of trial and error approaches can be found on pp 
36-40 of Computers in Mathematics: A Sourcebook of Ideas.) 

Three problem solving approaches have been discussed so far. Remember- 
ing them, how would you do this problem? A boy and his sister visited a farm 
where they saw a pen filled with pigs and chickens. When they returned 
home, the boy observed that there were 18 animals in all, and his sister re- 
ported that she had counted a total of 50 legs. How many pigs were there in 
the pen? 
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Quadratic Equation Solver 



For any values A, B, and C of a first degree quadratic equation (Ax^ + Bx 
H- C = 0), this program will compute the roots of the equation. The solution 
is based on the quadratic theorem which solves for roots with the following 
formula: 



Assuming that A, B, and C are real numbers, the following principles 
apply: 

1. If B^ - 4AC is positive, then the roots are real and unequal. 

2. If B^ - 4AC equals 0, then the roots are real and equal. 

3. If B^ - 4AC is negative, then the roots are imaginary and unequal. 
The program takes into account all these possibilities and correctly identi- 
fies the type of roots along with their values for any set of coefficients. 

Is this program useful by itself? Except for solving quadratic equations for 
algebra class, probably not. However, as a routine in a larger program to 
solve quadratic equations that might be encountered, it could be very useful. 



X = 



-b+Vb^^ac 

2A 



iO PRINT " QUftDRRT .1 C EQURTiaN 

'5 n 1 "J c D »• 



SOLVER 




SP E^'IH"^ '^PLEASE ENTER THE F 
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> R U M 

QUflDRRTIC EQUFITIDN SOLVER 
PLERSE ENTER THE FOLLDWING: 
VRLUE DF R=l 
VALUE OF B=2 
VALUE DF C=l 

THE ROOTS ARE 
-1 

*« DONE ** 

> R U N 

QURDRflTIC EQUATION SOLVER. 
PLERSE ENTER THE FOLLOWING: 
VALUE OF A=0 
VALUE OF B=2 
VALUE OF C=4 

THAT IS fl FIRST 
DEGREE EQUATION 

** IiDNE #* 

> R U N 

QUADRRTIC EQUATION SOLVER 
PLERSE ENTER THE FOLLGWINGs 
VALUE DF fl=2 
VALUE OF B=5 
VALUE DF C=2 

THE ROOTS ARE 

-i 

** DONE ** 

> R U N 

QUADRATIC EQUATION SOLVER 
PLEASE ENTER THE FOLLOWING: 
VALUE OF ft=4 
VALUE OF E=2 
VALUE DF C=4 

THE RDDTS ARE IMAGINARY 
-4 + 15.491933 39 *I 
-4 - 15.4919333S *I 

** DONE ** 
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Exponential Equation Solver 



Another general routine for solving a particular type of equation is this one 
to solve for an exponent in an exponential equation. 

Given the values A, B, m, and n, this program will solve for x in any 
exponential equation of the form: 

For example, the program will solve any of the following problems: 

1. 5^ = 40 

2. 53X+> = 7.6 

3. 17^-3 = 8.12 

4. 1P-2X = 247 

If you were to solve an exponential equation by hand, you would probably 
go through the following steps: 

5^ = 40 

log 5^ = log 40 

X log 5 = log 40 

X = log 40/log 5 

X = 1.6021/. 6990 = 2.292 
However, in more generalized form, the solution for x is: 

log B 

X = A 
M-(N/M) 

Note that the program to solve this problem is actually divided into two 
sub-programs. The first is a data loader program. It requires that data be en- 
tered in the following order: A, B, m, and n for each equation to be solved. If 
a coefficient is not present, it must be entered as a zero. 

The data for the four problems listed above were entered into the two 
data statements (50 and 60). 

The program as it is presented here can be improved in several ways. First, 
it always solves the same four equations. How can you generalize it to solve 
for other equations? Second, if you were to make use of this routine in an- 
other program, you would probably not be able to use a read statement; how 
could you get rid of it?. 
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10 PRINT "EXPDNEfilTIRL EQURTI 
N S D L V E R 

I5^PRINT"fl B M N 

2 ' R E ft D R , B , f1 , H 

3 X = < L □ G C B ) L □ G C R > > / M - < N . n > 
3 5 X = I W T < . 5 + 1 * X 1 

40 PRINT R? TRBCS) 5 B? TfiB< 10> 5 
M 5 TAB <: 1 5 > ; N 5 TRB c: 20 > ? X 
45 GOTO 20 

50 DftTR 5 . 40 , 1 r . 5 , 7 . 6 r 3 , 1 
60 DfiTft 17.8.12,1,-3.11,247, 
— 2 » 1 
99 END 



> R U N 

EXPDNENTIRL EQURTIDN lEiOLVER 

fl B M N X 

5 4 1 2.29 

5 7,63 1 ,09 

17 8.12 

1 1 247 -2 "l -^'hV 



* DRTR ERROR IN 20 
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Roots of Any Function 



This program will find the roots of a function, any function! The function 
may be linear, quadratic, cubic, trigonometric or any combination as long at 
it can be represented in the Basic language. The program as it appears here 
finds the roots between -20 and 20 although you can change these boundaries 
in Statement 120. 

The method used involves evaluating the function at small incremental 
intervals, finding places where the value of the function changes sign and 
then, by successive approximations, finding the zero point. This approach 
borrows from Newton's method in the final narrowing down but, unlike 
Newton's method, will not fail to converge in the event one makes an unlucky 
first guess. 

Before running the program you must first type in your function in State- 
ment 100. For example, 

DEF FNA (X) = 2*XT3+1 1 *XT2-31 *X- 180 

DEF FNA (X) = X-4 

DEF FNA (X) =SIN(X)-.5 

You may have to refer to the Basic manual with your system to see exactly 
how a function should be stated. 

The routine used in this program is very powerful and could possibly be 
used as a subroutine in many other programs. 

How can you use this program to help you solve this problem for x? 




26 



100 DEF FNfl X > =2*K ''3 + 1 1 2- 

.3 1 # X - 1 B 

110 21=-20 

120 FOR I=-19.9 TO 20 

130 IF SGNCFNflC I > ;'=SeN<:FNfl< I 

+ 1 ;0 THEN 300 

140 K=I 

150 J =1+1 

160 IF FNfl'::K)<FNl=K.J>THEN 200 

170 2=K. 

18 K = J 

19 J = Z 

20 Z=<K + .J>'^2 

210 IF FNRCZXO THEN 240 

2 2 J = Z 

230 GOTO 250 

240 K=Z 

250 IF flBS<:FNfl<;Z> ;■>. 000005 T 

HEN 200 

2 6 Z = SG N C Z ^ * I N T < RE S <; Z > * I 000 

+ . 5 ) r 1 

2?0 IF Z=21 THEN 300 

28 PRINT "F ( " 5 Z ; " ) = " 

290 21=Z 

300 NEXT I 

999 END 

R IJ N 

F - 5 > = 
F - 4 . 5 :> = 
F<: 4 >=0 

** DDNE ** 

100 BEF FNf=l<X. > = X-4 
RUN 

F'C 4 >=0 
*« DONE ** 

100 DEF FNfi':;X)=SIN(:X)-. 5 
RUN 

F < - 1 S . 3 2 6 > = 
F<- 16. 2316 >=0 
F-C-12.04 28 >=0 
F':-9.9484 >=0 
F<:-5.7596 >=0 
F<:-3.6652 :)=o 
F<: .5236 >=0 
F < 2.618 > = 
F 't 6.8068 > =0 
F < 8.9012 :) = 
F< 13.09 -j = 
FC 15. 1844 >=0 
FC 19.3732 >=0 

** DONE ** 
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Plot Any Function 



Here is a nifty program that will produce a plot of any function on the 
screen or printer. 

Before running this program, you must type in your function in line 200. 
Like the previous program, which finds the roots of any function, this one 
will plot any function. You must tell the program between what values you 
want the function plotted, i.e., a minimum and maximum value of the x co- 
ordinate. You also input the x plotting increment you wish. 

As it appears here, the program does not allow the user to select the y co- 
ordinates; the program plots y values between -30 and 30. 

Here are several functions you might want to try plotting: 



Function 




X Limits 


Increment 


DEF FNA(X) 


= 2*X 


- 15 15 


1 


DEF FNA(X) 


= 30*SIN(X) 


-5 5 


. 25 


DEF FNA(X) 


= X-X*X 


-5 6 


. 5 


DEF FNA(X) 


= 30*EXP( -X*X/ 100) 


-30 30 


1 . 5 


DEF FNA(X) 


= X*X-X 


-5 6 


1 


DEF FNA(X) 


= XT2-X-15 







The last function listed is the one plotted in the sample run with the pro- 
gram. Exponential functions are a great deal of fun and sometimes lead to 
unexpected and interesting results, particularly when combined with trigo- 
nometric functions. Experiment! Have fun! 
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2 D E F F N n < X ) = X ^- 2 - X - 1 5 
2 10 PRINT 

"INITIRL VnLUE 



220 INPUT 
• : X 1 

INPUT 

2 

I H P U T 

S 

PRINT 
PRINT 



230 
I) ■ V 
• ft 

24 

250 
290 
300 
3 1 
6 
32 
322 

32 5 

33 
340 

a n I 
r I 

350 
360 



□ F 

"FINAL VALUE OF K= 
"INCREMENTAL STEP 



. TAB < 1 2 > 5 " X = " ? X 1 

FOR X=X1 TO X2 STEP -S 
IF fi B S t: X )< . 1 THEN 



Y = FrvlAt:X>x2+14 

IF Y<:28 THEN 330 

IF'^Y>i4 THEN 410 
PRINT TRBC V> ? " ^ 



TAB t; 1 4 > 



370 
+ 
38 
390 
4 10 
? » * 
420 
430 
999 



GOTO 420 

PRINT 

+ Y " 

PRINT "-30 -2 -10 +1 
20 +30" 

"r^ S5 'x* — 3 

60TD 420 

PRINT T A B < 1 4 > S " ! " ; T A B < V > 



NEXT X 
PRINT TftB<i 
END 



> 5 " X = " 5 X 2 



INITIAL VALUE^DF X=7 
FINAL VIRLUE OF X = -6 
INCREMENTAL STEP =1 

X= 7 



-Y +Y 

-30 -20 -10 +10 +20 +30 
* 



X = -6 



** DONE ** 
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3 

Sets and Repetitive Trials 



For solving relatively simple problems, the computer may not be any help 
at all. In fact, it may take more time to write a program to solve a problem 
than it would to solve it by hand or with a calculator. This chapter should 
help you recognize problems that are suitable for computer solution and 
those that are not. 

In some problems, you may think that the computer will not be any help. 
However, one thing writing a computer program will always do is force you 
to reason out the approach to solving the problem logically and precisely. The 
computer can't solve problems unless it is told exactly how to proceed; hence 
you must understand a problem completely before you can program it for the 
computer. 

Several of the sections in this chapter discuss sets of data. While the sets 
used in the examples have relatively few elements or values, you should bear 
in mind that real world problems often have thousands or millions of pieces 
of data, and the only practical way to solve problems of this size is with a 
computer. For example, consider how you would most efficiently schedule 
the shipments on a railroad train leaving Boston with 4000 diverse cargoes 
bound for Phoenix and 780 points in between. Now consider that there are 
200 freight trains per day leaving Boston. Now add to that the 10,000 other 
trains leaving other cities every day that must use the same network of track 
and you can see that dealing with real world sets of data is no easy task. 
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Group of Girls and Boys 



In the previous chapter we said that there may be some problems for which 
brute force trial and error is appropriate. This would be the case if the prob- 
lems were relatively small and trying every possible solution would not tie up 
a great deal of valuable computer time. Here is a problem involving two lin- 
ear equations that lends itself to a trial and error approach. 

The problem is as follows: When 15 girls leave a group of boys and girls, 
there are two boys for every girl (lucky girls). Next, 45 boys decide to leave; 
then there are 5 girls for every boy (lucky boys!). How many girls were there 
in the group before anyone left? 

Before rushing to the computer, you must recognize that this problem re- 
quires the solution of two simultaneous equations. If G equals the original 
number of girls and B the original number of boys, then the two equations 
are: 

(G- 15) X 2 = B 



(B-45)x5 = (G- 15) 



10 1=0 

20 FDR G=l TO 100 

3 FDR B=l TO 10 

40 1=1+1 

50 IF 2*c;6-l.5><>B THEN 90 

6 IF 5 * < B - 4 5 > < >< G - 1 5 > T H E N 3 


7 PRINT G ? " G I R L S ^ " ? B 5 " B □ Y '5 
li 

8 G D T 12 

9 N E >i T B 

100 ne^:t 6 

110 print "no solution betwe 

EN 1 8< 10 0" 

12 PRINT I ? " C □ M B 1 N fl T I □ N S T R 

I ED " 

9 9 9 END 



^RUN 

4 GIRLS. 50 BOYS 

39 5 CDMEir^ATIGNS TRIED 



The computer program uses two FOR loops (Statements 20-100 and 30-90) 
to try every combination of values for B and G between 1 and 100 until a 
solution is found or until the program runs out of values. The variable I 
(Statement 40) is a counter which records the number of trials required to 
reach a solution. 
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The program is straightforward and finds a solution after 3950 trials. How- 
ever, it would have been a simple matter to substitute the value of B from the 
first equation in the second one and quickly solve the problem by hand or 
with the aid of a calculator. It is important to recognize that if a problem can 
easily be solved by other methods, the computer offers little or no advantage. 

Try this problem. You may or may not want to use your computer. If Mat- 
thew can beat Jeff by one-tenth of a mile in a two-mile race and Jeff can beat 
Steven by one-fifth of a mile in a two-mile race, by what distance could Mat- 
thew beat Steven in a two-mile race? (Hint: the answer is not 3/10 mile.) 
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Brown's Books 



The use of a trial and error approach can generally be improved signifi- 
cantly if the combinations to be tried can be narrowed down in some way. 
The solution to this problem illustrates how the speed of obtaining a solution 
can be improved well over 100 fold by combining equations and eliminating 
certain solution possibilities. 

Here is the problem. Brown sold 48 books at a flea market, some for $3 
each, some for $5 each and others for $8. He collected a total of $175. He 
remembered having an even number of $5 books. Can you determine how 
many of each kind of book he had? 

The equations for solution are (letting T equal the number of $3 books, F 
the number of $5 books, and E the number of $8 books): 

T + F + E = 48 

3*T + 5*F + 8*E = 175 

The first program was written simply to try all possible combinations of T, 
F, and E from 1 to 48. It yields three solutions for the problem, although the 
two solutions with an odd number of $5 books can be eliminated leaving just 
the one desired solution. 

10 PRINT "BOOKS PROBLEM" 
30 PRINT "'^^^ SOLUTIONS" 

50 PRINT 

60 FDR T=l TO 48 

70 FDR F=l TO 43 

SO FDR E=l TO 48 

?Sn ^?cr'T + F + E 2< >48 THEN 120 

EN 1^0 ^^■^•^ + ^*5 + E^f8)<>175 TH 

110 PRINT T? TflBc: 10 ) ? F; TRB.:20 

12 NEXT E 

130 NEXT F 

140 NEXT T 

,999 END 

>RUN 
BOOKS PROBLEM 
fiLL SOLUTIONS 

« 1 :$ 5 $ 8 

IP i i 

34 13 1 

** DONE ** 
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This program took approximately 22 minutes and 13 seconds to run on the 
TI 99/2 computer. A typical minicomputer (PDP-8/e) could run this prob- 
lem in about 7.3 seconds. In either case, this is a long time to tie up the 
computer. 

It is rather easy to combine the two equations into one by solving for T. 
The single equation is then: 

2*F +5*E = 31 

In this equation, the limits can be reduced (from 48 used in the first run) 
since F cannot possibly be greater than 31/2 or 15.5 and E cannot be greater 
than 31/5 or 6.2. Making the appropriate program modifications leads to the 
second program. 

>70 FDR F=l TO 
>S0 FDR E=l TD^ 
>90 IF CF*2+E*5 

MOO 

>105 T=*48-F-E 
>140 

> 

>RUN 
BDQKS PROBLECJ 
RLL SOLUTIONS 

ii:10 S5 

40 3 
37 8 
34 13 

** DONE ** 

Using this program produces a dramatic improvement in the time to solu- 
tion. On the TI 99/2, the time is approximately 2.5 seconds and on the PDP- 
8, about 0.16 seconds. 

Since the problem states that F must be even, a final modification which 
steps F by two in Statement 20, can be made. This version of the program 
takes only 1.25 seconds to run on the Timex and 0.06 seconds to run on the 
PDP-8. 

>20 PRU4T "SOLUTIONS, EVEN Vfi 

L U E D F F " 
>70 FDR F=2 TO 14 STEP 2 

Notice the enormous improvement in computing time required for a solu- 
tion, over 1000 fold on the TI 99/2 and 100 fold on a PDP-8. Brute force 
certainly is inefficient! It is generally worthwhile to think through most prob- 
lems, particularly big ones, before rushing to the computer. The computer 
may be fast, but we just improved its performance by 1000 times by using a 
Uttle common sense. 



> <> 3 1 



THEN 120 



S8 

5 
3 
1 
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Intersection of Sets 



Two sets of numbers can be combined to yield a third set by the operation 
of intersection. The intersection of two sets A and B is the set that contains all 
elements that belong to both A and B. It does not contain any other elements. 
The intersection is usually written A n B. 

For example if M = < 0,2,4,6 > and K = < 1,2,3,4 > , then M n K = 
< 2,4 >. 

This program finds the intersection of two sets of numbers. It has been 
written to find the intersection of the two repetitive sets described in State- 
ments 30 and 40. In the sample run, Statement 30 describes the set 
X = < 1,3,5, ... 19 > and Statement 40 describes the set 
y = < 2,5,8, . . . 29 > . 
Notice that successive values of x increase by 2 and y by 3. 

lOgPRINT "THE INTERSECTION □ 

20 PRINT "X RND Y ISs " 

30 FOR X=l TD 19 STEP 2 

4 FDR Y=2 TD 2 9 STEP 3 

50 IF X=Y THEN 100 

60 NEXT Y 

70 NEXT X 

80 STOP 

10 PRINT X 

110 GOTO 70 

99 9 END 

> 

1 1 
1 7 

* * DONE * * 



However, if the set cannot be so neatly described, it may be desirable to 
rewrite the program to examine any set of data. This is done with the read 
statement which reads into x the data points in the data statement.The pro- 
gram is set up to use the same y set as the first program, but the x set is de- 
fined in the data statement. 

You should be able to see from the first combination of sets that if there is a 
numerical pattern in the sets which intersect, then there is also a pattern in 
the resulting intersecting set. In the example, the x values increase by 2 and 
the y values by 3, hence the values in the intersecting set increase by 2 x 3 = 
6. Although the intersection of these sets could easily be calculated by hand, 
the computer can be an aid in evaluating more complicated sets. 



36 



10 PRINT "THE IHTERSECTIDN 
F SETS" 

2 PRINT "X RND Y IS: " 
30 RERD X 

40 FOR Y=2 TO 29 STEP 3 

50 IF X=Y THEN 100 

6 NEXT Y 

70 GOTO 30 

100 PRINT X 

110 GOTO 30 

2 00 DRTR 2,3,8,9,14»15,20,2l 

,26,27 

999 END 

> 
> 

'*> R IJ N 

THE INTERSECTION OF SETS 
X AND Y IS: 

2 

8 

1 4 

20 

2 6 

* DflTR ERROR IN 30 
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Prime Factors 



A prime factor is a positive integer that has no factor except itself and one. 
The first ten prime factors (or numbers) are 2, 3, 5, 7, 11, 13, 17, 19, 23, and 
29. The definition gives the basic method for determining whether a number 
is prime: divide by all smaller integers down to 2, testing whether the remain- 
der is zero for at least one of them. If not, the number is prime. 

But this is highly inefficient. It is obvious that a number is not prime if it is 
any even number greater than 2; hence only odd divisors need to be tried. 
Also, it is not necessary to try divisors greater than the square root of the 
number. 

Since the division method is inefficient, various schemes have been devised 
to avoid division. The basic idea underlying all such schemes is called the 
sieve of Eratosthenes (276 B.C.-195 B.C.). Imagine a list of odd numbers 
from 3 up. Strike out every third number after 3, every fifth number after 5, 
and so on. This will leave only prime numbers. 



1 

20 

3 

R I M 

40 

E6E 

50 

6 

STO 

90 

100 

1 1 

120 

130 

135 

140 

150 

160 

1 ?0 
180 
190 

2 

2 10 
22 
230 
24 
250 
260 
270 
280 
290 
300 

3 10 
320 
33 



DIM R 

c = o 

PRINT 
E " 

PRINT 
R " 

PRINT 
PR INT 
P " 

PR I NT 
INPLI 
X = 
M = I N 
IF M 
I F M 
PR I N 
GOTO 
P R I N 
I = I 
1 = 1 + 
IF I 
IF M 
X = X + 
fi C X ) 
M = M/ 
GOTO 
I F X 
FOR 
PR I N 
NEXT 
GOTO 
PR I N 
PR I N 
GOTO 
STOP 



C 1 40 > 

"PROGRflH COMPUTES P 
"FFICTDRS OF ANY INT 

"ENTER <2ER0> TO 

T "YOUR NUMBER ="sf1 
T < M ) 

>0 THEN 160 
=0 THEN 330 
T "MUST BE POSITIVE" 

1 
T M5 



>M THEN 250 
^I<:>INT<I1/'I :)THEN 180 

= I 
I 

2 00 

= I THEN 300 
L=l TD X 
T n < L > ; 
L 

3 1 

T "IS PRIME" 
9 
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PRDGRRM CDhPUTE'B PRl riE 
FACTORS OF RNV INTEbER 

ENTER (ZERO> TO STOP 



YOUR 
105 



NUMBER 

3 5 



= 105 
7 



YOUR NUMBER^=72 
72 2 2 2 3 

YOUR NUMBER =89 
89 IS PRIME 



YOUR NUMBER^=47 
47 IS PRIME 



YOUR NUMBER =0 
** DONE 



The program here finds the prime factors of any integer, or prints out "N is 
prime" if the integer has no proper divisors. 

Run this program for a large number of different integers and see if you 
can discover relationships between numbers and their prime factors. You 
should also try to figure out the method employed in the program to find the 
prime factors of any integer. To do this, you might want to draw a flowchart 
to show what is happening in the program. This will help you see the method 
used to find a prime factor and might help you in writing a program to gen- 
erate primes. 

In writing a program to generate prime factors, you can use the sieve 
method. However, as the numbers become very large, you will have to figure 
out a way to represent integers with more digits than your computer can han- 
dle at one time. (One approach is described on pp. 19-21 of Computers in 
Mathematics.) 

Goldbach was a mathematician who made a conjecture that every even 
number greater than 4 can be written as the sum of two prime numbers (16 = 
11 + 5, 30 = 17 + 13, etc.). No one has ever proved it but no one has dis- 
proved ii either. That is why it is called a conjecture. Can you write a pro- 
gram that will prove or disprove this conjecture? Or how about writing a 
program to prove Goldbach's conjecture for even numbers up to 50? You 
should be able to write this program with 12 or fewer statements. 

Here is another problem involving prime numbers. Assume a life span of 
80 years. In what year of the 20th century (1900-1999) would a person have 
to be born to have the maximum number of birthdays occurring in prime 
years? The minimum number? 
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Greatest Common Divisor 



The greatest common divisor of a set of numbers is, as its name implies, the 
greatest integer that will divide into a set of two or more numbers. For exam- 
ple, the set of numbers 12, 20, and 28 have a greatest common divisor of 4. 
Nothing larger than 4 will divide evenly into all three numbers. 

This program will find the greatest common divisor for any set of integers. 
To run it, you simply input the number of integers in your set, type them in 
when requested and let the program calculate the GCD. The heart of the 
calculation is in Statement 180. 

Do you know the meaning of a relatively prime set of numbers? Can you 
figure out the meaning from the third sample run of the program or from 
runs of your own? How is a set of relatively prime numbers different from a 
set of prime factors? Can you find a set of 10 integers that is relatively prime? 



I 

RER 

20 

30 

4 

N 

50 
60 
70 
80 
9 
100 
1 1 
120 
130 
140 
150 
160 
170 
180 
EN 
1 9 
200 
210 
220 
230 
24 
I ME 
250 
260 
270 
230 
N > 
290 
300 
3 10 



PR I NT 
TEST" 
PRINT 
PRINT 
INPUT 

N = I N T c: 

Din y. c 

PRINT 

FDR K = 
INPUT 
IF Xc: 
S = X < K 
NEXT 
PRINT 
G = 
FDR M 
FDR I 
IF X< 

2 1 
NEXT 
G = M 
NEXT 
PRINT 
IF 6> 
P R I N T 

GDTD 
PRINT 
PRINT 
PRINT 

Jl m 

INPUT 
IF R« 
STOP 



"PROGRftn CDMPUTES G 
" CatlMDN D I V I SDR . " 
"NUMBERS IN SET =": 
N > 

100 ) 

"ENTER NUMBERS" 

I TON 

X •:: K > 
K ) > S THEN 130 
> 
K. 



= 2 TD S 
= 1 TD N 

I / M < > I N T C X ( I > f1 T H 

I 

N 

NUMBERS "; 
THEN 26 
"RRE RELATIVELY PR 

27 

" G . C . D . = " 5 G 

"AND THE R SET <Y DR 

R S 

="Y" THEN 3 
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>RUN 
PRDGRRM 
COM M D N D 


COMPUTES GREATEST 
I V I S □ R . 


H 1 J ri B £ R S 
EMTER NU 
? 1 2 
? 36 
? 9 6 


IN SET 

til T> P P C 
1 1 JJ C C% 


= 3 


N UMBERS 


G . C . D . 


= 1 2 


fi N D T H E R 


SET CY 


□ R M > ? Y 


NUMBERS 
ENTER N U 
? 2 
? 36 
? 97 


IN SET 
MBERS 


= 3 


NUMBERS 


ARE RELATIVELY PRIME 


RNDTHER 


SET C Y 


□R N>? ■ 



In the last section we discussed prime numbers. Here is an interesting chal- 
lenge for you involving prime numbers. Until late 1982, the longest pro- 
gression of prime numbers in which all differed by the same number was 17. 
Prof. Paul Pritchard in the computer science department at Cornell Univer- 
sity wrote a program to determine if there was a longer progression. Using a 
DEC VAX-1 1/780, he found the string of 18 numbers shown below. He also 
discovered fourteen other 17-number progressions and ten 18-number pro- 
gressions, but none yet with 19 numbers. He beUeves there is at least one; can 
you find it? 



107928278317 197233324147 

117851061187 207156107017 

127773844057 217078889887 

1 37696626927 227001 672757 

1 476 1 9409797 236924455627 

1 57542 1 92667 246847238497 

1 67464975537 25677002 1 367 

1 77387758407 266692804237 

187310541277 276615587107 
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Cryptarithmetic Problems 



Cryptarithmetic or alphametic problems are arithmetic expressions in 
which the digits are replaced by letters of the alphabet. Each digit is asso- 
ciated with a letter to produce an interesting statement, for example: 

SEND 
+ MORE 

MONEY 

If the college student who sent this message to his father needed $106.52 
for plane fare home, this was the right message to send since this combination 
of letters has one unique solution, in particular: 

9567 
+ 1085 

10652 

However, if the student let things go to the last moment and was in more of 
a rush, he might have reworded the message: 

WIRE 
+ MOR E 
MONEY 

In this case, how much should his dad send? Earlier in the book, trial and 
error approaches to solving problems were discussed. It was noted that the 
brute force approach of trying every alternative was sometimes appropriate. 
Is it in this case? 

No! The number of possible alternative solutions is the factorial of the 
number of different letters in the alphametic expression, i.e., 8! or 40,320. A 
program to try out every one of these possibilities would run for a 1-o-n-g 
time. 

In this case it is much more efficient to apply some common sense to nar- 
row down the number of alternatives. The best approach to this process is to 
divide up the search space into large classes (or sets), according to a common 
property shared by members of each class, and then attempt to eliminate en- 
tire classes by the method of contradiction. 

Consider the "WIRE + MORE = MONEY" problem. Can E = 0? Since 
E + E = Y, Y must also equal 0, contradicting the fact that Y and E must be 
different digits. Thus, the entire class of solutions in which E = can be 
ruled out. 

Consider E = 3. Now Y = 6 and there is no carry to the next column. So 
in this column R + R=:EorE + 10 if a carry is involved. But in either case 
E must be an even number since 2R is always even; this contradicts the 
assumption that E = 3. 

By following this type of classificatory contradiction process for each of the 
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digits in the order E, R, I, O and N, the computer program will search out all 
possible solutions to the problem. Unlike the "send more money" problem, 
the "wire more money" problem has five possible solutions. A smart father 
would choose the lowest solution and wire his son $103.48. 





I 


R 


E " 


+ M 


□ 


R 


E " 

CI 


□ 


N 


E 


Y " 


TO 3 








THEN 


1 5 








10 PRINT " P R D 6 R R ri FINDS H L L 
SOLUTIONS" 

20 PRINT "TO THE FOLLOWING I 
R B L E M : " 
3 PRINT 
40 PRINT " 

5 PRINT " 

6 PRINT " 

7 PRINT " M 

8 PRINT 

9 M=l 
100 FDR E=2 
110 Y = E + E 
120 IF Y>tO 

13 C 1 = 

14 GOTO 170 
150 Cl=l 

184 IF R = Y THEN 470^^ 

190 IF R+R+C1=E THEN 220 

200 IF R+R+C1=E+10 THEN 240 

210 GOTO 470 

22 C2=0 

230 GOTO 250 

2 40 C2=l . 

250 FDR 1=0 TO 9 

260 IF I=M THEN 460 

262 IF I=E THEN 460 

264 IF I=Y THEN 460 

266 IF I=R THEN 460 

27 FOR 0=0 TO 9 

2 80 IF D=N THEN 4 50 

282 IF D=E THEN 450 

284 IF D=V THEN 45y 

286 IF 0=R THEN 450 

238 IF 0=1 THEN 450 

2 90 N=I+0+C2 

300 IF N>=10 THEN 330 

310 C3 = 

32 GOTO 35 

3 3 C 3 = 1 
340 N=N-10 

350 IF N=M THEN 450 
360 FOR W=0 TD 9 
362 IF W=M THEN 440 
364 IF W=E THEN 440 
36 6 IF W=Y THEN 44 

36 8 IF W=R THEN 440 
370 IF U=I THEN 44g 

37 2 IF U=D THEN 440 
374 IF U=N THEN 440 

330 IF 0+ 10< MJ + M + C3 THEN 440 
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3 9 PRINT " " 5 y ? I ? R 5 E 
400 PRINT •■ ";Mj05R;E 



4 10 PRINT 

420 PRINT M5D5N5E5Y 

4 3 PRINT 

4 4 NEXT W 

450 NEXT □ 

46 NEXT I 

470 NEXT R 

480 NEXT E 

45 STOP 



PROGRftM FINDS ALL SOLUTIDNS 
TD THE FOLLOWING PROBLEM: 



+ 


111 
M 


T 
X 

□ 


P 

Pv 

R 


UllU { 


M 


n 


N 


E 


Y 




Q 

1 


r 




O 

6 


2 
2 


1 





8 


2 


4 




3 
1 


2 



7 
7 


4 

4 


1 





3 


4 


8 




■9 
1 


5 



7 
7 


4 
4 


1 





6 


4 


8 




9 
1 


2 



8 
3 


7 
7 


1 





3 


7 


4 




■9 
1 


5 



8 
8 


7 
7 


1 





6 


7 


4 



DONE 



There are other approaches to solving cryptarithmetic problems, but all of 
them benefit greatly from reducing the search space as much as possible be- 
fore putting the problem on the computer. See if you can devise another 
successful approach and write a program to implement it. 

44 



Here are some problems for you to try. 

DONALD TWO 
+ GERALD X TWO 



ROBERT 

ABC 
X DE 

FEC 
DEC 

HGBC 



THREE 

THE 
EARTH 
VENUS 
SATURN 
+ URANUS 

NEPTUNE 



ABODE 
X 4 

EDCBA 

SPRING 
RAINS 
BRING 
+ GREEN 

PLAINS 



VIOLIN + VIOLIN + VIOLA + CELLO = QUARTET 
THREE + NINE = EIGHT + FOUR 



ONE 
TWO 
+ FIVE 

EIGHT 



FORTY 
TEN 
+ TEN 

SIXTY 



FIVE 
-FOUR 

ONE 
+ ONE 

TWO 
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Sailors and Monkey Problem 



There are many variations of the sailors and monkey problem. Here is one 
of them. 

Five sailors and a monkey were on an island. One evening the sailors 
rounded up all the coconuts they could find and put them in a large pile. Be- 
ing exhausted from working so hard, they decided to wait and divide them up 
equally in the morning. During the night, a sailor awoke and separated the 
nuts into five equal piles, but had one nut left over which he gave to the mon- 
key. He took one pile, hid it, and pushed the other four together and went 
back to sleep. He was followed in this action by the other four sailors, each of 
whom did exactly the same thing. Next morning the remaining nuts were di- 
vided equally with one remaining nut going to the monkey. What is the small- 
est number of coconuts with which they could have begun? 

Although there is an elegant algebraic solution to this problem, a more 
suitable approach for the computer is that of working backwards. A typical 
solution to a problem can be thought of as a path that leads from the given 
information to the goal. However, in this case the goal, or final state, is 
known, thus it is easier to start there and work backwards to the initial state. 

As mentioned at the outset, many sailor and monkey problems exist, in 
fact, an infinite number of them. For example, instead of five sailors, there 
could be three or six or 14. Thus it is desirable to devise a general solution in- 
stead of just one to solve one specific problem. 

In the flowchart and computer program, S is the number of sailors and A is 
the number of coconuts that each sailor received in the final division of the 
pile. Since one coconut was given to the monkey at each division, the total 
number of coconuts left in the morning must be S x A + 1. But this pile came 
from pushing together S - 1 equal piles. Thus, the key condition that must 
hold for (S X A + 1) / (S - 1) to be an integer K, which represents the num- 
ber of coconuts that the last sailor stole from a pile of S x K + 1 coconuts. 
But this pile is the result of pushing together S - 1 equal piles by the previous 
thief, so again (S x K + 1) / (S - 1) is an integer and so on back through all S 
raids on the pile. 

Note in the flowchart (and program) that the first trial value for A is 1 
(Statement 70). In Statement 1 10 this value is increased by 1 until the value of 
(S X K + 1)/ (S - 1) is an integer as tested for in Statement 100. This process 
is then continued until the counter for the second loop, N (nighttime pile di- 
visions) equals the number of sailors. 

Although the program will work for any number of sailors, it takes a fairly 
long time to run for more than five. Remembering what you have learned ear- 
lier in this chapter, can you devise a way to make the program more efficient? 



46 




LOOP 



1 


A" 


1 












^ 

• o 



S-l 



PRlt4T 



Loop 
I 
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10 PRINT ^' P R G R R M SOLVES SRI 
LOR'- fi N D " 

20 PRINT " MONKEY PROELEn BY 
W R K I N 5 " 

30 PRINT "BnCKWFlRDS. " 

50 INPUT "NUMBER OF SRILnRS= 

" ■ s 

6 0' PRINT 

7 R=l 
QO K = fl 

lOo'^IF CS*K+l)/-(:S-l> = INTCi:S* 

K+1 >><:S-1:''jTHEN 130 

1 10 fi=fl+l . . 

120 GOTO "SO 

130 N=N+1. ^ ^ _ 

14 K = < S * K + 1 > 'M. S - 1 > 
150 4F H=S THEN 170 

170 PRINT "THE FEWEST CnCDNU 
T S T H fl T " 5 S 

1130 PRINT "SfilLDRS CAN HAVE 
T □ E E S I N " 

190 PRINT "IS";S*K+1 
20 PRINT 

2 10 PRINT "IN THE MORNING, E 

FICH SAILOR" 

220 PRINT "eETS"5fi 

230 STOP 

";> 
> 

>RUN 

PRGGRfiM SOLVES SfllLDRS AND 
MONKEY PROBLEM BY WORKING 
BACKWARDS. 

NUMBER GF SftlLDRS^S 

THE FEWEST COCONUTS THAT 3 
SAILORS CAN HAVE TO BEGIN 

15 79 

IN THE MORNING, EACH SAILDR 
GETS 7 

* * DONE .* * 

> R U N 

PROGRAM SOLVES SAILORS AND 
MONKEY PROBLEM BY WORKING 
BACKWARDS. 

NUMBER OF SflILDRS=5 

THE FEWEST COCONUTS THAT 5 
SAILORS CAN HAVE TO BEGIN 
IS 15 62 1 

IN THE MORNING, EACH SAILOR 
GETS 1023 

** DONE ** 
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Super Accuracy 



Under normal circumstances, your computer performs computations to six 
or seven digits of accuracy. Double precision computations increase accuracy 
to 13 digits or so. 

However, it is possible to do computations one digit at a time and assign 
each digit to an element in an array. This will achieve virtually any desired 
accuracy. Any, up to the maximum array size that is. 

This example program performs the rather simple operation of successively 
doubling a number (which is the same as raising 2 to a power). 

If the number to be represented is 8192, then: 



To add this to itself, first the rightmost digits are added: A(l) + A(l). 
If there is a carry, the variable C is set equal to 1, otherwise C is 0. The total 
is put into B in Line 100. 

If B is less than 10, there is no carry (C = 0) and the new A(l) equals B. If 
B is greater than 10 there is a carry (C= 1) and the new A(l) equals B-10. 

This operation is continued for all the digits (D) of the number and, when 
it is finished, the new number is printed in Lines 220-240. 

If A(N) is printed followed by a semicolon(;) for tight packing, the Basic 
print routine would leave a space in front of each digit (for the sign) and a 
space after each digit (for readability). In the program here, these spaces are 
not wanted, hence the print routine in Line 230 is used which prints the string 
value of the ASCII value of each digit (which is the same as the digit itself) 
but without the spaces. 

This program, incidentally solves the challenge to calculate the number of 
moves in the Towers of Brahma problem (see "Change for Any Amount to 
$5.00). The approach is also used in the next section, "Palindromes." 



A(4) 
8 



A(3) 



A(2) 
9 



A(l) 



2 




PRINT '^COMPUTES 2 TD Nth 
RNY" : "RCCURRCY:": : 
DIM fi < 1 > 
M = 



50 



>RUN 

CgMPUTES 2 TO Nth TO RNY 
RCCURflCY . 



160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 



I 2 

? i 

3 8 

4 16 

5 32 

6 64 

7 128 
S 256 

9 512 

10 1024 

II 2048 

12 4096 

13 8192 

14 16384 

15 32768 

16 65536^ 

17 131072 

18 262144 

19 524288 

20 104S576 

21 2097152 

22 4194304 

23 3383608 

26 67108864 

27 134217728 

28 268435456 

29 536870912 

30 1073741824 

31 2147483648 

32 4294967296 

33 S589934592 

34 17179869184 

35 34359738368 

36 68719476736 

37 137438953472 



!16 

^32 



C = 

IF 
IF 



> = B 

KB 7HEN 

D=I 

FOR N=:I TO 
PRINT CHR« 
NEXT N 
PRINT 
GOTD 70 



100 
100 



I STEP -1 
C fl < N ) + 4 S > 5 



^^^5^599627 3 7 4 9^ 
53 900719925474tD992 
- ■ 18014398509481984 
360287970 18963968 
72057594037927936 
1441 1518S075855S72 
28823037615171 1744 
5764607 523034 23488 
1 152921504606346976 
2305843009 213693952 
461 1686018427387904 
9223372036854775808 . 
1344 6744073709 551616 
3 6893488 1474 19103 23 2 
73 78 6 976294838206464„ 
1 4 7573 9 525 89 6 7641 29 2 8 
2951479051793 52825856 
5902958103 587056 51712 
1 180 59 16207 1741 130342 4 
23 61 18 324 143 482260 6 84S 
4722366 482869645213696 
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56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
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Palindromes 



A palindrome is a word, verse, or number that reads the same backwards 
or forwards. For example, the words *'mom" and "eye" are palindromes. So 
are each of the lines in this verse: 

Egad, a base life defiles a bad age 

Doom an evil deed, liven a mood 

Harass sensuousness, Sarah 

Golf; No, sir, prefer prison-flog 

Ban campus motto, "Bottoms up, MacNab" 

Numeric palindromes are those numbers which read the same backward as 
forward. The examination of these numbers is a field rich with possibilities 
for creative computing. 

One conjecture concerning palindromes raises an interesting unanswered 
question. Begin with any positive integer. If it is not a palindrome, reverse its 
digits and add the two numbers. If the sum is not a palindrome, treat it as the 
original number and continue. The process stops when a palindrome is ob- 
tained. For example, beginning with 78: 

78 
+ 87 



165 
+ 561 

726 
+ 627 

1353 
+ 3531 

4884 

The conjecture, often assumed true, is that this process will always lead to 
a palindrome. And indeed that is just what usually happens. Most numbers 
less than 10,000 will produce a palindrome in less than 24 additions. But 
there is a real thorn in the side of this conjecture, the number 196. Can you 
determine if a palindrome will ever be produced with a starting number of 
196? 

The number 196 will produce 1675 after two reversals, but after 100 rever- 
sals the resultant sum has 47 digits and is still not palindromic. Why mention 
1675? Because ten other numbers under 1000 will also lead to the sum of 
1675 and thus may not become paHndromic. The first five of these numbers 
are 196, 295, 394, 493, and 592. What are the other five? 
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The program here will accept any number as a starting value and complete 
the process of adding the successive reversals and testing if the sum is a palin- 
drome. Try it with some numbers and see if you can identify any patterns. 



10 PRINT "PRQGRflfl TESTS C□^^J 
E C T U R E " 

20 PRINT "THAT REVERSING DIG 
ITS ft M n " 

30 PRINT "flDDING WILL fiLWflYS 

PROKUCE" 
40 PRINT "fl PiflL I NDRDME . " 
45 PRINT "ENTER CZERQ) TO 
QUIT " 

5 HIM B < 5 > 

6 PRINT 

70 IWPIJT " START NUMBER =":fl 
75 IF fl=0 THEN 440 

7 7 PRINT 

80 E = g 

90 E=E+1 

10 fl = fl .'• 1 

110 IF INTCfl)>0 THEN 90 

120 FOR C=E TD 1 STEP -1 

130 fl=fi*10 

140 B(;C> = INT(:fl-10*INT<fl,'10> > 

150 NEXT C 

16 D = 

170 FDR C=l TD INT<E/'2> 

180 IF B •:: C > =E «: E+ 1 -C ) THEN 200 

190 D=l 

200 NEXT 

210 FDR C=E TO 1 STEP -1 

220 PRINT CHRS < E < C > +48 .^^ 5 

230 NEXT C 

240 IF 11=1 THEN 270 

250 PRINT " IS fl PflLINBRDME" 

260 GOTO 60 

270 PRINT " IS NOT fl PRLINDR 
□ ME" 

2 8 IF E /' 2 < = I N T C E 2 > T H E N 3 

290 BC INT<;E-'2> + 1 )=2*B<: IWT^E." 

ioO^FDR C = l T0„INT<E/;2> 
310 E<:C> = BCC>+BCE+1-C> 
320 NEXT C 

330 FDR C=l TO INT<:E.'2> 
340 E c: E+ 1 -C> s=B < C> 
350 NEXT C 

360 b<:e + i> = o 

370 FOR C=l TO E 

38 B <: C + 1 ) = B C C + 1 .> + I N T < B < C ) I 
> 

39 B <; C > = B < C ) - 1 * I N T < B 'C C ) . '• 1 
4 00 NEXT C 

410 IF^B(:E+1><=0 THEN 160 

420 E=E+1 

430 GOTO 160 

4 4 END 
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>RUN 

ENTER ^iZERD) TD i3U I T 

START NUMBER =19 

^?«I?^"DT^fl PALINDROME 
lip IS NOT fl PALINDROME 
121 IS fl PALINDROME 

START NUMBER =96 

9 6IS^NDT^fl„PALINDR0ME 
lb5 IS NOT A PALINDROME 
?2b IS^NQT A PALINDROME 

i§ PRLIN DROME 

4y84 IS A PALINDROME 

START NUMBER =0 

Using this method, write a program that examines all the integers between 
1 and 10,000 excluding those that sum to 1675 at any point. What does this 
show? By the way, you will have to devise a way to deal with 14-digit integers 
which are larger than your computer can normally handle. 

Huh? Is this program "too hot to hoot?" 
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4 

Convergence and Recursion 



The computer is especially suitable for doing repetitive and tedious calcula- 
tions. Two mathematical approaches for solving problems that involve repet- 
itive calculations are convergence and recursion. 

Some problems can be reasonably easily stated in words or described with a 
few simple equations but there are many possible solutions. For example, how 
many ways can you make change for a dime? It is simply stated and the num- 
ber of ways can be enumerated fairly easily: two nickels, one nickel and five 
pennies, or ten pennies, three ways in all. But if you want to solve for all the 
ways of making change for a dollar or five dollars, it would be nice to have 
some help. 

Help on this kind of problem comes from a class of computer program that 
simply breaks the problem into smaller ones and counts up all the alternative 
solutions according to a set of rules. But an even more powerful technique is 
known as recursion. Using this technique, a simple solving algorithm or rou- 
tine is set up to solve the smallest subset of the problem. The unique power in 
a recursive routine comes from the ability of the routine actualUy being able 
to call itself. This is discussed further in the second program in this section. 

Another approach for solving problems that do not have an exact answer is 
that of successive approximations. For example, the exact value of pi, e or the 
length of an irregular curve cannot be precisely determined. But by means of 
increasingly accurate approximations, it is possible to approach the desired 
value from above or below or to converge on it from two directions. The last 
four programs in this chapter illustrate successive approximations and 
convergence. 
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Change For a Dollar 



Even though there is very Httle you can buy for a penny these days, the 
coin will probably be around for some time to come since it is needed to make 
change for odd amounts of sales tax and to fill up penny collections. 

Today U.S. coinage consists of five coins: penny, nickel, dime, quarter, and 
half dollar. How many ways can coins of these denominations be used to 
make change for one dollar? For example, one way is two half dollars, an- 
other is one half dollar and two quarters, and so on. Make a best guess now 
and write it down before you read further. 

There are several different ways to approach a problem of this kind. One is 
to break it down into smaller, more easily solved problems. In other words, 
how many ways can you make change for a quarter? For a dime? You would 
solve these subproblems and combine the answers to give the overall solution. 

If you were more mathematically inclined, you could write a series of equa- 
tions relating each piece of change to every other one and to the dollar and 
solve them. 

A third approach is to do the problem by writing down combinations until 
all the different possibilities are exhausted (or until you are exhausted) and 
then count them all up. This might be called solving the problem by exhaus- 
tion and is a method quite suitable for putting on the computer. 

Write a program that uses this approach to solve the problem. If you use 
loops and count by one, it could take a long time for the computer to run 
through all the possible combinations, possibly many hours. 

Also, if you want to print out all the possible combinations, be warned that 
the printing could also take quite some time and a fair amount of paper. 
There are more combinations than you might think! 

In fact, most people will not be able to guess the answer to this problem, or 
even come close. Ask several of your friends how many ways they think a 
dollar can be changed. Record all the responses and then tabulate them on 
your computer. What is the mean (average) of all the guesses? The extremes? 

The program included here uses the first method discussed to solve the 
problem, in particular, breaking down the problem into subproblems and 
then combining the solutions into one final answer. 

First, the main problem is broken into the next smaller one of making 
change for half dollars. There are three such problems: no half dollars (H = 
0), one half dollar (H = 1), and two half dollars (H = 2). The last problem is 
trivial since there is only one way, but the other two need to be broken down 
further. 

This is done by dividing the remaining money into quarters and consid- 
ering the subproblems on down to the lower denominations. As the number 
of subproblems is expanded, each one becomes easier to solve. In fact, 
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subgoals, which can be solved in only one way, are finally reached. For exam- 
ple, if H = 1, Q = 1, D = 2, and N = 0, then the pennies (P) must equal 
five in order that the total equal 100. 

Notice that at the quarter, dime and nickel stages, adjustments are made in 
the limits of the loops depending upon how much money there is left to 
change. For example, if H = 1, the only possible subgoals for quarters are 0, 
1, and 2, but not 3 or 4. Also notice that there is no need to test combinations 
of coins to see if they add up to 100, nor is it necessary to include the penny 
as a variable. Simply counting the number of subgoals is sufficient since each 
one can be solved in only one way. 



C = 
FOR 
FOR 
FDR 
FOR 



H = 
Q = 
D=0 
N = 



TQ 
TO 
TO 
TO 



4-2*H 

1 0-5*H- . 5*Q 
20-1 0*H-5*Q-2* 



C=C+1 

NEXT N 

NEXT D 

NEXT Q 

00 NEXT H 

10 PRINT "THERE 

ERENT WfiVS":"TD 

LflR BILL. " 



ARE" S C! 
CHANGE 



"DIF 
A DO 



>RUN 
THERE ARE 
TD CHANGE 



292 DIFFERENT 
A DDLLAR BILL 



WAYS 



** DONE ** 



Try to make some changes in this program or write a new one to solve the 
following problems. Say you want two quarters in your change to play some 
video games. In how many ways can a dollar be changed to provide at least 
two quarters? 

Visiting a small town, you find the parking meters still take pennies. In 
how many ways might you get change so that you had at least three pennies? 
Is this any different than the number of ways that would give you five pen- 
nies? Say you want to make a phone call also; in how many ways can you 
change a dollar to produce at least four pennies and one dime? 
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Change For Any Amount to $5.00 



Another way to attack the change problem in the previous section is by 
means of the programming technique called recursion. Get familiar with this 
one— it is very powerful! Donald Piele and Larry Wood described this 
method in an issue of Creative Computing. 

First, define the variables which represent the number of ways to make 
change for n cents using the coins specified: 
A Only pennies 
B Nickels and pennies 
C Dimes, nickels, and pennies 
D Quarters, dimes, nickels, and pennies 
E Halves, quarters, dimes, nickels, and pennies. 
Initially, there are two subproblems in making change for n cents. In the 
first, no half dollars are used, and D is the number of ways to change n cents. 
Second, when one or more halves are used, after one is paid, there remain 
n-50 cents to pay which can be done in En_5o ways. 

Since these two cases are mutually exclusive, it can be inferred that £„ = 
Dn + En_5o- Similarly, 

Dn = C„ + Dn_25 
= Bn + Cn_io 
Bn = An + Bn_5 

Now, begin with the simplest case and build up to Eioo- First of all, it is 
easy to understand why E^ = 1. From above, when n = 50, E^^ = D^^ + 
Eq, and it is possible to make change for 50 cents only one more way if half 
dollars are allowed. Therefore E^ = 1. Likewise, D^ = C^, = B^ =Ao = 1. 
It is also true that An = 1 for all values of n since there is only one way to 
make change using only pennies. Now the recursive relationships can be used 
to solve the original problem. 

This is the strategy used in the program. It also has the added advantage 
that it can count the number of ways of making change (with coins) for any 
specified amount. 

Now it is your turn. Can you modify the program here to include dollar 
bills so it could count the number of ways to make change for any amount up 
to $10.00? 

Using any method of change making you prefer, write a program to make 
change for one ruble. Russian coins come in denominations of 1, 2, 3, 5, 10, 
15, 20, 50, and 100 kopecks. There are 100 kopecks in one ruble. 
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10 PRINT "PRDGRnn CDMPUTES W 
RVS TD"!"MfiKE CHRHSE FDR RNY 

RtlDUNT" s "UP TD $5-00" 
2 PRINT 

30 DIN fi •■: 1 1 > 

40 DIN Be: 101) 

5 DIN C < 1 1 > 

6 DIN D < 1 1 > 

7 DIM E < 1 1 > 

80 INPUT "RNDUNT TO CHRNGE = 

90 'n= I MT < 20*X ;> + 1 

10 n < 1 > = i 

110 B <: 1 > = 1 

120 c t: 1 > = i 

13 D <: 1 ) = 1 

14 E •:: 1 :> = 1 

150 FDR .J = 2 TU M 

16 fl J ) = 1 

17 b c j ) = r c j > + b < j - 1 > 
180 C(:j)=bcj> 

190 IF J<=2 THEN 210 

200 C < J > =E C J rj +C < J-2 > 

2 10 D C J :) = C C J > 

22 IF JOS THEN 240 

2 3 D <: J > = C < J J + D < J - 5 > 

240 eCJ>=D<J> 

250 IF J<: = 10 THEN 270 

260 EC J>=D< J>+E< J-10> 

270 NEXT J 

280 PRINT 

290 PRINT "YOU CJ=lN MRKE CHAN 

»3E FDR J: " ; X S " I N " J E < t1 > 5 " D I F 

FERENT UlflVS" 

300 STOP 



> R U N 

PRDGRflM COMPUTES MAYS TD 
MAKE CHANGE FDR ANY AMOUNT 
UP TO S5,00 

flMQUNT TO CHRNGE =4.99 

YDU CAN riAKE CHANGE FDR 
f 4.99 IN 57420 
DIFFERENT WAYS 

* * DONE * * 

> R U H 

PROGRAM COMPUTES WAYS TO 
MAKE CHANGE FDR ANY AMOUNT 
UP TD :E5.00 

AMOUNT TO CHANGE =1.13 

YDU CAN MAKE CHANGE FDR 

ff 1.13 IN 384 DIFFERENT WAYS 

* * DONE * * 
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Perhaps the most famous problem used to demonstrate the principles of 
recursion is the Towers of Brahma. It is sometimes called the Towers of Ha- 
noi or Pharoah's Needles. Here is the problem in as close to original form as 
possible. You should be able to solve it with a relatively short program using 
recursion. 

In the great temple at Benares beneath the dome which marks the center of 
the world rests a brass plate in which are fixed three diamond needles, each a 
cubit high and as thick as the body of a bee. On one of these needles, at the 
Creation, God placed 64 discs of pure gold, the largest disc resting on the 
brass plate and the others getting smaller and smaller up to the top one. This 
is the Tower of Brahma. 

Day and night unceasingly, the priests transfer the discs from one needle to 
another, according to the fixed and immutable laws of Brahma. These laws 
require that the priest on duty must not move more than one disc at a time 
and that he must place this disc on a needle so there is no smaller disc below 
it. When the 64 discs shall have been thus transferred from the needle which, 
at the Creation, God placed them, to one of the other needles; tower, temple, 
and Brahmans alike will crumble into dust, and with a thunderclap, the 
world will vanish. 

If the priests were to effect one transfer every second, and work 24 hours 
per day for each day of the year, it would take them 58,454,204,609 decades 
plus slightly more than six years to perform the feat, assuming they never 
made a mistake — for one small slip would undo all the work. 

How many transfers are required to fulfill the prophecy? Try out your pro- 
gram with fewer discs than 64 to make sure you are on the right track. Here 
is a table of the first few transfers: 



Discs Moves Discs Moves 

1 1 6 63 

2 3 7 127 

3 7 8 255 

4 15 9 511 

5 31 10 1023 
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Converge on e and Pi 



An incredibly important mathematical constant is designated by the small 
letter e. This constant is both irrational and transcendental. Look up those 
terms in a dictionary or math book if you wish, or just plunge on to the next 
paragraph. 

The constant e was first derived by John Napier, also the inventor of loga- 
rithms, to whom we owe an eternal debt of gratitude. Why? If e had never 
been discovered, advances in mathematics, physics, and astronomy would 
have lagged a century or more, because e is the base of all natural logarithms 
and these logarithms are the basis for many branches of science and 
mathematics. 

How is e calculated? The constant e is the limiting value of this expression 
as n approaches infinity: 

e = (l + yn)i^ 

Its exact value can never be found, but to 15 places e equals 
2.718281828459045 . . . How can e be calculated? First take 1-% and square 
it; that equals I-Va- Then cube l-Va and you get 2.3686. Raising \-% to the 
fourth power gives 2.414, and so on. Write a program for this method and 
have it print out the initial value of e and each value after each additional 
fraction is added on. 

Another approach is to expand the expression above using the binominal 
theorem and, again, letting n approach infinity. The expression for the expan- 
sion is: 

^ ,1,1,1.1 JL 
^ = ^ ^ 27 3l T». 5l ' ' ' n! 

Now here is where the computer can again be of some assistance. Since 3! is 
3*2! and 4! is 4*3!, all the calculations need not be done for each additional 
fraction. Look at the program and particularly note the calculations in State- 
ments 70 and 90. 

10 PRINT '^CONVERGE ON E":"HI 
T 'FCTN^M^ TO HALT'* 
20 E=l 2 ^ 

3 1=0 2.5 
40 1=1+1 2-666S66667 
50 D=l , 2-708333333 

60 FOR J=l TO I 2-716666667 
70 D=D*J 2.718055556 
80 NEXT J 2.718253968 
90 E = E+1/If 2.71827877 
100 PRINT E 2.718281526 
110 GOTO 40 2.718281S01 

2. 718281826 
2. 718281S28 
3. 718281828 

61 



Pi is another important mathematical constant that is irrational (meaning 
its exact value can never be determined) and transcendental (meaning it is not 
the solution to any algebraic equation). Interestingly, pi was known and used 
by the ancients. Archimedes, who lived in the second century B.C., by using a 
regular polygon of 96 sides (nearly a circle), proved that the value of pi was 
less than 22/7 and greater than 3-^%i, a remarkable achievement for the 
mathematics of his day. 

Ptolemy in 150 A.D. used the value of 3.1416 for pi and in the middle of 
the sixteenth century the amazing fraction 355/113 was discovered, giving 
the value of pi accurately to six decimal places. 

Incidentally, in 1897, the General Assembly of Indiana passed a bill ruling 
that the value of pi was four. 

Several infinite series can be used to grind out increasingly accurate values 
for pi. One such series is + - % + %...). This series is called an 
arithmetic series and converges very slowly. The program here displays only 
every 500th value of the series. Don't be alarmed if the program does not 
seem to be running very fast; a fair amount of calculating is going on between 
each value that is printed. 

10 PRINT "CONVERGE □ N PI WIT 
tin^&";:UN^If'^I"rE SERIES^ 

20^print^^;;be prtient; this i 

SO^PRINT "HIT 'FCTN/4-' TO MR 
110 s=-s 

120 IF^Q<4 99 THEN 80 
1 30 Q = 

140 PRir^lT P*4 3.1435 9 666 

150 GOTO SO 3. 14059065 

3 , 1 4 2 6 6 5 6 
3 . 14 10 9 16 5 2 
3. 141993455 
3. 141258 6 52 
3 . 14 18 7 8 9 4 
3. 14 1342 153 
3. 14 1815321 
3. 141392253 
3. 141774836 
3- 14H^'5653 
3 . 14 17 4 6 8 8 
3. 14144 9 51 
3 . 1 4 17 2 6 2 5 4 
Actually, the approach used by Archimedes converges much more quickly 

and, with the aid of a computer, it is possible to go far beyond the 96-sided 

polygon used by Archimedes. 
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His approach was to construct inscribed and circumscribed polygons and 
measure the perimeters to approximate the circumference of a circle. Con- 
sider a polygon circumscribed around a circle of radius 1. The perimeter 
equals the length of one side times the number of sides. Since the tangent of x 
= AB/BC, but BC = 1, then tan (x) = AB and the length of a side = 2 tan 
(x). Since the circumference = 2 pi r and r = 1, then pi is the circumference 
(or perimeter of an n-sided polygon) divided by 2. 

Similar trigonometry leads to the perimeter of an inscribed polygon being 
equal to the number of sides times sin (x) * cos (x). 

The second program produces values for pi using inscribed and circum- 
scribed polygons. Unfortunately, there is one large flaw in the program, be- 
cause degrees must be converted into radians in Statement 50. This means, of 
course, that you must already know the value of pi, since the conversion fac- 
tor is 360 degrees divided by 2 pi. 

Setting this flaw aside, it is interesting to note how quickly this program 
converges on the value of pi compared to the preceding one. That is because 
this one converges geometrically rather than arithmetically. 

Can you figure out a geometric convergence to the value of pi that does not 
require that you know it (or a conversion factor) before you start? 

3 N = 6 

4 M=2*H 

5Q X = 360^<:N#57. 29578> 

f^oK*?.^. ^^T I " < :<> *cns < X ) , n 

70 GDTQ 40 

C □ N V E R G E □ H P I 

HIT •■ FC:TN.-^4 TD HFiLT 

INSCRIBED CIRCUMSCRIBED 

2.598076198 3.46410158 

i-?Slflll?I 3.215390281 

3.105828515 3.159 659915 



3.132628587 3.146086188 

3.139350176 3.142714573 

S- M 3. 141873023 

3.141452446 3.14166272 

3. 14155 7581 3.14161015 

3.141583865 3.141597008 

3.141590437 3.141593722 
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Convergence on Pi Revisited 



In answer to the question posed in the last paragraph of the preceding sec- 
tion, here is a way to converge on pi without knowing its value beforehand. 

As in the previous program, the basic approach is to add up the length of 
the sides on an inscribed polygon and divide by 2r to obtain a value for pi. 
The program starts with a square (four sides) and doubles the number of sides 
each time. If the old side length is S, then the length of S' of a side of a new 
polygon with twice as many sides is obtained by applying the Pythagorean 
theorem. In particular, 

-h (S/2)2 = R2 

(R-xy + (s/iy = (sy 

Thus, 

s' = V (r-Vr2.(s/2)2)2 + (s/iy 

It is easy to reduce this formula algebraically, but accuracy suffers if this is 
done. Also, you will find that S*S is slightly more accurate than St 2. 

Unfortunately, striving for maximum accuracy is somewhat moot on a 
computer that does not have double precision arithmetic. Notice that ac- 
curacy does not improve with more than8192sides and, indeed, as numbers in 
the calculations start to exceed the capacity of the computer, the accuracy 
starts to deteriorate badly. 

There is yet another method to compute pi by convergence. It uses discov- 
eries of Gregory and Euler. Gregory discovered the formula for arctangent: 
Arctan x = x - xV3 + xV5 - xV7 + ... 

Euler came up with a rather interesting formula for pi: 
pi = 4 (arctan CA ) + arctan (Vg)) 

See if you can combine these two formulas to calculate pi. If you are very 
clever, you can do your calculation to yield far more than the seven decimal 
place accuracy obtained by the programs presented here so far. 




10 PRIMT "CONVERGE CN PI »a) I T 
H" : " IHSCR lEED POLVeOHS":" " 
20 PRINT "WHEN PRDGRflM BREAK 
S, TYPE " s ' CDN ' TD CONTINUE" 

■ II II 

30 PRINT "SIDES" r "PERIMETER" 

40 R=10 

50 Z = 2 

6 N=2+Z 

70 S=R*S0R<Z) 

80 FOR K=l TD 28 

90 IF K014 THEN 110 

100 BREAK 110 

110 PRINT N,(:S*N>^'t:Z*R> 

120 V = S*S''<Z + Z> 

130 X=:R-SeR<R*R-Y> 

140 S=SQR<X*X+Y> 

150 N=N*Z 

160 NEXT K 



> R LI H 

CONVERGE QN PI WITH 
INSCRIBED POL Y 6 DNS 



W HEN 
'CON 



PRDGRflM BREAKS: 
TO CONTINUE 



TYPE 



SI DES 
4 
8 

16 

32 

64 

128 

256 

512 

1 024 

2048 

4096 

8192 

16384 



PERIMETER 
2 . 828427 125 
3. 061467459 
3. 121445152 

13 6548491 

14 3 3 115 7 
1 4 1 27725 1 
14 1513801 
14 1 57294 
141587725 
1415S1422 
14 1592 346 
141592577 
141592634 



3, 



3, 
3 . 
3 , 



* BREflKPDINT FlT 100 



CON 
32768 
6 5 536 



1310 7 2 
2 62 1 44 
5242S8 
1048576 
2097152 
41 94304 
8388608 
167772 1 
3355443L 
67108864 
1 34217728 
268435456 
536870912 

#* DONE ** 



6 
12 



1 41592649 
14 1592652 
141592 653 
141592654 
141592654 
141592654 
1 41592654 
14 15 9 2 6 5 4 
1 4 1 592654 
1 4 1592654 
141592654 
14 1592654 
141592654 
1 41592654 
141592654 
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Length of Any Curve 



The previous programs have demonstrated how it is possible to compute a 
very accurate value of pi by adding together the length of the sides of a poly- 
gon as it approaches a circle and dividing by 2r. 

Using a similar approach, it should be possible to inscribe a polygon, or 
portions of a polygon, inside any regular curve and thus determine the length 
of the curve. This program approximates the length of any curve as defined in 
Statement 110 by dividing it into an increasing number of subintervals and 
computing the sum of the secants (a straight line that cuts a curve at two or 
more points). 

To run the program, you must enter the formula or equation describing 
your curve in Statement 110 in the form: 

DEF FNA(X) = Any function of X 

You then type RUN and enter the end points of the curve you want to use 
in your calculation. These points are entered in the form of the abscissa, 
which means the horizontal (or y) coordinate of the point. 

The program is written to sum the successive secant lengths and to cal- 
culate the percent of change in each summation compared to the preceding 
one. The sample run uses the function 2x^ + -2x +3. Note the substan- 
tial improvements in the length calculation as the number of intervals in- 
creases from 2 to 16, but the rather slight improvements beyond 16. 

Try this program with different curves and functions. It might help to plot 
the function first (remember the program to do that?) and then compute its 
length. Is this method of length calculation more accurate for a function with 
no changes in direction of the curve within the interval selected or for a func- 
tion with one or more changes? Why? 



66 



iO PRINT "CURVE LENSTHS"!" •' 

20 PRINT "DEFINE YDUR FUNCTI 

□ W I N " ! " L I N E 110":! 

30 PRINT "ENTER flBSCISSflS OF 



PTS 



< SMALLER ONE FIRS 



P 7 Q 



NO 



OF SUM OF 

SECRNT 

LENGTHS 



•2-2 



END 
T ) " J 

40 INPUT 
50 PRINT 

PERCENT" 
60 PRINT "SUBIN- 
CHflNGE IN- 
FO PRINT "TERVfiLS 
LENGTH" 
80 PRINT 

90 S1=0 „ , 

100 FOR M=l TO 9 
110 DEF FNfl<X> = 2*X-"-3 + 3*X' 
*X + 3 
120 E = 2'^t:N-l) 
130 H=<G!-P>/'E 
■ - S=0 

FOR 1=0 TD E-1 

L'SQRC <:FNfl<P + I*H + H>-FNR< 

P+ I *H > > •^2 + H*H ) 
170 S=S+L 
180 NEXT I 

190 IF S1>0 THEN 220 
200 PRINT E5 TlFIBCS^ ;S5 TfiB< 18? 
; "NO PREV VfiL" 
210 GOTO 240 

— P5=< <:fiBs<si-s;o.'Si >*ioo 

PRINT E5 TlF!B<5;' ? S? TflBC 17> 



140 

150 
160 



22 

23 
5 P5 
240 
250 



S1 = S 
NEXT 



CURVE LENGTHS 

DEFINE VQUR FUNCTION IN 
LINE 110 

ENTER flBSCISSftS OF END PTS 
<SMflLLER ONE FIRST>? -1,6 



NO . DF 

S U B I N - 
TERVflLS 



1 

2 
4 
8 

16 
32 
64 



525, 
525 
529. 
531 , 

531 

532, 
5 32 , 



SUM DF 
SECANT 
LENGTHS 

0466646 
158262 
652247 
0171331 

9642557 

1 65657 

0416782 



PERCENT 
CHRNGE IN 
LENGTH 

NO PREV VRL 
. 0212547633 
. 855739171S 

. 2576947669 

. 1783600797 

. 0098333641 

. 0047202447 
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Converge on a Square Root 



Since most square roots are irrational, methods used to calculate them usu- 
ally involve successive approximations. Although you can simply call up the 
square root function in Basic or on many pocket calculators, it is interesting 
to explore various methods of calculating square roots without these built-in 
functions. After all, these built-in functions are nothing more than successive 
approximation routines already installed in the machine. 

Obviously, a square root is the inverse of the operation of squaring a num- 
ber. All of the methods of calculating square roots use this fact, but the way 
in which it is employed is quite different in various calculators and 
computers. 

The program here calculates an upper and lower limit for the square root 
of a number and, by successive approximations, pinches the root to within a 
smaller and smaller interval until it reaches the desired level of accuracy. 

The starting value for the lower limit is and for the upper limit the num- 
ber, Z, whose square root is sought. The program then divides this interval 
into ten steps by simply dividing the difference between the numbers by 10. 
The variable I is increased from the lower limit to the upper one by the value 
of the step, S. At any point, if I squared becomes greater than Z, a new upper 
limit is set to I and a new lower limit is set to I - S. 

This method converges very quickly and adds approximately one decimal 
place of accuracy with each pass beyond the third. What happens when you 
enter into the program a number that has an exact square root such as 25 or 
49? Why? 

Another approach to calculating square roots by successive approxima- 
tions is to start with a trial root, X. If X * X is less than the original number 
N, then increase the trial value by a 0. 1. If X * X is greater than N, return to 
the previous value. This is the first digit of the root. Now, start advancing by 

0. 01. Continuing in this way, one digit is developed at a time until the desired 
precision is reached. 

This method is quite suitable and fast for square roots of numbers less than 

1. A good first trial root value is 0.1. But is it suitable for larger numbers? 
How should a starting trial value be determined? In this method, especially 
for numbers greater than 10, the initial trial value for the root matters a great 
deal in determining the length of time it will take for the calculation to con- 
verge. Write a program using this method and compare the speed and ac- 
curacy with the program in the book. 
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10 PRINT "CflLCULf=lTE SQURRE R 
□ D T S " : " " 

20 PRINT "SET flCCURflCY IH LI 
N E 4 " 

30 INPUT "YOUR NUMBER ="!2 
40 £=.0 0001 

50 PRINT ' 'LOWER LIMIT"." 

UPPER LIMIT" 
60 fl = 
70 B = Z 

SO S=<B-fl>'''10 

lOo'^IF SBS<fl*B-Z><E THEN 190 

110 FOR I=fi TO B STEP S 

120 IF Z<I*I THEN 160 

130 NEXT I 

140 B=B*10„ 

150 GOTD 80 

160 B=I 

170 fl=I-S 

180 SOTO SO 

190 PRINT " "i"flCCURftCY LIMI 
T PEflTHED" 

200 PRINT "RVERflGE VALUE 
C ft + B > ^2 

> R U N 

CflLCULflTE SQUARE ROOTS 

SET flCCURfiCY IN LINE 40 
YOUR NUMBER =3 

LOWER LIMIT UPPER LIMIT 

}:Pi :?4 

RCCURftCY LIMIT RERCHEEi 
ftVERRGE VALUE = 1.73205 

** DONE ** 

'*CflLCULRTE SQURRE ROOTS 

SET RCCURflCY IN LINE 40 
YOUR NUMBER =54 

LOWER LIMIT UPPER LIMIT 
54 
5.4 10. Q 

7.344 

7. 344 7. 3494, 

7 34832 7.34886. 

7 348428 Z-§f|^i?o 

7.3484658 S-§484712 

7.3 4846904 7. 3484 6 958 

RCCURflCY LIMIT RERCHEIi 
flVERRGE VALUE = 7.3484fa931 

** DONE ** 
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5 

Compounding 



As with successive approximations and recursion, compounding requires 
many repetitive calculations. Some compound interest and growth situations 
can be represented by a formula, but for many problems, solving by repetitive 
calculations is an excellent approach. 

There is nothing magical about compounding; you generally start with an 
initial amount of money, number of animals, etc. This quantity then grows or 
diminishes by a certain percentage at set intervals. The new amount is the old 
amount increased or decreased by this percentage. Repeat this calculation 
over and over again, and you have solved the problem. 

Five different types of problems which involve compounding of some sort 
are described in this chapter. Try the extra problems that follow some of the 
programs; you may be surprised at the results. 
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Indians and Interest 



Here is a simple compound interest problem that produces an astonishing 
answer. The problem has to do with the sale of Manhattan Island to the 
Dutch for approximately $24 worth of trinkets and beads. 

If the $24 that the Indians received in 1626 had been deposited in a bank 
paying 5-y4% interest compounded annually, how much would it amount to 
in 1983? 

The program here solves this little problem by making use of the formula 
to calculate compound interest. In particular, if P dollars are invested at an 
interest rate of R (expressed as a decimal) and compounded N times, then the 
total amount A is given by the formula: 

A = P(l + R)N 

How much was gained in 1983 alone? How much was gained in the decade 
from 1974 to 1983? You can change the value of N in Statement 30 to get the 
answers to these questions. But a better way might be to enter the ending year 
with an input statement. 

Can you read a number expressed in the E (exponential format)? In a more 
conventional format, the number is $11,176,500,000 or $11.2 billion. 

There are many other ways you can improve the program and make it 
more suitable for general purpose compound interest calculations. Modify it 
to accept any starting and ending year, any rate of interest, and any starting 
principal amount. 

The problem as stated is somewhat unrealistic since banks were not paying 
interest rates of 5-y4% from 1626 to 1983. Change the program to calculate 
the total amount based on the following interest rates: 



Years 


Interest Rate 


1626-1830 


1.5% 


1831-1870 


2.0% 


1871-1910 


3.0% 


1911-1921 


3.5% 


1922-1929 


6.5% 


1930-1940 


2.3% 


1941-1945 


3.5% 


1946-1960 


5.3% 


1961-1980 


6.5% 


1981-1983 


9.5% 
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10 P=24 

20 R = . 05 75 

30 N=19837l626 

P,9 ?KI^'T "IF THE S24 WHICH T 

RD BEEN""""' "H'^TTflN ISLAND H 

60 PRINT "INVESTED RT 5 



>RUN 

IE„JHE S24 WHICH THE INDIfiMS 
RECEIVED IN 1626 FOR MflN- 
HfiJTflN ISLAND HftD BEEN 
INVESTED RT. 5. 75^., Iv 1983, 
ri'.'?y^?5g55fg RMDUNTED TO 

HOME 




Systematic Savings 

In the previous program, compound interest was calculated by means of a 
formula well known to bankers, money lenders and real estate agents. How- 
ever, if you did not know the formula, how would you calculate interest on 
money in a savings account by hand or with a calculator? 

You would probably begin by multiplying the principal amount P by the 
interest rate R and adding that amount to the original principal at the start of 
the second year. Doing that for all the years the money is in the bank will 
yield a final amount. Why not write a program to perform the calculations in 
this manner rather than use a formula? Here is such a program. 

It is set up for an initial principal amount of 100 (Statement 50), an interest 
rate of 10% (R = 0.1 in Statement 60) and ten years (N = 10 in Statement 
70). Naturally these values could be read in using input statements. The 
clever calculation in Statement 100 rounds off the amount to two decimal 
places (dollars and cents). 

In contrast with the program in the previous section, this one does not use 
a compound interest formula, but simply adds the interest each year to the 
growing principal amount in Statement 90. 



10 PRINT "CRLCULFITES INTERES 
T DN $100" s " INVEST EB RT 10^: 
FUR 10 Y R S " 5 " " 

20 PRINT '*flT END" , "CURRENT" 
30 PRINT "DP YERR*' . " BRLRNCE" 
40 PRINT 
50 P=100 
60 R=. 1 

7 N = 1 

8 FOR I=i TD H 
90 P=P+P*R 

iOO PRINT I , INT(:P*100+. S^-MO 


110 NEXT I 



"CRLCULRTES INTEREST ON SIQO 
INVESTED RT 10^: FOR 10 YRS 



RT 
□ F 



YEAR 



1 
2 
3 

i 

6 
7 

I 

10 



CURRENT 
BflLfllNCE 



110 
1 2 1 
133 
146 
161 
177 

1 94 

2 1 4 
235 
259 



. 1 
. 4 1 
. 05 
. 1 6 
. 87 
. 36 
. 79 
. 37 
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Now, how could this program be modified to allow for a plan of systematic 
saving? In other words, instead of letting the $100 lie around all lonely while 
it is compounding, each year you add another $100 to it. With this program, 
making the modification for systematic savings is easy: Statement 105 is 
added to add the new deposit each year to the ever growing principal. 

EVERY VEfiR":"FnR ID YEARS"!" 
gg^PglNT "RT END c 

iLflNCE""^ "OF VEFtR INVESTED E 

4 PRINT 

5 P=lCiO 
55 C=100 
60 R= . 1 

7 N=li:i 

80 FDR_^I = i TO N 

105P=P+C 
.110 NEXT I 

> 

> R U H 

fstiyvgs^^T^ys-^'^i?jR?M^g^ 

FOR 10 YEARS "^'tni rtMK 

□ P VpSp TMUir.-x.rn CURRENT 
UP- fERR INVENTED BRLRHCE 

T i'rS 231 

i 3 00 364.1 

i 400 SlO.Sl 

% 1 00 8 4 3. 72 
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In an effort to attract depositors, many banks over the last 20 years have 
started offering savings and investment accounts in which the deposits are 
compounded more often than annually. In the early 50's, many banks went to 
quarterly compounding; in the late 50's, to daily; and in the early 60's, some 
"competitive" S&L's went to continuous compounding. 

How much does more frequent compounding really mean? Modify either 
of the two programs to compute the interest for more frequent compounding. 
What is the difference in interest for one year for annual, quarterly, monthly, 
daily, and continuous (every second) compounding on a principal amount of 
$1000 invested at 8%? How about for ten years? 

Incidentally, if you want to use the formula method in the previous section 
for this calculation, the formula for P principal invested at R rate com- 
pounded N times per year is: 

A = P(l + R/N)N 

Try this problem which uses the same principles of compounding. Con- 
sumer prices rose an average of 8.8% in 1980. While the government keeps 
trying to bring inflation under control, they don't seem to be meeting with 
much success. Assuming that prices continue to go up this much (8.8%) ev- 
ery year, how much will a $6000 economy car (1980 dollars) cost in the year 
2000? How much will it cost when you are 65 years old? 
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Systematic Savings Revisited 



Using the formula for compound interest and systematic savings, this pro- 
gram will calculate the amount accumulated after a given period of time. 

When you run the program, it will ask how much you wish to save each 
month, the number of compounding periods in a year, the interest rate, and 
the length of time you wish to continue your systematic savings program. The 
program will then calculate the total amount at the end of that period. 

From the preceding programs, you should be able to see that a systematic 
savings program is a very effective way to accumulate a nest egg for the 
future. 

10 PRINT "CRLCULfiTES INTERES 

T RND ^' 2 " BRLRNCE FDR R SYSTEM 

R T I C " : " S R V I N G S P R □ G R R M " : " " 

20 INPUT ^'MONTHLY DEPOSIT =" 
• 

30 INPUT "CnMPDUHDINC PERIOD 
S Y R = " ! B 

40 INPUT "INTEREST RATE <XX. 
y. 'i> = " : R 
50 R«=R.'-I0O 

60 INPUT "NUMBER OF YERRS =" 
s N 

70 PRINT 

80 PRINT "RT END C 
URRENT " 

Ip^PRINT "OF YEAR INVESTED B 

R L R N C E 

100 PRINT 

110 C=12*R 

120 P=12*R 

130 FOR 1=1 TO N 

14 P = P * <: 1 + R / p. > B 

150 PRINT I?TflBc:il.');i*C!iTRBc; 

20> J INT<P*100+. 5>x'100 

lt)0 P = P + C 

170 NEXT I 



>RUN 

CRLCULRTES INTEREST RND 
BRLRNCE FOR fl SYSTEMRTIC 
SAVINGS PRDGRRM 

MONTHLY DEPOSIT =10 
CDMPnUNDINS PERIDDS-'YR =12 
INTEREST R R T E < X H . X J; ) =15 
NUMBER DF YEARS =2 

RT END CURRENT 
OF YEfiR INVESTED BRLRNCE 

1 120 139.29 

2 240 300.97 

** DONE ** 
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CflLCULflTES INTEREST^RND^ 
BRLRNCE FDR^R^SYSTEMRTIC 
SflV IN5S PRDGRfiM 

MONTHLY DEPOSIT =10 
CDMPDUNDING PERIODS/YR =1 ^ 
INTEREST RRTE^^XX.X > =3.5 
NUMBER OF YERRS =8 

RT END CURRENT 
□F YEftR INVESTED BRLRNCE 

1 120 13 0.2 

2 240 271.47 

3 360 424.74 

4 4 80 5 9 1.04 

5 600 771.48 

6 720 967.26 

7 840 1179.6 3 

8 960 1410.15 

* * DOME * * 

CRLCULRTES INTEREST RND^ 
BRLRNCE FDR R SYSTEMRTIC 
SRVINGS PROGRflM 

MONTHLY DEPOSIT =10 
CDMPDUNDING PERIDDS/YR =365 
INTEREST^RRTE^ <XJ<. X ) =8.5 
NUMBER DP YERRS =8 

RT END CURRENT 
DP YEAR INVESTED BRLRNCE 

1 120 130.64 

2 2 4 2 7 2.88 

3 3 60 427.7 3 

4 480 596 . 32 

5 6 00 7 7 9.36 

6 7 20 9 7 9.68 

7 840 1197.23 

8 960 1434.08 

** DONE 

Try to combine the things that you have learned from the programs in 
these sections to write a completely generalized program for systematic sav- 
ings. It should accept the following information as input: 

• Initial deposit 

• Frequency of periodic deposits 

• Amount of each deposit 

• Interest rate 

• Interest compounding frequency 

• Length of time of savings program 

Your program should produce output in tabular form showing years (or 
other time periods), amount invested without interest, total amount with in- 
terest, and the interest alone. 
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Loan Payments 



Although saving money in a systematic way is a noble goal, many people 
frequently find themselves talking to a different bank officer, namely the one 
in charge of loans. 

This program will calculate the payments for a loan for a period of one 
year or longer. The program asks you to enter the key facets of the loan in 
question: amount borrowed, annual rate of interest, interval between pay- 
ments and term of the loan in years. 

The sample run shows a relatively short-term loan (2 years) of $3000 at a 
bargain basement interest rate of 8%. The monthly payment of the loan is 
computed to be $135.68 and the total interest $256.34. Why is the total in- 
terest not $240 (8% of $3000)? Instead it is 8.54% Is this a mistake in the 
program? 

Try the program for some longer term loans at real world interest rates. 
For example, run it for an automobile loan of $8000 at 12% over a 5-year 
period. Perhaps you can see from this that systematically saving your money 
and paying cash for an item makes more sense than time payments. 

Run the program to calculate the mortgage payments on a $150,000 house 
with a $40,000 down payment. Try it with a 16% interest rate stretched over 
a 30-year period. Ouch! Look at all that interest! 

At the heart of this program are Statements 150-180. Can you see what is 
happening in these calculations? 



□ yTPLIT flS TRELE iEWTFP 
TDTflLS DNLY ENTER 2> 



OR 



JSTfiL INTEREST = 414 

AND TDTflLS S 539153. 
* * n □ N E * * 



1 5 3 

2 50 
797 
& 1 



= 1 7 

5 = ^ 
1 ) 
? 2 

. 61 


. 1 8 
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10 PRINT "PROGRnn CflLCULflTES 
PRVMENTS " 5 DN n LONG-TERM L 
UfiH " : " " 

20 INPUT "flMDUNT BQRRDWEB =" 
: R 

30 INPUT "INTEREST RATE <X>;. 

4 ■ I N P UT " P R V M E N T I N T E R V RL < 
M □ S . > = " : P 

5 INPUT ''TERM DF LDRN <YERR 

igi ^1 ^ II ■ y 

6b PRINT "DUTPUT flS TABLE <.E 

NTER 1>":"DR TQTflLS ONLY EMT 

E R 2 > " ; 

7 INPUT V 

80 PRINT 

90 IF V=2 THEN 150 

too PRINT " " ! " T fl B L E I N " 5 P ? " 

riDNTH I NTERV'RLS " : " " 

110 PRINT T fl B < 6 ) 5 " P R N C P L " 5 T fl 

B<22> 5 "PRNCPL" 

120 PRINT "PER- AT BES INT 
DUE REPAID" 

130 PRINT "IDD OF PER END 

PER END PER" 

140 PRINT 

150 2=t:Y*12)/P 

160 K=c: I#<P^12> ■J.'lOO 

17 E = A * K z' < 1 - 1 / < 1 + K > Z > 

18 E=INTt:E*100+.5).'100 
190 C=A 

200 F=0 

210 D1=0 

22 T1=0 

230 T1=T1+1 

240 IF Tl>2 THEN 360 

250 B=T1 

260 C=C-F 

270 D=C*K 

2 3 F=E-D 

290 C=INT<C*100+.5;'/10 
300 D=INTCD*100+.5)/'100 
310 F=INT(:F*100+. 5)''100 

3 20 D1=D1+D 

330 IF V=2 THEN 230 

?§?TS§?^Lfr^'''''^'^«^^^^^ 

350 GOTO 230 

360 IF y=l THEN 410 

370 Dl = IMT«CDl*100+.5>/100 

380 PRINT "TOTAL INTEREST =" 

390. PRINT "PRINCIPAL REPAID 

400 GOTO 430 

4 10 PRINT 

420^PRI NT. "TOTAL"; TAB<14>SB1 
;TABC22.>!!fl 

430 E5=INT<<Dl+A>*10 + .5)-'10 

4 4 E 6 = E 5 . ' < t: V * 1 2 ) ^ P ) 

450 E6=INTC100*E6+.5>/100 

4^0^PRINT "MONTHLY PAYMENT = 

470 PRINT "AND TOTALS S"5E5 
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>RUH 

PRDSRflM CRLCULFITES PAYMENTS 
DN R LONG-TERM LDflN 

RflDUNT BDRRQWED =3000 
INTEREST RATE CKX.X Ji > =S 
PAYMENT INTERVAL M S . > =1 
TERM OF LOAN CYEARS) =2 
OUTPUT AS TABLE CENTER 1) 
□R TOTALS nMLV ENTER 2>? 2 

TOTAL INTEREREST = 256.34 
PRINCIPAL REPAID = 3000 
MONTHLY PAYMENT =« 135.68 
AND TOTALS S 3256.34 

** DONE ** 

> R U N 

PROGRAM CALCULATES PAYMENTS 
DN A LONG-TERM LDflN 

AMOUNT BORROWED =8000 
INTEREST RATE CXX.X y. ) =13 
PAYMENT INTERVAL <MOS.> =3 
TERM OF LOAN CYEARS) =5 
OUTPUT as TABLE < ENTER 1> 
OR TOTALS ONLY ENTER 2:)? 1 



TABLE IN 3 MONTH INTERVALS 

PRNCPL PRNCPL 

PER- AT BEG INT DUE REPAID 

IDD DF PER END PER END PER 



260 
250 



2 9 0.23 
299. 66 




6 
7 
8 
9 

10 
1 1 
1 2 
13 
14 
1 5 

if 

18 
1 9 
20 

TOTAL 



645 1 
6 110 
575 9 
5396 
5021 
4634 
4234 
382 2 
3396 
2956 
2502 
2033 
1548 
1049 



. 4 1 

. 85 

. 2 2 

. 1 6 

. 31 

. 2 7 

. 65 

. 5 

. 04 

. 1 8 

. 3 

. 1 2 

. 9 7 

. 03 
95 



209 
198 
187 
175 
163 
150 
137 
124 
1 1 
9 6 . 
8 1 . 
66 . 
50 . 
34. 



. 67 
. 6 
. 17 

. 3 8 
. 19 
. 6 1 
. 6 3 
. 22 
. 37 
OS 
3 2 
08 
34 
1 

3 2 



34 

35 1 

36 3 

37 4 
387 
39 9 
412 
426 
43 9 
454 
468 
484 
49 9 
516 
532 



:il 
:ii 

. 04 

. 62 

. 6 

. 01 

. 86 

• i5 

• ?1 
. 1 5 
. 39 
. 1 3 
. 91 



3 4.64 

8000 

MONTHLY PAYMENT =S 550.23 
AND TOTALS S 11004.64 



*« DONE ** 
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Interest on Credit Purchases 



All too frequently, the rate of interest to be charged on a loan or credit pur- 
chase is hidden in the small type. After all, the bank or car dealer or finance 
company wants to convince you that you can afford that new car or home 
improvement or whatever. 

This program calculates the interest rate on a loan given the principal value 
of the loan, the number of payments, and the amount per payment. To make 
things even easier, the program will accept the cash purchase price of the arti- 
cle and the down payment and automatically compute the principal value of 
the loan. 

The sample run shows a rather unrealistic loan for an item costing $88.99. 
The down payment was $10.00 and the loan is over a period of 18 months; 
each monthly payment is $4.85. The program run shows that the actual in- 
terest rate is a modest 5.69%. 

Now try the program with a more realistic loan. A termite company in 
New Jersey advertises a total treatment for your house for only $200; just 
$29.95 down and 24 monthly payments of $11.95 per month. Is this the bar- 
gain that it seems? What is the annual rate of interest? 
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10 PRINT 

INTEREST^' 
URCHftSES " I 
20 INPUT " 
" " INPUT " 
INPUT " 
: N 

INPUT " 
■ S 

INPUT " 

70 B=R*N+D 
80 I=B-C 
90 M = M-'''.:S* 
100 T=INT< 
. 5 :w 1 
110 PRINT 
I S " ; T ; 



3 
40 

50 

60 



PRDGRfttI CfiLCULFiTES 
:"RfiTE DN CREDIT P 

PRICE DF ITEM =''!C 
DOWN PRY ME NT ="!D 
NUMBER DF PAYMENTS 

PflVriENTS PER MONTH 

flMDUNT PER PRYMEMT 



12 > 

ioo*<i*ioo^<b«m:>> + 
" "'."interest rrte 



►RUN 

PRQGRflM CflLCULRTES INTEREST 
RRTE ON CREDIT PURCHRSES 

PRICE OF ITEM =88.99 
DOWN PAYMENT =10 
NUMBER OF PAYMENTS =1S 
PAYMENTS PER MONTH =1 
AMOUNT PER PAYMENT =4,85 

INTEREST RRTE IS 5.69 'i 

** DONE ** 
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Population Growth 



So far, all the compounding problems in this chapter have involved 
money — interest, savings, and loans. But many fascinating compounding 
problems do not involve money. Consider the following problem. 

In 1960, the population figures for the United States and Mexico were 180 
million and 85 milHon respectively. The annual population growth rate for 
the United States was 1.23% and for Mexico 2.23%. If these growth rates 
remain steady, in some distant year the population of Mexico will exceed that 
of the United States. In what year will this occur? 

The program uses the method of accumulating the principal amount rather 
than a formula, and applies it to both populations. The sample run reveals 
that the population of Mexico will overtake that of the U.S. in the not too 
distant future. 

10 U=180000000 
20 M=35000000 
30 Rl=l . 0123 
40 R2= 1 . 0223 
50 y=1960 
60 U=U*R1 
70 M=H*R2 
8 Y=Y+1 

90 IF M<U THEN 60 

100 PRINT M'ERR ="5Y5" POP 

U L fl T I □ N " : " " 

110 PRINT U . S . n , , U 

120 PRINT ^'tlEXICO" , M 

> 

> R U N 

VERR = 2037 PGPULRTIQN 

U . S . . 46 1406536. 3 

MEXICO 464 46 46 9 3. 9 



But what if the Mexican people were diligently trying to bring their 
increasing population under control and their annual growth rate was 
increasing by only 0.001% per year, compared to the increase of the U.S. 
growth rate of 0.01% per year. If that were the case, would the population of 
Mexico ever exceed that of the U.S.? 

Insert Statements 60, 70, and 85 into the first program, run it and find out 
the answer to the question. 



>60 IJ«U*C 1 • 01 *R1 > 
>70 M*N*< 1 • 001*R2:) 
>S5 IF Y>9999 THEN 



100 
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Going back to the first program, run it with the population data fi*om the 
U.S. (180 million) and California (15.7 million) from 1960. At that time, the 
annual growth rates were 1.23% and 3.7% respectively. Running the pro- 
gram will indicate the year that the population of California will exceed that 
of the United States. This, of course, is nonsense. Where is the discrepancy? 

This program might make more sense if it were restated to ask in what year 
the population of CaHfornia would exceed that of the remainder of the U.S., if 
ever? 

Compounding can also be used to solve other kinds of population growth 
problems. Consider the bristleworm. The bristleworm can reproduce by split- 
ting itself into 24 segments, each of which grows a new head and tail. What is 
the maximum number of bristleworms that could be obtained in this fashion, 
starting with only one worm, after ten splittings? Assuming a splitting occurs 
every 22 days, how many offspring will one worm produce in a year? When 
will the earth be overrun with bristleworms? 

Here is another problem for which the principle of compounding is useful. 
It takes nature about 500 years to produce one inch of topsoil. Many years 
ago, the United States had an average depth of almost nine inches of this 
good dirt, but as of 1975, the country was down to about six inches. This type 
of soil is necessary, of course, for growing food. 

Careless management of our soil causes about 1% per year to erode away; 
it is then lost forever. Once soil depth reaches three inches or less, it is impos- 
sible to grow crops on a large scale. Write a program to calculate the year in 
which the U.S. will have less the 3" of topsoil, assuming that it continues to 
erode away at 1% per year. 

Will you be alive then? Will your children be alive? Will anyone be alive? 
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6 

Probability 



Statistics and probability are subjects that bring up all kinds of images. 
Some people who have not had a pleasant time in a required college statistics 
course keep as far away from the subject as possible. To other people, statis- 
tics is something with which devious manufacturers can mislead you about 
their products. 

Not long ago, one foreign automobile manufacturer boasted that 90% of 
their cars sold in the United States in the last seven years were still on the 
road. This sounds like excellent reliability. But consider the fact that this 
manufacturer was rapidly expanding in the U.S. market and 65% of the cars 
they had sold in the U.S. had been sold in just the previous two years. You 
would expect that all of these would still be running. 

If the manufacturer had sold 10% of their 7-year volume in Years 1 to 3, 
and most of these cars were not running, then the advertising claim loses 
much of its meaning. 

When approached in a logical, step-by-step manner, statistics and probabil- 
ity are not difficult. In fact, they can be a great deal of fun. 

Some of the programs use a formula while others simulate the event in 
which we are interested. The results are the same, but you may find that the 
simulations help you to understand exactly what is happening. 
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Pascal's Triangle— Calculated 



Pascal's Triangle is quite fascinating. As you can see from the portion of 
one reproduced below, each row is symmetrical. Each row also contains the 
coefficients for a binomial expansion. The sums across the ascending diag- 
onals form the Fibonacci sequence. The sums across the rows are all powers 
of two. Each row corresponds to the digits of a power of 1 1. Every element is 
the sum of the two above it. And, in case you care, all the elements in it are 
identities in combinatorial theory. 

The ways this marvelous triangle can be generated are as varied and in- 
teresting as its properties, though perhaps more difficult to figure out. Here is 
one way to use a computer to generate the triangle. 

Any element can be found be adding together the two elements immedi- 
ately above it. The program uses this principle to produce a triangle eight lev- 
els deep. Lines 10-30 set the rightmost diagonal to 1. Each element is stored 
in the two-dimensional variable P(R,C) with R denoting the row and C 
denoting the "column." The variable T (Line 50) simply leaves some blank 
spaces so the output resembles a triangle. The crux of the calculation is in 
Line 70 in which the value of each new element is calculated. 

There is another interesting way to calculate Pascal's Triangle in even 
fewer statements than the program here. It generates the triangle one element 
at a time and does not use any arrays or two-dimensional variables. Can you 
figure out how to write such a program? 
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10 FDR C=l TO 9 

2 P>: 0,0 = 1 

3 NEXT C 

40 FOR R=1„T0 9 
50 T=17-R*2 

6 FDR C=2 TO R+1 

7 P < R + 1 J C > = P < R , C 5 + P t: R , C ■ 
80 IF P<R,CO=0 THEN 110 

9 PRINT TfiB t; T > ; P < R , C > ? 
100 T=T+4, 
110 NEXT C 

12 PRINT 

13 PRINT 
140 NEXT R 



►RUN 

1 

1 1 
1 2 1 
13 3 1 
1 4 6 4 1 
1 5 10 10 5 1 
1 6 15 20 15 6 

^7 21 35 3 5 21 7 

* * n D N E * * 
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Pascal's Triangle— by Probability 



This program simulates the dropping of balls through a triangular array 
shown below. At each level, a ball is equally likely to fall either to the left or 
right. At the bottom of the array a cup is placed at each end point; these cups 
collect the balls. After each group of balls has been dropped, the number of 
balls in each cup is tallied and displayed. 




u u u u u u 



B Y " 
20 
3 

II ■ If 

LS " 

5 
60 

7 

8 

9 
100 
i 10 
120 
13 
140 
150 



Pf^SCRL'S TRIRNGLE 
i . P R 9 B R E I L I T Y " ! " " 
INPUT "BALLS T □ DROP = " : M 
INPUT "NUMBER QF LEVELS = 



PRINT " 

■ II II 

N = i 



" : "POSITION", " B R L 
TO M 



FDR 
T = 

FOR L=i TD K 

I F^.RND > .5 THEN 90 

NEXT L 

b<:t+i> = b<t + i:) + i 

N E X T N 

FDR L=l TD K+1 
PF?INT,L, B<L) 

B <: L ;> = 

NEXT L 



PgSCgL'S TRIRNGLE BY 
PRDBRBILITY 

BRLLS^TO DROP =1000 
NUMBER OF LEVELS =2 



POSITION 

1 

2 
3 



B R L L S 

234 
5 IS 
248 
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■^PRSCfiL'S TRIANGLE BY 
PROBABILITY 

BRLLS TD DROP^=i00g 
NUMBER DF LEVELS =4 



PDS I T I QN 

1 
2 
3 
4 
5 



BALLS 

51 

229 
40 7 
240 
7 3 



DONE 



The number of balls landing in the various cups at each level, when re- 
duced to the lowest common denominator, should approximate the numbers 
in Pascal's Triangle. Do you know why they should? 

Try a few runs of the program with a different number of balls. Do the 
numbers obtained really approximate those in Pascal's Triangle? If, at each 
dividing point in the array, every other ball went in the opposite direction, 
then Pascal's Triangle would be produced. Try running the program with 4 
balls at Level 2; with 8 balls at Level 3; and with 16 balls at Level 4. Do your 
results look like those of the previous program which calculated the triangle? 
Does the approximation come closer as the number of balls is increased? 

How can you determine how close your results from this program are to 
the exact value of Pascal's Triangle? One way is to take the number of balls 
that dropped into each cup on a given level and divide that by the total num- 
ber of balls divided by the theoretical sum of the row. So, in the sample run at 
Level 4, you would divide the balls in each cup by 62.5 (1000/16). You then 
compare that to the "correct" number and compute the percent difference. 

Here is the result of this procedure for Level 4 (actually the fifth row) in 
the sample run: 



Cup 


Balls 


■i-62.5 


Correct 


Deviation 


1 


51 


0.82 


1 


18.00% 


2 


229 


3.66 


4 


8.50% 


3 


407 


6.51 


6 


8.50% 


4 


240 


3.84 


4 


4.00% 


5 


73 


1.17 


1 


16.80% 



Why are the outside cups "ofT' by a greater percentage than those closer to 
the center? Does this also happen on other levels with different numbers of 
balls? 
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Common Birthdays 



Here is an interesting little problem in probability. In a group often people 
selected at random, what is the probability that any of them will share the 
same birthday? How about a group of 20 people? Of 50 people? 

Conversely, how many people would you need in a group such that there is 
a 50% probabiHty that at least two of them have the same birthday? How 
many people would be needed for a 90% probability of an overlap? 

Try to answer these questions, either by guess or by calculation, before you 
look at the output from the program. 

This program provides a painless introduction to the world of statistics. 
The calculation is actually quite trivial. The probability that any person in a 
group has a birthday on January 1 is 1/365. If our group has only two people 
in it, the probability that both of them have a birthday on January 1 is 1/365 
times 1/365, or a very small number indeed. The probability that they have a 
common birthday is 365 times the very small number just obtained, or about 
0.27%. 

However, if there are three people in the group, the probability goes up 
slightly. Call the three people Betsy, Ken, and Larry. From the reasoning 
above, we know that there is a 0.27% probability that Betsy and Ken have 
the same birthday; also a 0.27 % probability that Betsy and Larry share a 
birthday; and finally, a 0.27% chance that Larry and Ken have a common 
birthday. Hence, the total probability for a group of three people is three 
times the probability for just two people. 

A group of four increases the probability over that of two people by a fac- 
tor of six, five people by a factor of 10, six people by 15, seven people by 21, 
and so on. There is a progression here, but the probability can also be cal- 
culated by the formula: 



P = 1 . 365-N 
365 

Can you see why this formula produces the same result as the description 
in words and associated progression above? 

Consider all the presidents of the United States. (How many have there 
been to date?) Two of them, James Polk and Warren Harding, were bom on 
November 2. Is this to be expected given the size of the group? 

Since birthdays can be predicted, at least statistically, it ought to be pos- 
sible to predict deaths as well. It is interesting too that John Adams, James 
Monroe, and Thomas Jefferson all died on July 4th. Millard Fillmore and 
William Taft both died on March 8. What is the probability of that set of 
events? 
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10 PRINT T R B 1 1 > 
Y THAT'* 

2 PRINT ''ND. DP 
T TWO HRVE" ^ 

3 PRINT "PEOPLE 
E BIRTHDRY" 

4 Q = 3 6 4 3 6 5 

5 FDR N=2 TD 4 
60 P=100^^<: 1-Q> 



'PRDBRBILIT 
RT LEflS 
THE SRM 



b U — I U U V. i ~ W . » 

7 PRINT N ; T R B 1 3 ) 5 I N T < P * 1 
+ . 5>^' iOOl 'y^i ^ ^ 
" ~ ' ^ " N > y 3 6 5 



SO 
3 



Q = Q* < 36 5 
NEXT N 



> R U N 

ND . UF 
PEOPLE 



4 
5 
6 
7 
S 
9 

1 
1 1 
12 
13 



PRDBRBILITY THRT 
RT LEAST TWO HAVE 



THE 



SAME BIRTHDAY 
2 7 X 

82 

.64 

.71 y. 
. 5 >: 
. ^ z 
.43 >: 
. 4 6 y. 
1.69 y. 

4.11 ^. 
6.7 

9.44 ;i 



30 
3 1 
32 
33 
34 
3 5 
36 
37 
3 8 
3 9 
40 



7 , 
73 , 
75 , 
77 , 
79 , 
81 , 
83, 
84 , 
86 
87 
89 



63 
05 
33 

i3' 

44 
22 
87 
41 
82 
12 



Now it is your turn to write a program. Here is a game to be played among 
your friends. You make a bet that they have similar preferences in colors. 
Each person in the group puts up one penny and you, because you are so sure 
of yourself (and probabihty), put up an amount in cents equal to the number 
of people in the group. If you win, you get to keep all of the money; if your 
friends win, each of them gets double his original bet, or two cents. 

Then each person selects a color (from a larger list than there are people, of 
course). If any two have picked the same color, you win; if all have picked 
different colors, they win. 

Since you would like to win, you need to know how many colors should be 
on the Ust for different size groups to give you a better than 50-50 chance of 
winning. You might Hke the probability to be up around 70% or so. You can 
produce the necessary table with a five-line program. Go to it! 
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Coins in a Poclcet 



The next two programs were originally written by Glenda Lappan and 
M.J. Winter and appeared in Creative Computing magazine. They are marvel- 
ous simulations for illustrating various aspects of probability. 

Coins in a Pocket is a simulation of the following situation. A newspaper 
costs 5 cents. A customer has 5 pennies and a dime in his pocket and offers to 
pay for his paper by letting you, the vendor, select at random two of the six 
coins. If you and this customer repeat this procedure for the 20 working days 
of a month, how much more or less than $1.00 (20 days x 5 cents) are you 
likely to have collected? 

The program below solves this little problem. The first random coin is se- 
lected from a group of six (Line 70 in which X = INT(RND*6)). The value 
of X can be 0, 1, 2, 3, 4, or 5. If it is 0, we assume the dime was selected and a 
second pick is not made, because it will surely be a penny. Thus 1 1 cents is 
added to the running total in Line 130. 

If the first coin is a penny (X = 1, 2, 3, 4, or 5) then we make a second pick 
from the five remaining coins. Again, if the dime is chosen (Y = 0), 1 1 cents 
is added to the total; otherwise 2 cents is added. 

As you can see from the sample runs, after a great number of trials, the 
average amount of money collected each day seems to be close to 5 cents. If 
you would like to convince yourself that the answer is, in fact, 5 cents, con- 
sider the following. Assign a letter to each coin, A - F. Let A be the dime and 
B through F be the pennies. Since all combinations are equally probable, here 
are all the possible combinations: 

AB 

AC BC 

AD BD CD 

AE BE CE DE 

AF BF CF DF EF 
There are five combinations with a dime (5x11 cents) and ten combina- 
tions of only pennies (10 x 2 cents). Add it up and you have 55 cents plus 20 
cents divided by a total of 15 combinations which equals an average value of 
75/15 or 5 cents. 
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II 

40 INPUT "NUMBER OF TRIRLS = 
" : N 

50 U = 

60 FDR K=1,T0.N 

70 X=rNT<RND*6.:' ^„ 

SO IF X=0 THEN 130 

•30 Y=rNT<RND*65 

100 IF Y=0 THEN 130 

110 U=U+2 

1-20 GOTO 140 

130 U=U+11 

140 NEXT K 

iio PRINT "RVERflGE VRLUE OF 
COINS IS" ! U/N5 "CENTS" 



PRQGRflM SiriULRTES TfiKING TWO 
COINS AT RflNDOM FROM R 
DIME AND 5 PENNIES 

NUMBER OF TRIRLs =1000 

AVERAGE VRLUE OF COINS IS 
4.934 CENTS 

** DONE ** 

^PROGRAM SIMULATES^TAKING TWO 
COINS AT RANDOM FROM H 
DIME AND 5 PENNIES 

NUMBER OF TRIALS =100 

AVERAGE VALUE OF COINS IS 
4.8 8 CENTS 

** DONE »* 
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Baseball Cards 



In statistics and probability, it is frequently necessary to predict the num- 
ber of trials on average until one is successful. For example, if you are trying 
to roll a four with one die, on average how many tries will it take until you do 
so? 

Since a die has six sides, the probability of rolling any number between 1 
and 6 is p = Ve. Thus, on any given roll, you have a Ve chance of rolling a 
four and chance of rolling something else. This failure is designated q. So 
to be successful in rolling a four, we have a Ve chance on the first roll. On two 
rolls, the probability is two trials x the probability of one failure, then success 
(2 X Ye xVe). On the third roll, the probability of success is three trials x the 
probability of two failures followed by one success (3 x x Ve x Ve). 

Continuing this reasoning leads to the formula for the expected number of 
trials to success: 

E=lpH-2qxp + 3qxp-h4qxp + ... 

This series can be reduced and solved for E which leads to: 

E= -L = -L 

l-q p 

So now the answer to the original problem can be calculated. The expected 
number of trials to roll a four is 1/p = 6 tries until success. 

But there is another way to approach this type of problem with the assis- 
tance of the computer. Consider the problem. Assume there are ten different 
prizes in Crunchies cereal boxes. How many boxes of cereal would you have 
to buy to obtain the complete set? The formula above can be expanded to 
solve this problem: 

N ^ N-1 ^1 ~^^^+T^"3"^ ^N" 

For a value of ten, you can solve this formula by hand. However, let's say 
you would like to solve this problem for baseball cards as found in packs of 
bubble gum. If there are 50 cards in a complete set, how many packs of gum 
must you buy, on average, to get a complete set. Write a computer program 
using the general formula above to solve for a set of any size. The answer for 
50 cards is 224.96 packs of gum; how many would you have to buy for a set of 
100 cards? 

There is another way of approaching this problem. This method is similar 
to that used to randomly select coins out of the pocket. In this case, the ran- 
dom group is set equal to the total number of cards; this value is accepted as 
input in Statement 20. Each purchase is set equal to a random number be- 
tween 1 and N (Statement 150). This card is then put into its proper place in 
the collection (Statement 160). However, if there is already a card there from 
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a previous purchase, we simply increment the purchase counter (Statement 
190) but do not get any closer to obtaining a complete set. After each pur- 
chase, we test to see if the set is complete (Statement 200), otherwise we go on 
buying more packs of bubble gum. 

When the set is complete, the number of packs of gum are tallied up and 
printed out. After a set of trials, the average is computed. This averaging 
value should be reasonably close to the value obtained by the formula al- 
though it will take a great number of trials before the two numbers are within 
1% of each other. 

IG PRINT ^'PROGRRH SIMULRTES 
Buy ING" : "BUBBLE GUM WITH BRS 

2g^inplit^^ "crrds . crrds in se 
30 'if""n<ioi then so 

40 PRINT "SORRY, 100 IS Mfl^CI 
MUM" 

5 sum 2 

6G_INPUT "NUMBER DP SIMULRTI 

N S = : K 

7 DIM C < i > 

SO PRINT " 2 TR I BL " , " PRCKS " 
9 S = 

100 FOR 1^1 TO K 
110 n = 

120 FOR J=l TO N 

13 C <: J ) = 

14 HEKT J 

15 X = I N T < R N D * N ) + 1 

16 C X > = C C X ) -f 1 

17 IF C < X = 1 THEN 190 
180 GDTD 150 

ISO I) = D+1 

200 IF^D=N^THEN 220 
210 GOTO 150 

220 T=^0 

230 FDR J=l TD N 
240 T=T+C<J> 
250 NEXT J 
260 PRINT I,T 
2 7 S = S + T 
280 NEXT I 

2 9 PRINT " fl V E R R G E " ^ S .-^ K 

> R I J N 

PRCIGRRM SIMULRTES BUVIKG 
BUBBLE GUM WITH BRSEBRLL 
CRRDS. CRRDS IN SERIES =50 
NUMBER OF S I M U L R T 1 □ N S = 5 



TRIRL PRCKS 

1 2 8 4 

2 246 

1 5 3 

4 22 2 

5 29 
RVERRGE 239 

* * D □ ^^ E * * 
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1^ 

PROGRfitI SIMULATES BUYING 
BUBBLE GUM WITH ERSEBflLL 
fi!iSP§A CRRDS IN SERIES =10 
NUMBER OF SIMULflTlDNS =10 



TR I RL 

2 
3 
4 

t) 
7 
S 
9 

1 

flVERFIGE 



P R C K S 

1 4 

28 
27 

1 9 
28 
25 

2 
26 
22 



DONE 



If you run this second program for a set of 100 cards, it will sometimes run 
for a very long time before arriving at the answer. Try it with a set of 500 
cards as are used by some real bubble gum manufacturers. You can go have 
your dinner while the program computes just the first value. Be sure to 
change the dimension statement in Line 70 to C(500). 
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System Reliability 



As more and more people in the world come to depend upon mechanical 
and electronic devices in a myriad of different ways, it is important that these 
technological devices continue to function. 

Some years ago, the military came up with a measure that could be applied 
to all kinds of systems, big and small, to measure reliability. It is called 
"mean time between failures." What this means is the length of time, on av- 
erage, between breakdowns. For a tank, this may be 100 hours, while for a 
spacecraft, the MTBF must be considerably longer than the planned mission. 

As we saw in an earUer section, it is frequently desirable to break down a 
large problem into smaller subproblems. To calculate the MTBF for a space- 
craft would be quite impossible. Instead, it is necessary to start with smaller 
systems and build up to the whole. 

Consider one of the electrical subsystems of a spacecraft. It uses five 
components as shown below, two parallel systems A and B, arranged in se- 
ries. The parallel subsystems are said to be "redundant." This is one method 
of increasing reliability since the system Will continue to work if at least one 
of the components works. 



D 



a 



A B 

If the manufacturer of the components stated that each one has a 60% 
probability of lasting 1000 hours, what is the chance that the entire system 
will last 1000 hours? Take a guess before reading on to the solution below. Is 
your guess greater than 60% or less? Why? 

The program below is a simulation of this system. Remember in subsystem 
A, it will continue to work if any of the three components works and in B if 
either of the two components works. However, the system will fail if all three 
of the A components or both of the B components fail. 

In Statement 40, the program is set up to make 500 trials of the system. 
Statement 60 selects a random integer between 1 and 10 for each component 
for each trial. If the value of this integer is 6 or under, the component is okay 
(60% probability of working at the end of 100 hours). If it is 7, 8, 9, or 10, 
this means the component has failed, and the program prints "Bad!" 
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1 

2 
USE 
URS 
30 

4 

5 
€.0 
70 
SO 
90 
100 
1 10 
120 
130 
140 
1 50 
160 
170 
ISO 

1 9 
200 

2 1 
220 
230 
" ; I 
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RANDOM I ZE 

PRINT "PROGRAM S I M U L fl T E £ 
^OF" : "CIRCUIT FOR 1000 HO 

S = 

FOR 1=1 TO 500 
FOR .J = l TO 5 
C C J > = I N T C R N D * 1 > + 1 
IF C< J)> >6 THEN 100 
(J C J > = 1 
GOTO 110 
C <: J ) = 
NEWT .J 

IF C':;i) = I THEN 170 
It Qt:2)=l THEN 170 
i E , P. i -3 :> = 1 T H E N 17 
PRINT " B fi II ! " J 
GOTO 220 

JE '<<.i'> = l THEN 210 
IF C<5)=1 THEN 210 
PRINT "BfiD ! " 5 
t5DTD 220 
S = S + 1 
NEXT I 

^ PRINT. "RELIftBILITY ABOUT 
NT<. 100U0*CS''500> + . 5>''100; 

END 



P R □ G R R M 
C 1 R C U I T 



SIMULATES USE 
FOR 1000 HO U R : 



□ F 



BRD ! 

BRD ! 

BRD ! 

BRD ! 

BRD ! 

BRD ! 

BRD ! 

BRD ! 

BAD ! 

BRD ! 

BRD ! 

BRD ! 

BRD ! 

BRD ! 

BRD >. 

BRD ( 

REL I 



BRD 

BRD 

BAD 

BRD 

BRD 

BRD 

BRD 

BRD 

BRD 

BRD 

BRD 

BAD 

BAD 

BAD 

BRD 

BAD . 

AB I L 



BAD 
BAD 
BRD 
BRD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
I T V 



* * DONE * * 



BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BRD 
BAD 
ABQ1, 



BAD ! 
E AD ! 
B AD ! 
BAD ! 
B AD ! 
BAD! 
BAD ! 
BRD! 
BAD! 
BAD ! 
BAD ! 
BRD ! 
BAD ! 
BAD ! 
BAD! 
BAD ! 
T 77 



BFlD 
BAD 

BAD 
BAD 
BAD 

BAD 

BAD 
BAD 
BAD 
BRD 
BAD 
BAD 
BAD 
. S 



BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BAD 
BRD 
BRD 
BRD 
BAD 
BRD 



Statements 120-140 tests whether subsystem A works and Statements 170- 
180 checks out subsystem B. If both work, a success is counted in Statement 
210. 

Now it is your turn to use what you have learned in this chapter and write 
a program to determine the mean time between failures of the system above. 
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7 

Geometry and Calculus 



In this chapter several of the problem solving approaches from earlier 
chapters are used to solve geometric problems. You will find that many prob- 
lems can be solved in a variety of different ways — by applying a formula, by 
trial and error, by successive approximations, and, in some cases, by common 
sense. Perhaps the most important thing to learn from these problems and 
programs is how to analyze a problem to reach the solution quickly and 
painlessly. 
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Crossed and Slipping Ladders 



Here is a simple problem of a slipping ladder. A ladder 25 feet long is 
placed so its foot is 7 feet from the base of a building. The base of the ladder 
slipped on some loose gravel so that the top is 4 feet lower than where it was 
to start. How far did the foot of the ladder slip? 




The diagrams show the two positions of the ladder. By the Pythagorean 
theorem we know that a^ + b^ = c^, hence, the equations needed to solve the 
problem are: 

X = VIV' - V 



y = V25^ - (x-4)2 

and the amount of slippage of the base is z = y - 7. These three equations are 
put into a computer program which quickly calculates an answer of 8 feet. 

While the aritmetic in the above problem is not particularly messy, it is still 
no great joy to solve by hand. However, the computer is just as happy to do 
the problem with really messy dimensions, say a ladder length of 27.83 feet 
and a distance from the wall of 7.62 feet. 

Id PRINT "SLIPPING LRDDEf^:" 

20 C=25 

3Q B = ? 

4 X = S Q R < C ^ 2 - B ^'^ 2 ) 

5 Y = S Q R < C 2 - < X - 4 ) ^ 2 :> 

6 2 = y - B 

70^PRINT "BRSE S L I PPED " 5 Z 5 " F 
>RUN 

SLIPPING LfiDDER 
BRSE SLIPPED 8 FEET 

•* * D f J E * *■ 
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>20 C=2 7.83 
>30 B=?.62 

> R U N 

SLIPPING LADDER 

EASE SLIPPED 8.386131551 

FEET 

•** DONE 



Let's consider an old problem found in many classic collections. Two lad- 
ders, one 25 feet long and the other 35 feet long lean against buildings on op- 
posite sides of an alley as shown below. The point at which the ladders cross 
is 12 feet above the ground. How wide is the alley? 




X y 

By using similar triangles twice we find that 

12/a + 12/b = 1 
Then, by applying the Pythagorean theorem and reducing, we obtain: 

a^ - b^ = 600 

Using one of the methods described in the Problem Solving chapter, you 
can solve these two simultaneous equations. Or they can be combined into 
one equation in which z, the width of the alley is 

z = V352 . ^2 

By eliminating b, the following fouth degree equation is obtained 
. 24a^ - eOOa^ + 14400a - 86400 = 

Again, this can be solved using one of the methods in the chapter on Prob- 
lem Solving. This is the traditional approach, but there is another. The solu- 
tion is the intersection of the curves described by the two original equations. 
By using successive approximations, say in steps of 1 for a, you could solve 
for bj and bg in the following restated original equations 
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b2 = 1/600 - 

When b, falls below \, reduce the step to 0. 1 to obtain a closer approxima- 
tion. By continuing this method of successive approximations, it is possible to 
obtain a very accurate solution. 

Is there another approach? Yes, there is and it also avoids the quartic equa- 
tion. It uses the original equations in a trial and error procedure as described 
in the chapter on Sets and Repetitive Trials. See if you can write a program 
using this approach. 

Here is another classic problem for you to solve in any way that you like. 
What is the longest ladder that can be carried in a horizontal position around 
the corner made where a 12-foot wide alley meets one that is 8 feet wide. The 
diagram shows the problem. 
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Distance Between Coordinate Points 



This program solves for the distance between any two points in three- 
dimensional space defined by their x, y, and z coordinates. 

The formula for the distance between two points in three-dimensional 
space is: 

d = (X, - x,)^ + (y, - y,)^ + (z, - z,r 

It would be desirable for the program to be able to solve any problem of 
this kind. One approach would be to use input statements to accept the point 
coordinates. Another is to use a data statement to define the points. Then 
only this one statement has to be changed for a new problem or set of 
problems. 

The program is written to calculate the distance between three sets of 
points. The points used were: 

0, 0, and 3, 4, 5 

3.5,-4.7,6.2 and -0.9,3.0,4.4 
67, 36, 82 and 54, 25, 90 

This calculation is used extensively in aerospace navigation. Can you deter- 
mine the angle or "compass heading" of the resulting flight path? Better start 
with just two dimensions. 
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^ DibTfiNCE": '^FETUIEEN POINTS I 
N JHREE-" : ■^DiriENSIDNRL SPACE 

2 PRINT POINT'* 

3Q PI? I NT "CDDRDINRTES D I 

5S F-^SSr.^ ' B , C , D , E. F 

5 J. = SQR <. < R-rO -^2 + < E-E ) "^2+ < C- 

60 PRINT R 5 B ? C 

'5 G □' T □ 4 

log SQTfl 0,0,0,3,4,5 

110 DfiTfi 3. 5, -4. 7, 6. 2,-. 9,3, 

.120 DRTfl 67, 36^ S2, 54, 25, 90 

> 

> R U N 

DIMENSIDHflL SPflCE 
POINT 

COaRDIHflTES DISTANCE 


3 4 5 7.071067812 

3.5 -4.7 G.2 
-.9 3 4.4 9.049309366 

6 7 3 6 82 

54 25 90 18.81488772 
* DRTfl ERROR; IN 40 
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Area — by Calculation 



It is a simple matter to calculate the area of regular geometric figures by 
using the usual formulae. However, the problem becomes more difficult when 
it is necessary to calculate the area of two combined regular shapes or the 
area of irregular shapes. This program shows the method of analysis for solv- 
ing a problem of the first type, while the next program demonstrates four 
methods of dealing with irregular areas. 

The problem is to solve for the shaded area of the figure below for any 
value of the radius, R. 




The first step is to recall the formulae for calculating the area of a circle 
and square: 

A (circle) = pi * 
A (square) = Side^ or (2 * R)^ 
The difference in area between a square and a circle inscribed within its 
borders is: 

A (difference) = A (square) - A (circle) 
and the area of one comer is the difference divided by 4. The program below 
will calculate the area for any radius. 

Now it is your turn. Write a program to calculate the difference in area 
between a circle and square in which the square is inscribed within the circle. 
Now, change your program to calculate the difference in area for a triangle, a 
hexagon and an octagon inscribed within a circle. 



10 PRINT "PROe 
CRIIeD WITHI 



RM 
J fl 
•\ fl 




INPUT 'RQDIUS 
§s c; 2#R > -^'2 



. RINT 
PRINT 



CftLCU 



CTRCL1 
SQUflr 

= " 8 R 



TRAPPED AREA 
DNE CORNER = 



.flTES 
INS 



5 Dl 
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>RUN 

PRgeRfiM CPlLCULRTES flREfl 

INSCRIBED 

RfiDIUS =10 

TRRPPEB^ flREfl = 85.841 
□NE CORNER = 21.46025 

** DONE »♦ 

>RUN 

PRgSggM CRLCULftTES RREFI 

S!?Sf^^^;giggk^ Inscribed 

?8§iy§^=3A 14159 
IR9P£gE.,9P'Eft = 472152802 
ONE CORNER = 2.1180382 

DONE ** 



Extend your program to calculate the difference in area for any regular fig- 
ure inscribed within a circle of radius 1.0. Set up a table or results as follows: 
Number of sides Difference in Area 

3 
4 
5 
6 
7 
8 

What can you conclude from these results? Do these areas follow in some 
sort of progression? 
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Area — by Integration 



In many cases it is necessary to determine the area of an irregular figure or 
the area under a curve where an exact formula is not available. The approach 
most commonly used is to divide up the enclosed area into small regularly 
shaped pieces and sum up the areas of all of these pieces. 

The easiest shape to use in these calculations is a rectangle. A group of 
rectangles can either be inscribed within the irregular figure or curve, or 
circumscribed around it. The first two diagrams show these two methods be- 
ing used to find the area under a curve. A third method is to use trapezoids 
which, depending upon the direction of curvature, will either be inscribed or 
circumscribed automatically. 




Inscribed Circumscribed Trapezoid 

Rectangle Rectangle 
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A fourth method, known as Simpson's Rule, essentially fits a series of 
parabolas between the points of the curve and calculates the average area. It 
requires that the area be divided into an even number of parallel slits. Let us 
call the total number of divisions 2m which are h distance apart. The point on 
the curve where the first line intersects is yo, the second yi, and so on until 
y2nT The area is then given by the formula: 

A= VMiyo + y2m) + 4(y, + yj + ... 

+ y2n,-l) + 2(y2+ y4 + • . • y2m-2)] 

The area as calculated by Simpson's Rule converges extremely quickly 
compared to the other methods. Nevertheless, as the number of intervals in- 
creases, all the methods approach the same limit. As might be expected, the 
trapezoidal approach converges more quickly than either of the approaches 
using rectangles. Compute the average of the two methods using rectangles. 
What do you get? Does this suggest another method? 

The methods used in this program involve the calculus. And you thought 
the calculus was difficult! Now you know otherwise. 



1 
REfl 
MET 
RflT 

2 
HE 

3 

n y 

4 

ST 

5Q 

6 
SCR 
RCU 

7 
SO 
90 

1 
1 1 
120 
130 
140 
150 
160 
170 
1 3 
190 
200 
210 
220 



PR I NT 
U N D E 
HDDS 
I D N " t 
PRINT 
100" 
PRINT 
FILUES 
INPUT 



" P R Q G R fi M C D M P U T E ^3 Fl 
R " 5 " fi CURVE BY FOUR 
of;; ! "NUtlERICRL INTES 

"DEFINE CURV/E IN LI 

/'INPUT STRRT AND EN 

"FOR X <SMFILLER FIR 



PRINT 
PRINT 
R E C T 
M REC 
PRINT 
M«-2 
S = 
DEF 
11 = M + 
FDR 
C = S 

S = 
P = 

D = 2'*' 

FDR 
X = FI + 
P = P + 
Q = Q + 
NEXT 



B 



„"Na. INTERVRLS" s " IN 
S. TRRPEZQIDS" s "CI 
TS PflRFIBDLRS" 



FNfi<X) =X'-^3 
NsM TO M+2 



N 

-fl > 

1=0 TO D-1 

H * F N fK X J 
H * F N fl < X + H > 
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T = C P + Gl > ' 2 

FOR J=2 TD <D-2>STEP 2 
L) = U + 2*FHfl<ft + J*H) 
NEXT J 

FOR K=l TD <B-1>STEP 2 
V = V+4*FNn<ft + K*H:) 

RI,NT DS P, T: ©.» V:; 

r'i5<§4 THEN 110 
F flBSC (:c-s>/< <C. + S>/2) ) 
THEN 110 



>RUN 



PROGRfitl CaMPUTES^I 
ft CURVE Sy FniJ5„M 
NUMERICAL INTEGRA 



DEFINE CURVE 
j||uT STRR- 



STfiRT BN 
CSMSLLER 

MTERVfiLB 



lSsC^"RECtS' 
fc^Rcfin RECTS 



IREfl^UNDER 
THDDS OF 
■ ION 

_IME 100 
ND VALUES 
RST:) = 1 , 10 



TRfiPEZDIDS 
PARRBGLAS 



?53 . 1 875 



524S 
4 

150 1 - 
374S , 
8 

1969, 
3 9 3, 
1& 

2226 ■ 

2788 

32 

2 3 6 1 
2642 
64 

2429 
2570 



6875 

17 1375 
92137 5 

136719 
011719 

6 12 3 5 
54 930 5 

223389 
192139 

, 997253 
, 481 628 



3000 . 
2499 . 

2625. 
2499, 

■ 

2 4 9 9, 

2507, 
2499, 

2 5 1 
2499 

2 5 
2 4 99 



9375 
75 

046875 
7 5 

7 4 219 
75 

58 1055 



■76 4 



70 
75 



23944 1 
75 



** DONE ** 
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8 

Science 



Using techniques and approaches presented in the previous chapters, this 
chapter contains five programs in the area of science. One uses a formula to 
solve simple gas problems, two are drill exercises on the gas laws of Boyle and 
Charles, and the last two are simulations. You will see that the simulations 
draw upon many previous techniques such as progressions and repetitive 
calculations. 



Ill 



Gas Volume 



Here is a simple program to produce a not-so-simple table of values for gas 
volumes. 

The volume of a gas varies directly with the absolute temperature T (Kel- 
vin) and inversely with the pressure P. If a certain quantity of gas occupies 
500 cubic feet at a pressure of 53 pounds per square foot and an absolute tem- 
perature of 500 degrees, what volume will it occupy at 600 degrees absolute 
temperature and pressures from 100 to 1000 pounds per square foot in in- 
crements of 50 pounds? 

The original conditions are used to solve for the constant K in Statement 
40 (K = V*P/T). Then new volumes are computed for varying pressures 
with T equal to 600 degrees. The formula used is V = K*T/P. 

In the second part of the program, Lines 70 to 100 are replaced to produce 
a plot of the gas volume for the various pressures. 

How would you modify the program to deal with a more general case (i.e., 
other gasses and different temperatures)? Second, can you write a program 
that produces a table of values for a gas at different pressures and 
temperatures? 





RS VOLUMES 



TEMPERATURE c:K ) =<$00 



PRESSURE 
1 



2 D 
2 50 
30 
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4 50 

5 00 
5 50 
600 
650 
700 
7 50 
800 
850 
•3 
9 50 

1 00 

■>r * DONE * * 



70 . 66666667 
63.6 

57. 818 ISl 82 

■4I. 92307692 
4 5. 42857143 
4 2.4 
39.75 

3 7.41176471 
35. 33333333 
33 . 47368421 
3 1 . S 



10 V=500 

20 P=53 

30 T=300 

40 K = S'*P.'T 

50 PRINT "PRDGRftM COMPUTES G 

as VOLUMES" : "AT DIFFERENT PR 
eaSURES " s " " 

60 INPUT "TEMPERATURE CK> =" 

70 PRINT "PRES.'IOO PLOT 
OF VOLUME" 

80 FOR P=100 TO 1000 STEP 50 
90 V = K*T1.'-P 

100 PRINT P/IOOJ TftB'':4 + INT<V-' 
I 3 ) > S " * " 
110 NEXT P 



■^pHBsRflM COMPUTES GfiS VOLUMES 
fif DIFFERENT PRESSURES 

TEMPERFITURE K ) =600 

PRES^lOO PLOT OF VOLUME 

1 * 
1.5 * 

2 * 
2.5 * 

3 * 
3.5 * 

4 * 
4.5 * 

5 * 
5.5 * 

6 * 
6.5 * 

7 * 
7.5 * 
S * 
8.5 * 

9 * 
9.5 * 

10 * 
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Charles' Law Drill 



In the previous program, gas volumes were calculated using Charles' Law 
and Boyle's Law. Here is a drill and practice program (remember Chapter 1?) 
that produces problems relating the volume and temperature of a gas. When 
pressure is constant, the volume/temperature relationships can be stated as 
follows: 

Vq = Vj_ 
To T, 

The program presents four problems, one to solve for each of the four vari- 
ables in the equation above. 

How can you make the program more efficient? More interesting? 



10 RflMBOniZE 

20^^PRINT "DRILL ON CHARLES' 

30 'print "vglumes in MILLILI 

;; TEMPERATURES IN DESRE 

4 V = i N T < R N D * 5 + 5 > * 1 

50 V1=INT<RND*50+50>*100 

60 T=INT<RND*125+125> 

70 T1=INT<RND*125+125> 

y c = c + 1 

90 ON C GOTO 100,140,180,220 

? ^ 6 U 

100 v=o 

118 g?iPlT;!^¥^^ 

170 GOTO 270 
ISO T=0 

ilO^PRINT "CALCULATE THE VAL 

200 Ql=V*Tix'Vl 

210 GDTD 270 

220 T1=0 

230 PRINT "SOLVE FOR Tl SIVE 

240 Q1=T*V1/'V 

;^50 eOTD 270 

i'bO STOP 

270 PRINT " V VI T 

3 10 PRINT 
320 GOTC 40 

114 



•RUN 

DRILL □ N 



CHARLES 



LAW 



VOLUMES IN MILLILITRES 
TErfPERflTURES IN DEGREES K 



SOLVE FOR V WHEN 

V VI T 

9 8 16 9 

YOUR FINSI..JER =6600 
CORRECT VFILUE = 6900 

WHAT IS VI GIVEN 

V V 1 T 
6 8 15 2 

YOUR ANSWER =6500 
CORRECT VALUE = 653 1 

CALCULATE THE VALUE 

V V 1 T 
6800 S900 

YOUR ANSWER =161 
CORRECT VALUE = 158 



T 1 

240 

S33333 



T 1 
146 

573947 

OF T 
T 1 
208 

9213483 



SOLVE FDR Tl GIVEN 
V VI T 

S600 6300 248 
YOUR ANSbJER =200 
CORRECT VALUE = 196 



Tl 



0930233 



** DONE ** 
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Boyle's Law Drill 

Boyle's Law describes the behavior of gases under ideal conditions when 
pressure and volume are varied. When the temperature is constant, Boyle 
found that: 

Po * Vo = P, * V, 
Pressure is normally measured in centimeters of mercury while volume is 
measured in millihtres. As with the drill on Charles' Law, this program 
presents four problems, one to solve for each of the four variables in the equa- 
tion above. 

Unlike many other drill and practice programs, these two do not compare 
your answer with the correct one. Instead, they leave that up to you to do. Is 
this desirable? Why or why not? If you feel it is undesirable, change the pro- 
grams so they do compare answers and calculate a score. 



10 RfiNDDMIZE 

20 PRINT "DRILL OH BOYLE'S L 
flW" : " " 

30 PRINT "VOLUMES IN MILLILI 
TRES" ! "PRESSURES IH CM OF ME 
R C U R Y " ! " " 

40 V=INTCRND*50+50>*100 
50 V1 = INTCRND*50 + 50.^*100 
60 P=INT<RND*150+150> 
70 P1=INTCRHB*150+150> 
8 C=C+1 

30 ON C SOTO iOO, 140, ISO, 220 
, 260 
100 V=0 

110 PRINT "SOLVE FDR V WHEN" 
120 ei=Pi*Vl/p 
130 GOTO 270 
140 V1=0 

150 PRINT "WHAT IS VI GIVEN" 
160 iai'=P*V-^Pl 
170 GOTO 270 
180 P=0 

1|0^PRINT "ClRLCULflTE THE VftL 

200 Q1=P1*V1/'V 

210 GOTO 270 

220 P1=0 

230 PRINT "SOLVE FDR Tl GIVE 

240 Qi=P*V./Vl 

250 GOTO 270 

260 STOP 

270 PRINT " V VI P 

280 PRINT V; TflB(;7> ? VI J TRB'.: 14 
> ; PJ TflB<21 >5 PI 

"VDUR ANSWER ="!0 
3y0 PRINT "CORRECT VftLUE 
W 1 

3 10 PRINT 
J20 GOTO 40 
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>RUN 

DRILL DH BDYLE'S LRU 

VOLUMES IN MILlsIUJTBES,,,^,. 
PRESSURES IN CM DF MERCURY 

SOLVE FOR V WHEN 

V VI P PI 

7400 251 249 

YOUR ANSWER =7350 
CORRECT VALUE = 7341.035857 

WHAT IS VI GIVEN 

V VI P PI 
9400 232 215 

YOUR ANSWER =10060 

CORRECT VALUE = 10143.25581 

CALCULATE THE VALUE OF T 

V VI P PI 
8800 6700 193 

YOUR ANSWER =144 

CORRECT VALUE = 146.9431818 

SOLVE FOR Ti GIVEN 

V VI P PI 
6300 9700 265 

YOUR ANSWER =170 

CORRECT VALUE = 172.1134021 



** BONE ** 
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Photoelectric Emissions 



When light of a short wavelength falls on a metal surface, electrons are 
emitted from the metal. According to the description of this phenomenon by 
Einstein, there is a maximum wavelength for every metal above which no 
electrons are emitted. This is called the critical wavelength of the metal. 

This program simulates a laboratory experiment in which a metal is placed 
in a vacuum and bombarded with soft X-rays. The number of electrons emit- 
ted is collected and measured with a microammeter. The program simulates 
three trials at each of nine wave lenghts. 

After each set of experimental data, the program asks if you would like an- 
other run at a higher light intensity. The reason for doing this is that some- 
times at low light intensities, not enough electrons have been emitted for 
meaningful measurements. 

You can increase the precision of the experiment by increasing the number 
of wavelengths at which it is run. This value can be changed in Statement 60; 
note that the variable L is divided into 1000 in order to express the wave- 
length in Angstroms. One Angstrom (A) equals 10"^ centimeters or 10 
microns. 

Here are the coefficients for several metals: 



Silver .308 

Bismuth .338 

Cadmium .318 

Lead .340 

Platinum .385 



It is a rare physics laboratory in a high school or college today that has the 
experimental apparatus to run this experiment, yet with a small computer the 
equipment can be simulated. There are many other things with which you 
would not normally be able to experiment that can be simulated with a com- 
puter. Things such as a nuclear power plant, a malaria epidemic, an urban 
mass transit system, and a bicycle factory. 

Can you write a simulation for a real world system? You will find sections 
on simulations in the books Computers in Mathematics and Computers in Sci- 
ence and Social Studies. These, along with articles in 99'er and Creative 
Computing magazines, might be of some help in writing a simulation of your 
own. 
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iO PRINT '^PROGRfl'M SIMULflTES^ 
METAL" : "BDtlBflRDMENT WXTH SOF 
T X - R ft Y S " ! " " 

20 INPUT "CnEFFICEENT OF MET 

hiO K = INTC1+2*RND) 

40 PRINT "OUTPUT IN MICRORMP 

E R E S " s " " 

Bo PRINT "IJftVE, TRIAL TRIRL 
TRIAL"! "LENGTH 1 2 

60 FbR"L=.42 TO .25 STEP -.0 

70 M=INT< 1000/L) 
80 PRINT M? „ ^ 
90 FDR J=l TO 3 
100 IF L>V THEN 130 

110 i«SQR< int<:25*rnd:o 

120 GOTO 140 

130 I=SQRCK*K*100h-INT<35*RND 

ilio N^INTClOsI + .S-J/lO 

.50 PRINT TABC J*7'-» J N; 

,60 NEXT J 

:,70 PRINT 

,80 NEXT L 

:,90 PRINT " ":"INCREfiSE LIGH 
■ INTENSITV" s "<V OR N>"5 

2 00 INPUT AS 

210 IF AS="N" THEN 270 

220 PRINT "BY WHAT FACTOR <1 

TO 1 > " J 

230 INPUT F 

240 K=K*F 

250 PRINT 

260 GOTO 50 

270 STOP 
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>RUM 

PROGRAM SIMULlhlTES METAL 
BOMBARDMENT WITH SOFT X- 

CDEFFICIENT OF METAL = 34 
OUTPUT IN MICROAMPERES 



RAYS 



WAVE TRIAL 
LENGTH 1 



2380 
2500 
26 31 
2777 
29 41 
3125 
§333 
3571 
3846 



20 . 
20 . 
20. 
20. 
2. .3 



TRIAL 
2 

20 . 4 
20. 5 
20. 6 
20 . 4 

1 - y 



1 

4, 
3 , 
4 



7 
4 



TRIAL 



?9- 3 

i?5^ 

4 . 4 
4.6 

3 . 5 



BY WHAT FACTOR 



INTENS I TV 
< 1 TD IQi 



1 



'aIAVE,, TRI AL 
LENGTH 1 



2380 
2500 

294 1 
3125 
3333 
357 1 
3846 



200 
4 



4 

2 

3 , 
3 



TRIAL TRIAL 



2 00 . 1 

200. 1 
200. 1 

3 . 6 
3. 5 
3 . 3 

3 . 7 

4 . 5 



20 
20 
200 
20 

2. 2 

3. 9 
4 . 4 
2. 4 
3 . 7 



INCREASE LIGHT 
<Y OR N>? N 

** DONE ** 



INTENSITY 
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Mutation of Moths 



This program is a simulation of the growth of a colony of pepper moths. 
The program allows a genetic mutation to be introduced in some year be- 
tween 1 and 30. The mutation can favor either dark or Ught colored moths. 

The program as it appears starts with a total colony of light moths; you 
could add an initial group of dark colored moths in Statement 35 as PI. 

The sample run shows a mutation which favors dark moths occurring in 
Year 3. The mutation affects about 2% of the light moths each year and 
causes them to become dark. The program displays the number of moths of 
each color over a 30-year period. 

Obviously, because of the long time period involved it would be difficult to 
carry out this experiment in school, so the computer again is of real benefit. 



i P 
MUTfil 
THS " 
20 I 
1-10 
30 P 
40 Z 
50 P 
CHRN 



R " . 
€.0 I 
70 L 
8 D 
■9 P 

□ R 
100 
1 10 
HS 
120 
130 
140 
150 
160 
170 
130 
1 90 
200 
210 
220 
> 5 PO 
230 



R 1 NT 
TION 

■ n II 

N P IJ T 

> = " ! 

0=10 
= P0 
R INT 
6 E " ! 

NPUT 
-1 
= 2 

Rl NT 

DRRK 

INPU 

PR IN 

LIGH 

FOR 

IF T 

Pl=0 

GQTD 

IF E 

P0 = I 
IF P 
Pl=2 
P0 = 
PRIN 



"PRDGRflM SIMULATES 
s " DF ft COLONV OF MO 



"RATE 

M 

000 



□F MUTATION < 



"IN WHAT YEAR DDES 
IN ENVIRONMENT OCCU 

X 



"CHANGE FAVOR 
" : "MOTHS <L OR 
T ES 

T " " : "YEAR 
T MOTHS" 
T=l TD 30 
>=X THEN 160 



LIGHT 
D) " ? 



DARK MOT 



220 

SO'-D" THEN 140 
NT <P 1 + . 01 *M*PO + . 5> 
NTtZ-Pi+. 5> 
KZ THEN 220 



T T5 TAB<:8> ; PI ; TABC 19 
NEXT T 
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IN ENVIRONMENT DCCUr" 3 



VEflR LARK MDTHS 
A 

3 2 00 

4 396 

5 588 
b 776 

7 960 

8 114 1 

10 1492 

11 1662 

12 1S29 

13 1992 

14 2152 
2309 
2463 
26 1 4 
2762 
2907 

20 3049 

2 1 3 18 8 

22 3324 

23 3458 

24 3589 

25 3717 
2b 3843 

28 4087 

30 4321 



LIGHT, MOTHS 
1 
1 0000 
9800 
9604 
9412 
9224 
9040 
8859 
8682 
8508 
8338 
8171 
8008 
7848 
7691 
7537 
7386 
7238 
7093 
6951 
6812 
6676 
6542 
64 1 1 
6283 
6157 
60 3 4 
5913 
5795 
5679 



The so-called "killer bees" that were introduced into South America in the 
mid 70's were supposed to help honey production because they were much 
more energetic than the lazy honey bees in Brazil. The idea was that they 
would mate with the existing honey bees and produce a more productive 
strain. However, they went on a rampage killing people and terrorizing the 
country. Then they started to migrate north. U.S. agricultural officials be- 
came alarmed that they would invade this country and bring death and 
destruction. 

Over the years, the killer bees have started to mate with the more docile 
South American honey bees, but only at the rate of 3% per year. Assuming 
they pose a danger to the U.S. until their numbers are reduced to 15% of 
their original quantity, how many years will it take for this to occur? The 
present migration patterns will bring them to the U.S. by 1989; will they still 
be dangerous (assume that the year of introduction was 1974)? 
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Projectile Motion 



The path followed by a projectile is called its trajectory. The trajectory is 
affected to a large extent by air resistance, which makes an exact analysis of 
the motion extremely complex. We shall, however, neglect the effects of air 
resistance and assume the motion takes place in empty space. 

In the general case of projectile motion, the body (bullet, rocket, mortar, 
etc.) is given an initial velocity at some angle above (or below) the 
horizontal. 

If Vq represents the initial velocity (muzzle velocity), the horizontal and 
vertical components are: 

V^^ = V^cos0, V,y = V^sin^ 

Since we are neglecting air resistance, the horizontal velocity component 
remains constant throughout the motion. At any time, it is: 
= = Vq COS0 = constant (1) 

The vertical part of the motion is one of constant downward acceleration 
due to gravity. It is the same as for a body projected straight upward with an 
initial velocity Vq sin 6. At a time "t" after the start, the vertical velocity is: 

Vy = VQy.gt = VQSin^.gt (2) 
where "g" is the acceleration due to gravity. 
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The horizontal distance is given by: 

^ = ^oxt = (V^s^)t (3) 
and the vertical distance by: 

y=v^yt-i/2gt^ 

= (VQsin^)t- l/2gt2 (4) 

The time for the projectile to return to its initial elevation is found from 
Equation (4) by setting y = 0. This gives 

2 V.sin e 
t = (5) 

The horizontal distance when the projectile returns to its initial elevation is 
called the horizontal range. "R." Introducing the time to reach the point in 
Equation (3), we find: 

2 V^^sin e cos 6 

R=— (6) 

g ^ ^ 

Since 2sin 9 cos. 6 = sin 2$, Equation 6 becomes: 

sin 26 

R = — (7) 

g ^ ^ 

The horizontal range is thus proportional to the square of the initial veloc- 
ity for a given angle of elevation. Since the maximum value of sin 2^ is 1, the 
maximum horizontal range, Rmax is VoVg. But if sin 2^ = 1, then 26 = 90" 
and 6 = 45°. Hence the maximum horizontal range, in the absence of air 
resistance, is attained with angle of elevation of 45°. 

From the standpoint of gunnery, what one usually wishes to know is what 
the angle of elevation should be for a given muzzle velocity Vq in order to hit 
a target whose position is known. Assuming the target and gun are at the 
same elevation and the target is at a distance R, Equation (7) may be solved 
for 0. 



e = 1/2 sin 
= 1/2 sin 




(8) 



Provided R is less than the maximum range, this equation has two solu- 
tions for values of 6 between 0° and 90°. Either of the angles gives the same 
range. Of course, time of flight and maximum height reached are both greater 
for the high angle trajectory. 
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For example, say the maximum range of our gun is 10,000 yards and the 
target is at 5,900 yards: 



Try the computer game Gunner which appears below. Use trial and error 
to try and destroy the target (see sample run). You get five chances per target. 
How many shots did it take to destroy all five targets? Did you ever fail to de- 
stroy a target in five trials? 

From the discussion above, you should realize that 45 degrees of elevation 
provides maximum range with values over or under 45° providing less range. 

The maximum range of the gun will vary between 20,000 and 60,000 yards 
and the burst radius of a shell is 100 yards. 

You can also determine the correct firing angle from sine tables, or a slide 
rule, a scientific calculator or a computer. You should be able to destroy ev- 
ery target with just one shot. What happens when the target is very close? 
Can you always use whole angles? 

Now write a computer program to accept the maximum range of your gun 
and the range to the target and then calculate the correct firing angle. You 
will have to solve two problems to write such a program: 

1. The Basic language does not have an ARCSIN function. However, the 
following formula may help. 



e = 1/2 sin" 



-1 



5,900 
10,000 



= 1/2 36° 

= 18°, or 90°- 18° = 72* 




2.Y0U must convert from radians to degrees. 



125 



iO PRJMT "YOU ARE COMMfiNDINS 
."^..^UN" ! "CREW. HIT WITHIN 10 
§o^ggS§ni;§E THE TRRGET":: 
20 RfiHDDMIZE 

30 R= I NT < 40000*RND+20000 > 
40 PRINT "MAX RANGE "J R J" YDS " 
50 Z = 
60 S1=0 

?Q GOTO 250 

PRIWT.'iMIN = 1 DEGREE" 
110 GQTD 260 

JjgO PR IHT^L'MflX = 89 DEGREES" 
130 GDTD 260 

iag„gRINT "OVER BY" ; flBS<:E> I " 
f n R D S 

ISO GDTD 260 

"?SrDS" "SHDRT BY"JflBS<ES? 
170 GDTD 260 

180 PRINT "TARGET DESTROYED! 



II 



190 PR INT^I? "ROUNDS USED":: 

210 IF~274 THEM 410 
220 2=Z+1 

230 PRINT "THE FDRUfiRD DBSER 

VER^HflS"^;;SIGHTED MORE ENERy 

240 GDTD 70* 

250 PRINT "TARGET AT"ST5"YflR 
260 INPUT "GUN ELEVATION ~" t 

129 JE 120 
i§9 JEF„B<1 THEN 100 
290 S=S+i 

300 IF S<6 THEN 340 
310 PRINT TAB< 12>? "BDDM ! " s TA 
?<5) 5 1BD0M! "!TAB<15>"^BaDM " 
! BDOh ! " : : 

320^,PRINT "ENEtlY GOT YDU FIR 

330 GDTD 450 
340 B2=.035*B 

370 E = INT<>«!> 
400 GDTD 140 

|10^PRINT "TOTAL ROUNDS US 

420~IF*S1>18 THEN 450 

liR ESiHT "NICE SHDDTIN6! " 

440 GDTD 460 

^^pO^PRINT "GO BACK TO FT. SI 
G CDUR^^e"" "REf^'^ESHER TRAIN IN 
460 PRINT 

470 INPUT "TRY AGAIN <Y,N> " 
480 IF ZS="Y" THEN 30 
E^CnKp?'-"^ RETURN TO BBS 
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YOU ARE CaMMflNDING„fi gUN^^ 
CREW. HIT WITHIN 100 YflRIiS 
□F THE TRRGET 

MAX RRNGE 41167 YBg^^ 
TARGET AT 17004_YflRDS 
GUN EL,EVflTIDH^ = 77 
□VER BY 776 YARDS 
SUN ELEVflTIDH =78^ 
SHORT BY 534 YARDS, 
GUN ELEVATION =77.6 
TARGET DESTROYED! 
3 ROUNDS USED 

THE FORWARD OBSERVER HAS , 
SIGHTED MORE ENEMY ACTIVITY 

TARGET AT 21714 YARDS 
GUN ELEVATION =73^ 
OVER BY 1074 YARDS 
GUN ELEVATION =74 
SHORT BY 140 YARDS ^ 
GUN ELEVATION =73.9 
TARGET DESTROYED! 
3 ROUNDS USED 

THE FORWARD OBSERVER HAS 
SIGHTED MORE ENEMY ACTIVITY 

TARGET AT 16942 YARDS 
SUN ELEVATION =77.6 
TARGET DESTROYED! 
1 ROUNDS USED 

THE FORWARD OBSERVER HAS 
SrlHTED MORE ENEMY ACTIVITY 

TARGET AT 12531 YARDS 
GUN ELEVATION =82^^ 
SHORT BY 1487 YARDS 
GUN ELEVATION =81 
SHORT BY 106 YARDS 
GUN ELEVATION =80.9 
TARGET DESTROYED! 
3 ROUNDS USED 

THE FORWARD OBSERVER HAS 
SIGHTED MORE ENEMY ACTIVITY 

TARGET AT 22631 YARDS 
SUN ELEVATION =73. 1 
TARGET DESTROYED! 
1 ROUNDS USED 



TOTAL ROUNDS USED = 11 
NICE SHOOTING! 

TRY AGAIN CY-N^ N 

OK. RETURN TO BASE CAMP. 

** DONE ** 
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9 

Potpourri 



Here are five programs that didn't seem to fit anywhere else. The first and 
the last are games, although you may come to think of the lunar landing 
simulation as a game also. One is a nifty simulation of smog, and the other 
calculates depreciation by three different methods. 
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Number Guessing Game 



Here is a computer program that plays the popular number guessing game. 
In it the computer picks a secret number between 1 and 100. You attempt to 
guess that number in as few tries as possible. 

There are many ways to go about guessing the secret number. Let some of 
your friends play this game and see what approaches they use to find the se- 
cret number. One approach is to start with a guess of 10. If this is too low, in- 
crease each guess by 10 until the computer says that a guess is too high. 
When this point is reached, start from the previous guess and increase each 
guess by 1. This is the method used to solve some of the problems in the chap- 
ter on convergence. 

Another approach is to try to bracket the number between upper and lower 
limits and reduce the limits by steps until the number is finally found. Two of 
the convergence programs used this approach. 

Is one of these the best way? Well, these methods are not bad, particularly 
compared to starting with 1 and simply counting to 100 until the solution is 
found. But there is a better way. It is known as binary search. 

This technique involves dividing the search domain, in this case 1 to 100, in 
half, and then in half again, and so on until the secret number is found. Play 
the game many times using different approaches. In the long run you should 
find that the binary search approach is the most efficient. 

In Line 210, the program contains the statement that "you should not need 
more than 7 guesses." Why? If the secret number was between 1 and 128, 
what is the maximum number of guesses that would be necessary to find it? 
What if the number range were 1 to 130; then what would be the maximum 
number of guesses? 

Revise the program to choose a secret number between 1 and 10,000. Now 
the upper limit is 100 times the 1 to 100 game here which requires a maxi- 
mum of seven guesses to find the secret number; how many guesses will now 
be required? 

Can you write a program in which the roles of the computer and player are 
reversed? In other words, the computer will try to guess your secret number 
between 1 and some upper limit. After each guess, you enter L for low, H for 
high, or C for correct. Can you write this program so the computer can tell if 
you are cheating, i.e., giving it inconsistent clues? 
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10 RflMIiDMIZE 

20 PRINT "I HAVE CHOSEN fl SE 
CRET" :" INTEGER BETWEEN 1 AND 

10Ci":"TRY TO GUESS UNFIT IT 
IS."! "I WILL GIVE YOU CLUES" 

3 PRINT 

4 C = 

5 N=INT<RND*100>+1 

6 C = C + 1 

70 INPUT "YOUR GUESS ="!G 

SO IF N=6 THEN 140 

90 IF G>N THEN 120 

100 PRINT "TDD LDW" 

110 6DTD 60 

120 PRINT "TDD HIGH" 

130 GDTD 60 

140 PRINT "CORRECT IN"!!C;"6U 
ESSES" 

150 IF C>3 THEN 180 

160 PRINT "YOU WERE LUCKY!" 

170 GOTO 220 

180 IF OF THEN 210 

1 90 PR I NT "GOOD JOB . " 

200 GDTD 22 

210 PRINT "BUT YGU SHDULBN'T 
,,NEED MDRE":"THftN 7 GUESSES. 

220 PRINT "PLRY PlGFlIN <Y DR 
N> " ; 

230 INPUT flS 
240 PRINT 

250 IF flS="Y" THEN 20 
260 STOP 



■RUN 

I HfiVE CHOSEN fl SECRET 
INTEGER BETWEEN 1 AND 100 
TRY TO GUESS WHAT IT IS. 
I WILL GIVE YOU CLUES 

YOUR GUESS- =50 
TOO LOW 

YOUR GUESS =75 
TOO LOW 

YOUR GUESS =87 
TOO LOW 

YOUR GUESS =93 
TOO LOW 

YOUR GUESS =97 
CORRECT IN 5 GUESSES 
G D n JOB. 

PLAY AGAIN <Y OR N)? N 



* * DONE * * 
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Depreciation— Three Methods 

This program shows how a piece of capital equipment depreciates accord- 
ing to three commonly used methods of depreciation: straight line, sum of the 
year digits, and double declining. 

The program asks for the original cost of the item, its expected life in years 
(the period of time over which it is to be depreciated), and its expected scrap 
(or sale) value at the end of that time. A table showing the annual depreci- 
ation for each of the three methods is then displayed. 



ETHDdI";" "^^"^^^^ THREE M 

20 INPUT " □ R I G I N ft L COST = " : f 
30 INPUT ''LIFe'dF ITEri c:VEflR 

^8 IrVnT "^'^^^P VftLUE =":S 

70 V=C-S 
80 D1=V/L 

10 FDR X=l TO L 

IS f^^^VP X5TfiB<4>S 

- f u y = D 1 

180 GDSUB 2 80 

138 TflBC12>; 

230 PR INT TftB < 2 i > 5 

lis 

280 Q=INTCQ*100+.5>/100 

Ii8 fliS?*«°..T"^" ^'3° 

138 TMH 350 

118 SIIBSn*^' 



CI ■ 

? 
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PRQGRRM CftLCULflTES DEPRE- 
ClfiTION BY THREE METHODS 

□RieiNflL COST =3000 
LIFE OF ITEM <YEfiRS> =8 
SCRAP VALUE =500 



STRAIGHT 
YEAR LINE 



1 

2 
3 
4 
5 
6 
7 
3 



937 
937 
937 
937 
9 3 7. . 
9 3 7.5 
937 . 5 
937 ~ 



5 
5 
5 
5 
5 



SUM OF DOUBLE 

DIGITS DECLINIMG 

1666.67 2000 

14 58.33 1500 

1250 1125 

1041.67 843.75 

833 .33 632.81 

625 474.61 

416.67 355.96 

2 8.33 266.97 




Smog Simulation 



This program is an adaptation of the smog model originally written by 
Herbert Peckham. The model assumes that vehicular traffic is the sole pro- 
ducer of smog, a somewhat poor assumption. It also assumes that most auto- 
mobile traffic occurs during the daylight hours and that traffic volume is very 
low (actually, zero) at night. The smog generated by the cars is dissipated by 
atmospheric conditions which vary depending upon sunlight, temperature, 
and weather. All of these conditions may be specified by the user. 

The model could be improved significantly by taking into account other 
sources of smog, by varying vehicular traffic according to the hour of the day, 
and by allowing daily variation of weather factors. Nevertheless, even in this 
rudimentary form it is interesting and instructive. 

A plot of the smog level is produced in Statement 370. Under some con- 
ditions, the smog level reaches a value that cannot be plotted because it is 
greater than the width of the screen (and printer). For runs with conditions 
like this, you might want to delete the plotting routine. A more elegant solu- 
tion would be to estimate the maximum value of the smog level from the in- 
put factors and calculate an appropriate plotting multiplication factor. 
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10 PRINT "SMOG CITY SHIULlRTI 
N " s " " 

O T 

:|o PRINT "CITY CARS PER^RO 
flD MILE": "LOS flNSELES 200" 
""DETROIT iaO":"TULSR 

25" ! "SMDS CITY " 5 
40 INPUT C 
s n P R T W T 

60 PRINT " SMDS GENERRTID 

N CDNSTRHT " ! " I 950 RUTD .00 
5":"li7-5 flUTD .001"8"BUS 
.010" : "SMDS CITY " 5 
70 INPUT LI 
on PRINT 

"?0 PRINT " DAYTIME SMOG BRERK 
DOWN RRTE IN PERCENT PER H 

n ij p i» 

100 PRINT/' CLEAR., „ ,^02" 

: " CLQUDY .01 SMDS 

C I T V " 5 

110 INPUT Rl 
120 PRINT 

130 PRINT "NIGHTTIME SMOG BR 
ERKDDWN" : "RATE„C-/. PER HDUR>;; 
140 PRINT " HOT .05" 
! " COOL . 10" : " SMDG 

CITY"; 

150INPUTR2 

170 Print "smds dispersion < 

V PER H □ U P > " 

i'SO PR INT ^"HIGH WIND„ftND^ Rfll 
N 1.00": "NO WIND AND DRV 
0.01"! "SMOG CITY 
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190 INPUT R3 

200 PRINT " ":"HDUR SMOG LEV 
EL PLOT" 
210 R«R1 

lis h'WAi]un'im,THEu 27 

25 Kl=0 

260 R=R2 

270 S=S+K1#C*10-R*S-R3*S 

28 IF S>0 THEN 290 

290 T=sT + l 

300 T3=INT<T-^12) 

310 T2=T-T3#12+1 

320 X = I NT <; 1 000*S+ . 5 J 1 000 

330 IF T3.''2=INTCT3/'2) THEN 39 

340 PRINT T2; "PM" 5 TRB<7> ; }<!! 

350 IF X<14 THEN 370 

36 X = 1 4 

370 PRINT TAB<13+X>; "*" 

380 GOTO 210 

39 PRINT T2; "AM"? TAB(;7> SJ<J 

400 GOTO 370 
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> R U N 

SriDG CITY SIMULRTIDN 
SMQG CITY ? 100 

iSr ""^^ :§?4 

SnOG CITY ? .00 1 

?rr^fig!Ef1!?°^Ei'^^g5g°'^^ ^^^^ 

C L E fil R .02 
CLOUDY :oi 
SMOG CITY? loi 

SMOG CITY? :iO 

r3HDirg?T9'^° ^•'^ ?°:8I 

HDLIR^SMDG LEVEL PLOT 

10 AM 

11 AM 

12 fiM 
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Lunar Lander Simulation 



This program is one of the most popular computer simulations around. It 
is available in many versions; this one is adapted from the original program 
written in 1969. 

The program represents an exact simulation of an Apollo lunar landing 
module during the final descent. This portion of the descent would normally 
be controlled by the on-board computer backed up by another computer in 
the lunar orbiter, and still another computer on Earth. However, to exercise 
your knowledge of physics and to make an interesting game, we will assume 
that all three computers have had a simultaneous malfunction. Hence, it is up 
to you to land the spacecraft safely. 

To make a soft landing, you may change the burn rate of the retro rockets 
every ten seconds. You have a choice of not firing at all (burn rate of 0) or of 
firing at a fuel rate of between 8 and 200 pounds per second. Engine ignition 
occurs at 8 pounds, hence values between 1 and 7 pounds are not possible. 
You have 16,500 pounds of fuel. This is 500 pounds more than an actual 
LEM has, which will give you a little margin for error. When you get pro- 
ficient, change Statement 80 to N = 16000 to simulate the real thing more 
closely. The capsule weight is 33,000 pounds. 

Not that this is the way to come in, but if you did not fire the rockets at all, 
the estimated time for a free fall descent is 120 seconds to impact (and a huge 
splat). 

Good luck! 
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10 PRINT "COMPUTER MfiLFUNCTI 
□ N.! " : "flSSUME MflNUftL CONTRDL. 
20 PRINT "CAPSULE WEIGHT 32, 
?500 LBS f, ; ?''''fiII-flBLE FUEL 16 

3 print' "SET RETRD ROCKET B 
yRN„R9Ii"-,:'EVERY 10 SECONDS 

IB i ^si'^e^o^.f^^^B^'^h^f ii^wf 

'Ru^E^^^6feN"ff^ ^'^^^ "^^^ 

5 ft = 1 2 

6 V = 1 

70 M=33000 
SO N = 1 6 5 O 
?Q . 00 1 
100 2=1.8 

110 P=INT t:3600*V + . 5> 
120 IF P>=1 THEN 140 
130^P=INTt:r:3600*V + .5:-*1000>/ 

140 IF fl<10 THEN 170 
150 p=INT(:fi+.5> 
IbO GOTO 180 

|l8.?RSAr^^8??^§ti!h^^J?gg, 

^g8'?g^6^ TBB<24>, 

lis IE }<>200 THEN 250 
220 IF K>7 THEN 270 

250^PRINT "BftD VlFILUE--flGFlIN 

260 GOTO 200 

270 T=10 

280 IF M-N<.001 THEN 390 

U% l^J^<■''^^ THEN 110 

330 GQSUB 770 

3 40 IF I<=0 THEN 6:^'0 

350 IF V<=0 THEN 370 

3b0 IF J<0 THEN 690 

370 SOSUB 570 

3 B GOTO 280 

f?O^PRINT "OUT OF FUEL flT"5L 

420 L=L+S 
430 W=3600*V 

440 PRINT "ON MOON HT"5INT<L 

? W f??Ulf1PfiCT" : "VELOCITY =" 

5 I N T 1 * W + . 5 > X' 1 ; " M P H . " 

450 IF IJ>1.2 THEN 480 

470 STDp"^ "PERFECT LflWBING!" 

480 IF IJ>10 THEN 510 

490 PRINT "ODDDMPHi THfiT'S P 

□UGH ON": "THE SHOCK RBSDRBER 

500 STOP 
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510 IF U>60 THEN 540 

520 PRINT "SEVERE CRAFT DflMFI 

iSE. VnU ' RE " s " STRflHDED UNTIL 

HELP COMES" 

530 STOP 

540 PRINT "SORRYj BUT NO SUR 
V I V D R S ; " : Y □ U BLEW IT. IN F fl 
CT, YOU " : "BLflSTED R NEW LUNR 
R C P fl T E R " 

550 'PRINT I NT < W* . 273 ) 3 " FEET 

DEEP ! " 

560 STOP 

570 L=L+S 

5S0 T=T-S 

590 t1=»M-S*K 

60 fi=I 

6 10 V=J 

620 RETURN 

630 IF S<.005 THEN 430 

6 40 D = V + SQRt:V*V + 2*fl*<G-Z*K/'M 

650 S = 2*ft.''D 

660 GQSUB 770 

670 GDSUB 570 

680 GOTO 630 

690 U=< 1-M*G.''<Z*K> 

700 S = M*V/CZ*K*<W + SO.R<lJ*W + V/ 

2 . ■>>) + . 5 

710 GDSUB 770 

720 IF I<=0 THEN 630 

730 GQSUB 570 

740 IF J>0 THEN 280 

750 IF V>0 THEN 690 

760 GOTD 280 

7 7 Q = S * K M 

780 J==V + G*S-Z*(a**:i+Q*<.5 + Q3f-< 
1 ^ 3 + Q * < . 2 5 + Q 5 > > > 5 
790 I=R-6*S*S'-'2-V»S + Z#S*i3*<:. 
5 + Q*":i/-6 + Q*< l''12 + Q/20?:O 
900 RETURN 



COMPUTER MflLFUNCTIQN ! „^ 
ASSUME MflNUflL CONTROL OF 
LUNAR CAPSULE... 

CAPSULE WEIGHT 32^500 LBS. 
ftVfllLRBLE FUEL 16,500 LBS. 

SET RETRO ROCKET BURN RflTE^^ 

EVERY 10 SECDNDS TD^O^CZERD? 

□R VALUE BETWEEN 8 AND 200 
LBS PER SEC 

SEC MILES MPH FUEL BURN 

120 3600 16500 ? Q 

10 110 3636 16500 ? Q 

20 100 3672 16500 ? 

30 90 3708 16500 ? 

40 79 3744 16500 ?0 

50 6 9 3780 165 00 7 0^^ 

60 58 3816 16500 ? 200 
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70 48 344? 14500 ? 200 

80 39 3051 1250 0? 20 

?9r> 31. 2624 1050 ? 

100 24 2660 10500 ? 99 

}iS \7 2454 9510 ? 99 

120 10^^ 2238 8520 ? 99 
130 4.321 

2013 7530 '> 99 

BLflSTED^fl NEW LUNAR CRRTER 
506 FEET DEEP ! 

** DONE ** 



COMPUTER MfiLFUNCTIDN ! 
ASSUME, MANUPL CONTRDL DF 
L U N fi R CAPSULE... 



CAPSULE WEISHT 
AVAILABLE FUEL 



32, 500 
16> 500 



LBS 
LBS 



SET RETRD ROCKET BURN RATE 

EVERY 10 SECONDS TO <ZERD 

□R VALUE BETWEEN 8 AND 200 
LBS PER SEC 



SEC MILES MPH 



120 
1 10 
100 
90 
79 
69 
58 
48 
37 
28 
20 





10 
20 
30 
40 
50 
60 
70 
80 
9 
1 00 
1 1 
120 
130 
140 
150 
160 

1 7 
180 
190 

2 00 
210 
2 20 
2 30 

□N MDQN 



1 3 
8 . 



3600 
3636 
3672 
3708 
3744 
3780 
3816 
3852 
3483 
3087 
2660 
2197 



. 07 

.125 

. 796527 
766 212 



FUEL 

16500 
16500 
16500 
16500 
16 5 
16500 
16500 
16 5 
14500 
12500 
10500 
8500 



2500 
1500 
1000 

560 
460 

2 380 

3 290 
1 190 . . 
3i«SES', IMPACT 

□PPDMPH! THRT'S ROUGH DN 
THE SHOCK RBSORBERS 



37 
24 . 
1 76 
139 
1 1 7 
091 
056 
023 
AT 



66 

3^ 

9' 
7 
1 
1 

1 

•-> 



9 
•7 

? 

9 
9 
1> 
9 
9 
9 
9 

9 
9 
9 
9 
9 



BURN 

8 
1 

200 
200 
200 
200 
200 
200 
200 
1 00 
50 
19 

10 
8 
9 

1 
10 



VELOCITY = 
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Hammurabi 



Hammurabi is one of the all-time favorite computer games. On the one 
hand, it may be considered a game, but on the other it is an intriguing simula- 
tion of barter and management. 

Hammurabi is your servant as you try to manage the ancient city-state of 
Sumeria. The economy of the city-state revolves around just one thing — the 
annual crop of grain (probably soybeans). 

Each year, you must determine how many bushels of grain you wish to 
feed to your people (you'll quickly discover how much a person needs to sur- 
vive), how much you wish to use as seed in planting crops for the coming 
year, how much you wish to use for the purchase of additional land from 
your neighboring city-state, and how much you wish to put in storage. 

Of course, if you have a bad harvest or if rats overrun your grain storage 
bins, you may have to sell land in order to get enough grain to keep your peo- 
ple from starving, or to plant the land for the coming year. Unfortunately, 
disasters always seem to strike, forcing you to sell land, when the price is at 
an all-time low; but that's not any different from the real world. 

Most people start to play this game with noble ambitions. However, before 
long, they start longing for a plague to trim their growing population. Or they 
deliberately starve some people to keep things in balance (gosh, maybe these 
zero population growth people have something, after all!). 

Over the years, this game more than any other, has spawned a host of look- 
alikes, extensions, and modifications. Indeed, several manufacturers have 
taken my original with no changes whatsoever, put it in a fancy box, and 
charged a handsome price for it. Accept no imitations! Here is the original 
game (with the dialog shortened slightly) for you to run on your computer. 

If you want to experiment with changes, here are some suggestions. In the 
existing game, plagues randomly occur 15% of the time; lower this to 10% or 
5%. People now require a fixed amount of food; vary this amount slightly 
from year to year. Permit the construction of a rat-proof grain bin, but this 
must cost a fair amount. Introduce a mining industry as well as agriculture. 
How about fishing or tourism? Let your imagination run wild. Experiment! 
Have fun! 
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10 CftLL CLEAR 

IR &'iit!T XfiB<:9>J "HflMMURftBI" 
flHL"?! DflVID H. 

SO^ERINT "TRY TD SUCCESSFULL 
SrGDVERNJ: ;'flNCIENT SUMERIfi^F 

50 D1=0 ■ 
60 P1=0 
29 2 = D 
80 P=95 
?9»S=2800 

log y=30oo 

110 E=H-S 

120 Y=3 

130 RsHz-Y 

140 1=5 

ISO Q=l 

160 B=0 

170 PRINT 

180 Z=Z+1 

190 IF 2=11 THEN 1120 
■fgOpPRINT j;HRMt1URflBI: I BEG 

210 print" "YFfiR" 57!" ""n-i'DC 
gPLEl-|K^VES^?^RNS'''5i5.'.8ER^^ 
220 P=P+I 

^60 PRINT "PDPULftTIDN IS NOW 
J PR I NT "YOU OWN" ? fl! ACRES 

Ll°PER^flCPE'''"^'-^^°"''^' 
its ''?JtjI "^'^"'■^ fiTE"5E5 "BUSH 
SOO^gRINT., ;'HflVE" S S; "BSHLS IN 
310 Y=INtc i7+RND*10> 

shels^rcre"'"'^*^" COSTS"; y;"bu 

IfO/fNPUf "BUY HOW MRNY ACRE 

340 IF Q<0 THEN 1100 

350 IF Y*Q<=S THEN 380 

3^9 SS'SUB 1040 

370 eOTD 330 

380 IF Q=0 THEN 43 

390 R = fl + i3 

400 S=S-Y»Q 

410 C=0 

420 GOTO 510 

eI? "SELL HOW MRNY flCR 

f ^9 it <3< THEN 1 100 

450 IF C!<R THEN 480 

1§9 GOSUB 106.3 

470 GOTO 430 

480 fl=fl-Q 

490 S=S+Y*Q 
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500 C=0 

510 PRINT 

520 INPUT "BUSHELS TO FEED P 
E P L E ? " : Q 

530 IF C<0 THEN 1100 

540 REM INDT ENOUGH SRfilN? 

550 IF CK=S THEN 580 

560 GQSUB 1040 

570 GOTO 510 

580 S-S-Q 

590 C=l 

6 00 PRINT 

610 INPUT "PLANT HDW MANY PIC 
PES?" ■ D 

620 ' IF D=0 THEN 760 

630 IF r<0 THEN 1100 

640 REM IWDT ENOUGH LRND? 

650 IF D<:=fl THEN 690 

660 GOSUB 1060 

670 GOTO 610 

680 REM lENOUGH SEED GRAIN? 

690 IF INT':D.^2)< = S THEN 720 

700 GQSUB 1040 

710 GOTO 610 

720 IF D<10*p THEN 750 

730 PRINT "YOU- DO NOT HAVE E 

NDUGH" ! "PEOPLE TO TEND THE C 

ROPS . " 

740 GOTO 610 

750 S = S-INT<D''2> 

760 GOSUB 1080 

770 REM Ifl BOUNTIFUL HARVEST 

780 Y=C 

790 H=D*Y 

300 E=0 

810 GOSUB 1080 

820 IF INT<C.'2><>C.''2 THEN 85 

830 REM IRATS ARE GOING WILD 

140 E=INT<S/C> 

850 iaS-E+H 

860 GOSUB 1080 

870 REM iLET'S HAVE BABIES 

880 I = INT<C*<20*A + S ) /-P^lOO + l 

890 REM IPEDPLE WELL FED? 

900 C=INT(Q/'20> 

910 REM I15>: CHANCE OF PLflGU 

920 Q=INT< 10*<2*RMD-. 3> > 

930 IF P<C THEN 160 

940 D«P-C 

950 REM ! STARVE ENOUGH FDR 

IMPEACHMENT? 

960 IF D>.45#P THEN 1010 

970 Pl = < <2-l >*P1+D*100/'P)>'Z 

980 P=C 

990 Di=Dl+D 
1000 SOTO 170 
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fl.VE^BEEN ■' ! " r MPEfiCHED AND DEC 

*g^9rRRIHI. ."HflMMURftEI t THENK 
BSSCJ^OF^gKSi^i^^^^^ DMLV";S;" 
1050 RETURN 

^9§P PRIHI. "HfiMMIJRi=IBI : THINK 
ORES »■ ^^'-^ » °WN ONLV"?fi;''fl 

1070" RETURN 

1818 ^iit}^^^*'^^^^*^ 
Sii\S^JgT:?.h;gg^s?IEE'^fKUil 

IjjP^GET;;: "YOURSELF flNDTHER ST 
nip §STP 1260 

1120 PRINT "HRMMURflBI REPORT 
?n,XW9T 10":"YEflRi, YOU ST 

fiRVED" J Dl s "PEOPLE. THBT^S"?! 
NT< 100*P1 >/100; "■/. PER YEAR! " 

1140 PRINT "YOU STARTED yiTH 
10 ACRES": "PER PERSON AND E 

NDED WITH" 5 INT< lC*L>"iO! ! s 

1150 IF Pl>3 3 THEN 1020 

1160 IF L<7 THEN 1020 

1170 IF Pl>10 THEN 1230 

11 30 IF L<9 THEN 123 

11 9y IF Pl>3 THEN 1250 

1200 IF L<10 THEN 1250 

1220 GOTO 1260 

1230 PRINT "YOU RULE WITH A 
HEAVY HAND."! "THE PEOPLE HAT 
E^YDUR GUTS": "AND ARE REBELL 
INS. BETTER GET OUT QUICK! 

1240 GOTO 1260 

1 ^ 5 U PRINT "YOUR MANAGEMENT 

tiki ^g^^jDiS??"""""'' ^° 

12 70 END 
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HRMMURnB I 
BY DfiVID H, flHL 

TRY TD SUCCESSFULLY SDVERN 
ANCIENT SUMERIfl FOR 10 YRS 



HRMtlURflBls I BEG TO REPORT, 
YERR 1 , PEOPLE STARVED 
RND 5 IfJERE BORN. 
POPULATION IS NOW 100 
YOU OWN 1000 ACRES PJHICH 
VIELDED 3 BUSHELS PER RCRE. 
RATS ATE 200 BUSHELS. YOU 
HAVE 2800 BSHLS IN STORAGE. 

LAND COSTS 22„BUSHELS-'' ACRE 
BUY HOW MANY ACRESJ^O^ 
SELL HOW MANY ACRES? 

BUSHELS TO FEED PE0PLE72000 

PLANT HOW MANY ACRES?999 

HAMMURABI: I ^BEG^TD^REPDRT , 
YERR 2 s PEOPLE STARVED 
AND 5 WERE BORN. 
pdpulrtiBn IS„NDU„105, ,^ 
YOU OWN 1000 ACRES WHICH 
YIELDED 2 BUSHELS PER ACRE. 
RATS ATE BUSHELS. YOU 
HAVE 2299 BSHLS IN STORAGE. 

LAND COSTS 22„BySHELS/' ACRE 
BUY HQW MftNV^ ACRESJ^O^ 
SELL HOW MANY ACRES? 

BUSHELS TO FEED PE0PLE72010 

PLANT HOW MANY ACRES?999 
HAMMURABI: THINK RSfllN. YDU^, 
HAVE ONLY 289 BSHLS OF GRAIN 
PLANT HOW MANY ACRES?577 

?§Bg^l''?H ^EgP^E^§T§l?i§"' 

S^SdIr4'II!I PPggClE STRUCK! 
HALF THE POPULATION DIED. 
POPULATION IS„NgW^54,, ^ 
YDU OWN 1000 flCRES„WHICH^^ 
YIELDED 2 BUSHELS PER„ACRE. 
RATS ATE BUSHELS . ^YDU^ _^ 
HAVE 1155 BSHLS IN STORAGE. 

LAND COSTS 26„BySHELS/ACRE 
BUY HOW MANY ACRES? 0„, 
SELL HOW MANY ACRES? 50 

BUSHELS TO FEED PE0PLE?1080 

PLANT HOW MANY ACRES?530 
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HflMMURRBI: I_BES TD REPORT, 

SSSS^ERi ?SgR!r^ ^^^"^^'^^ 

Sc.."i''?il'-i„Pfip6UE STRUCK! 

¥?^Lg^S r§uiS§EI ^g^^ScRE. 
RftTS fiTE BUSHELS. YOU 
HAVE 2700 BSHLS IN STORAGE. 

LAND COSTS 25 BUSHELSx-fiCPE 
BUY. HOW MftNY ACRES? 
SELL HOW MANY ACRES? 50 

BUSHELS TO FEED PE0PLE?620 

PLANT HOW MANY ACRES?309 

HSHtl^-'Bf'B I ! I EES TD REPORT, 

YIELDED 4 BUSHELS PER ACRE 
RATS ATE BUSHELS; YOU 
HAVE 4412 BSHLS IN STORAGE. 

LAND COSTS 24 BUSHELS. ''ACRE 
BUY, HqW MANY ACRES? 
SELL HOW MANY ACRES? 

BUSHELS TO FEED PEDPLE?1200 

PLANT HOW MANY flCRES-'S'^g 

HAMMURABI: I BEG TO REPORT, 

S'^nds^Ff^ """"" 
Sc^'?'<?il'-f ^KkfSWE STRUCK! 

i^y, 20Q ACRES WHICH 

^^?b°^?EMH"'g§k§ECi'' 980'- 

HAVE 51S0 BSHLS IN STORAGE. 

LAND COSTS 20 BU.«;HELS''ArPP 
BUY HOW MANY ACRE§? 100 ^ 

BUSHELS TO FEED PEDPLE?720 

PLANT HOW MANY ACRES?359 

HAMMURABI: I BEG TO REPORT, 

YIELDED 4 BUSHELS PER ACRE. 
5§X§ 9IE BUSHELS. YOU 
HAVE 3717 BSHLS IN STORAGE. 

LAND COSTS 23 BUSHELS. -'ftCPF 
lyy. I^RIi! MANY ACRES? 
SELL HG»J MANY ACRES? 
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BUSHELS TO FEED PEnPLE?860 

PLANT HDW MfiNY filCRES?429 

HflMMURflBis I^EEjS^TQ^gEPgRT , 
YEAR 8 7 PEOPLE STARVED 
AND 12 WERE BDRN. 
PDPULATIDN IS NOW 55 
YOU OWN 1000 ACRES„WH1CH ^ 
YIELDED 4 BUSHELS PER ACRE. 
RATS ATE BUSHELS. YOU 
HAVE 4359 BSHLS IN STORAGE. 

LAND COSTS 25 B U S H E L S AC R E 
BUY HOW MANY. ACRES? 0^^ 
SELL HOW MANY ACRES? 50 

BUSHELS TO FEED PEDPLE?1100 

PLANT HDW MfiNY ACRES?549 

HAMMURABI: 1 ^BES^TQ^REPDRT , 
YEAR 9 , PEOPLE STARVED 
AND 17 UERE BDRN. 
A HORRIBLE PL-plGUE^STRyCK ! 
HALF THE POPULATION DIED. 
POPULATION IS NOW 36 
YOU OWN 950 ACRES WHICH ^ 
YIELDED 1 BUSHELS PER ACRE. 
RATS ATE 1053 BUSHELS. YOU 
HAVE 3726 BSHLS Hi STORAGf 



HAMMURABI! I BES TO REPORT, 
YEAR 10, PEOPLE STARVED 
AND 24 WERE BORN. 
POPuTaTION is now 60 
YOU OWN 950 ACRES WHICH „^ 
YIELDED 3 BUSHELS, PER ACRE. 
RATS ATE 1413 BUSHELS. YOU 
HAVE 2491 BSHLS IN STORAGE. 

LAND COSTS 2 3, B U S H E L S z' A C R E 
BUY HOW MANY ACRES?^0^ 
SELL HOW MANY ACRES? 

BUSHELS TO FEED PEDPLE?1200 

PLANT HDW MANY ACRES?599 

HAMMURABI REPORTS THAT IN 10 
YEARS, YOU STARVED 5 
PEOPLE. THfiT'S .47 
■ ' P F R 'y' E R R ' 

YOU STARTED WITH 10 ACRES 
PER PERSON AND ENDED WITH 
12.6 

YOUR 10-YR PERFORMfiNCE„WflS..^. 
FANTASTIC. WHY DON'T YDU^RUN 
FDR GOVERNOR OF NEW JERSEY? 



HAMMURABI: SO LDNG FDR NOW. 
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Magazines referred to in the text include: 

• 99'er, This magazine focuses on Texas Instruments home computers exclu- 
sively. It carries articles and reviews of peripherals, software and other acces- 
sories. Write for current subscription information: 

99'er 
Box 5537 

Eugene, OR 97405. 

• Creative Computing, This is the leading magazine of software and appUca- 
tions for all small computers. It carries articles, tutorials, how-to applications, 
and extensive in-depth evaluations. 

Books referred to in the text include: 

• Computers in Mathematics: A Sourcebook of Ideas, Hundreds of classroom- 
tested ideas for using computers to learn about mathematics. 

• Computers in Science and Social Studies, Scores of simulation programs in 
biology, ecology, physics and management of real world systems. 

All of these books and magazines (except 99'er) are available from Creative 
Computing, Write or call for the current price: 

Creative Computing 

39 E. Hanover Ave. 

Morris Plains, NJ 07950 

(800) 631-8112 

In NJ (201) 540-0445 
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This Ideabook contains dozens of ways to make the most out of your 
computer for solving practical, everyday problems. The 50 ready-to-run 
programs demonstrate scores of different techniques for solving probleriis in 
mathematics, science and business. 

The ten chapters deal v^ith solving problems by formulas and repetitive 
trials, convergence, recursion, compounding, probability, geometry, science, 
simulations, and drill and practice. 

Some of the problems demonstrate the capabilities of the computer; 
others identify its shortcomings. It is important to be familiar with both the 
strengths and weaknesses of your tools so you con recognize the types of 
|obs for which they are suitable. 
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